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Harald Anlauf 


Abstract 

High energy scattering processes of charged particles are accompa¬ 
nied by radiation of hard photons. Emission collinear to the incident 
particles, which leads to a reduction of the effective beam energy, 
and the possibility to directly measure these photons at the HERA 
electron-proton collider provides important physics opportunities. For 
deep inelastic scattering, the measurement of radiative processes ex¬ 
tends the kinematic range accessible to the HERA experiments to 
lower Q^, as well as helps in separating the proton structure functions 
without the need to run at different collider energies. QED corrections 
to these radiative processes are discussed, and the calculation of the 
model-independent leptonic corrections is described in some detail for 
the complete one-loop contributions as well as for the higher order 
leading logarithms. 
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Chapter 1 
Introduction 


In our quest for the understanding of the structure of matter, scattering 
experiments have always played an outstanding role. Important landmarks 
are the classic experiments by Rutherford leading to the discovery of the 
atomic nucleus, the determination of the size of nucleons and atomic nuclei 
from their electromagnetic form factors by Hofstadter and McAllister p|, 
and the establishing of the parton model of the nucleon in deep inelastic 
lepton-proton scattering at SLAG [§]. 

Supported by an increasing amount of data from high precision scat¬ 
tering experiments at high energy electron-positron, proton-antiproton, and 
electron-proton colliders, the electroweak Standard Model (SM) including 
Quantum Chromodynamics (QCD) is currently accepted as the simplest 
quantum held theory describing the observed phenomena of elementary par¬ 
ticle physics. The fundamental particles of the SM, the leptons, quarks, 
gauge helds and the (experimentally not yet established) Higgs are point¬ 
like. The very good agreement of many theoretical predictions and present 
experimental results from LEP, Tevatron and HERA indicates that any sub¬ 
structure of these particles must be smaller than 10“^® m, corresponding to 
a compositeness scale above several TeV iii- 

On the theoretical side, still one of the major challenges is a complete 
understanding of hadronic interactions and hadron structure. The startup 
of the lepton-proton collider HERA at DESY with a center of mass energy 
of ~ 300 GeV (27.5 GeV electrons/positrons on initially 820 GeV, cur¬ 
rently 920 GeV protons) marked the beginning of a new era of deep inelastic 
scattering (DIS) experiments, extending the kinematic domain for structure 
function measurements covered by earlier hxed-target experiments in the 
Bjorken variable x and momentum transfer by several orders of magni¬ 
tude. Besides the inclusive (w.r.t. the hadronic hnal state) electroproduction 
structure functions, HERA also provides increasingly precise data for the 
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study of hadronic final states in DIS, in particular multi-jet events and event 
shapes. 

The large accessible kinematic region at HERA allows the investigation of 
different aspects of strong interaction dynamics. In the perturbative regime 
of QCD, much interest is devoted to the study of the transition from the tra¬ 
ditional Altarelli-Parisi (DGLAP) evolution of structure functions at large 
X to the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 0 , |] behavior at small x, 
and to a determination of the gluon density which is expected to be strongly 
rising for x —>■ 0. Unitarity requires that the rapid rise predicted by the per¬ 
turbative evolution must eventually be damped. There are indications from 
the DIS data at low x that this unitarization manifests itself in a phenomenon 
called geometric scaling [Q, where the cross section is not a function of two 
independent variables x and but rather a function of a single variable 
r = Q^/Qs(x). The function Qs{x) is interpreted as a saturation scale in a 
description of non-linear QCD evolution [^, ^ . 

Of the structure functions of the proton, F 2 {x,Q‘^) and Fl{x,Q‘^), the 
longitudinal one, Fl, is much more difficult to access. The measurement 
of Fl at low X is particularly interesting due to its tight connection to the 
gluon distribution in the proton. There exist several ways to separately 
extract F 2 and F^ from the experimental data. One possibility is to run the 
collider at different —usually lower— center-of-mass energies, which may 
not be desirable from the point of view of other parts of the physics program 
looking for possible hints of physics beyond the SM at the highest accessible 
scales . 

Indirect methods for the determination of Fl can be used at fixed collision 
energy but usually require substantial input from theory or other assump¬ 
tions, and they depend more or less on the modeling of the hadronic final 
state, like extrapolations or QCD fits (see e.g., 1300,113), or the mea¬ 
surement of the azimuthal angle distribution of final state hadrons 
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Krasny et ah suggested a direct method that utilizes radiative events 
with an exclusive hard photon registered (tagged) in a forward photon detec¬ 
tor (PD). Such a device is actually a central part of the luminosity monitoring 
system of the HI and ZEUS experiments [^, ^ that measures the radiative 
elastic process ep —>■ ep + j ; it was further improved during the HERA 
luminosity upgrade . The idea of this method is that emission of photons 
in a direction close to the incident electron corresponds to a reduction in the 
effective beam energy. The effective electron energy for each radiative event 
is determined from the energy of the photon observed in the PD. 

Besides measuring Fl, radiative events extend the accessible kinematic 
range to lower values of Q^. This is important for a closing of the gap between 
the fixed-target experiments and HERA. The potential of this method is sup- 
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ported by preliminary results from the HI collaboration of an analysis at low 
Q‘^ for F 2 @ |2^, |26| (for earlier analyses that did not take into account 

m 


QED radiative corrections see 
determination of was studied in 


]). The feasibility of the corresponding 
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However, with currently analyzed 
data sets it is not yet possible to compete with F^ from extrapolations or 
QCD hts 13 . 


A precise analysis of experimental data requires the inclusion of radiative 
corrections. In the case of the exclusive radiative events, even at HERA en¬ 
ergies only the kinematic region of momentum transfers far below the masses 
of the W and Z bosons has a signihcant cross section. As a consequence only 
the QED subset of the electroweak interaction is relevant, so that we can 
safely restrict ourselves to the calculation of the QED corrections. A minor 
complication is the well known fact that the full theoretical control of the ra¬ 
diative corrections in non-radiative deep inelastic scattering is aggravated by 
the dependence of the corrections on the knowledge of the proton structure 
functions even in kinematic regions which may be difficult to access. How¬ 
ever, this is almost precisely where a measurement of the radiative process 
can help, and why it is important that this measurement is performed. 

From the calculation of radiative corrections to deep inelastic scattering 
it is well known that the model independent QED corrections on the lep¬ 
ton side are most important. For the tagged photon reaction in the HERA 
kinematic regime, the leptonic QED corrections have been discussed at the 
leading logarithmic level |^, ^ and taking into account next-to-leading 
logarithms [^, Q. Since the leading logarithmic corrections lead to 


compact and transparent expressions, we use them for a qualitative study of 
the dependence of the QED radiation effects on the kinematic reconstruction 
method. 

A central part of the present work deals with the complete leptonic O (a) 
corrections^ to this process which were outlined in . We will present here 
all essential details of this calculation. Since it is well known that also higher- 
order corrections are very important in a quantitative description of deep- 
inelastic scattering, we shall describe in some length a systematic method to 
obtain the leading logarithms to the tagged photon process to all orders. 

The organization of this report is as follows. Section 2 introduces our 
notation in the context of deep inelastic scattering to lowest order. We 
recapitulate some basic knowledge of the proton structure functions and their 
connection to QCD. After providing a brief review of the general features of 
QED radiative corrections to deep inelastic scattering, section 3 emphasizes 


^In this work corrections are always counted relative to the leading contribution for a 
given final state, which is already proportional to for radiative deep inelastic scattering. 
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the importance of the exclusive radiative processes that are the main subject 
of this work. In section 4 we calculate the QED radiative corrections to the 
tagged photon process, beginning with the leading logarithmic approximation 
and extending to the complete leptonic corrections. Section 5 considers the 
leading logarithmic contributions at higher orders. We conclude with a short 
summary, recent experimental results on the structure function F 2 obtained 
with the tagged photon method and an outlook on the application of photon 
tagging to e’''e“ collisions. Finally, the appendices explain several details 
of the more involved parts of the analytic calculations and collect several 
formulae that are needed for a numerical implementation. 
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Chapter 2 

Deep Inelastic Scattering 


2.1 Kinematics 

We start by considering the general deep inelastic electron proton scattering 
process 

e{p)+p{P)^l{p')+X{Px), (2.1) 

where e and p represent the incoming electron and proton, and I and X the 
scattered lepton and the hnal hadronic system, while the respective fonr- 
momenta are given in parentheses. Althongh the HERA ep collider is de¬ 
signed to also rnn with positron beams instead of electrons and even col¬ 
lected more integrated Inminosity with the former, we nse the latter as the 
generic denotation of the lepton beam below, bnt we shall explicitly indicate 
the dependence on the lepton charge where needed. 

Eqnation 0) describes the so-called nentral cnrrent (NC) processes, 
where the scattered lepton eqnals the incoming one, I = e, as well as charged- 


e I 



Fignre 2.1: Deep inelastic lepton-proton scattering 
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current (CC) processes, where the outgoing lepton is an electron neutrino 
Vf, for electron beams and an electron anti-neutrino Vf. for positron beams, 
respectively. 

Let us introduce the notation used in this report. We choose coordinates 
in the HERA laboratory system in which the incoming electrons move in the 
positive z-directionj^ whereas the protons move in the negative ^-direction. 
Furthermore, the direction of the a;-axis is chosen such that the momentum 
of the scattered lepton lies in the x-z plane: 

p = (Ee,0,0,Pe), 

P = {Ep, 0, 0, -pp) , 

P' = {E'l, p[ sine, 0,p'i cos 6) . (2.2) 

The HERA ep collider is typically run at an energy of 27.5 GeV for the 
electron beam, while the nominal energy of the proton ring, initially 820 GeV, 
now reaches up to 920 GeV. 

The kinematics of the scattering process (|2.1| ) is described by several 
Lorentz invariants. The square of the total center of mass energy available 
is given by 

s = (p + py = m^ + M^ + Sz^ AEeEp , (2.3) 

where we introduced the related invariant 


S = 2P-p. 


(2.4) 


At high energies, the electron mass m and the proton mass M can be mostly 
neglected, and s zi S. We shall therefore use them almost synonymously 
below. For the HERA beam energies mentioned above, \fS ~ 300 GeV and 
a/ 5 ~ 318 GeV, resp. 

The invariant momentum transfer from the lepton to the hadronic system 
reads 

= -q^ = -{p- p'f ~ AEeE'i sin^ ^ , (2.5) 

assuming that the energy and angle of the outgoing lepton can be measured. 
Further commonly used invariants are the dimensionless Bjorken variables 


X 


y 


Q2 ^ EeE/sin^f 

2P-q~ Ep (Ee - E'l cos2 f) ’ 


P ■ q 
P ■ p 



1 - 


Ei 




^ cos^ 


0 

2 ■ 


( 2 . 6 ) 

(2.7) 


^As a consequence, one should always keep in mind our definition of the forward di¬ 
rection being along the positive z-axis as the initial lepton direction, in contrast to the 
historical decision of the HERA experiments to use the direction of the incoming proton 
beam for that purpose. 






The quantity y corresponds to the relative energy loss of the lepton in the rest 
frame of the incoming hadron. It is easy to show that x and y are restricted 
to the range 

0 <a;, 2 /<l. ( 2 . 8 ) 

Finally, the total mass of the hadronic hnal state X is obtained as 

W‘^ = Ml = {P + qf = M^ + = M^ + (1- x)yS . (2.9) 

X 

Elastic ep scattering corresponds to a; = 1. 

2.2 Determination of kinematic variables 

Because of four-momentum conservation, the kinematic invariants x, y and 
= xyS can in principle be determined in different ways from the measured 
hnal state. Instead of using only the information from the scattered lepton, 
one can also take into account the measured hadronic hnal state. In the 
case of charged current scattering, where the hnal state neutrino escapes 
detection, this is actually necessary. 

In the description of deep inelastic scattering in the language of the naive 
quark-parton model (QPM), the lepton scatters elastically with a quasi-free 
quark within the proton. The two-body hnal state of this subprocess is 
completely determined by two variables, which can be taken to be e.g., the 
energy Eg and the polar angle 6 of the electron in the lab frame, or quantities 
constructed from the hadronic hnal state. 

Assuming local parton-hadron duality, one identihes the energy and di¬ 
rection of the hadron jet with the energy and direction of the scattered 
quark. Thus, denoting the four-momenta of the produced hadrons by ph = 
{Eh;p±^h,Pz,h), one can construct the following quantities:^ 


Sh = '^{Eh+Pz^) , 

h 



( 2 . 10 ) 


the latter being the total transverse momentum of the hadronic hnal state. 
Neglecting masses, one can dehne an inclusive angle 7 of the hadronic system 
which then corresponds to the angle of the scattered quark in the QPM, 

tan — = —— . ( 2 . 11 ) 

2 PT,h ^ ^ 

^Remember that the initial proton moves in the negative z direction. 
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Obviously one can define similar quantities for the scattered electron, 

6 S 

Eg =+ cos 0 ) , pTe = E'eSm.9^ ^ tan - = —( 2 . 12 ) 

2 pT,e 

Starting from these four variables 6 ^, S/^, 7 ], among which only two are 
independent, one derives the following “basic methods”, that use only two 
variables: 


For neutral current events, the electron method (e), which uses E'^ and 
9, see eqs. 0 - 0 . and 

the double angle method (DA), which uses 9 and 7 [^, whereas 


for charged current events only the “hadrons only” method (h), also 
known as Jacquet-Blondel method (JB) |^, is available. 


The Jacquet-Blondel method defines: 


yjB — 


2Eg 


Qjb — 


PT,h 


yjB 


XJB = 


q: 


JB 


yjBS 


(2.13) 


For neutral current processes, there are possibilities and reasons to also con¬ 
sider methods that employ more than two of the measured variables to deter¬ 
mine the kinematics, either analytically or using the full event information 
for kinematic fits. Real detectors are not perfect; they possess finite energy 
and angle resolution, and they usually do not cover the full solid angle, so 
that particles may escape in the beam pipe. Therefore e.g. the energy of the 
incoming electron may not be well known as it might emit a collinear photon 
before the hard collision. 

As an example for an analytic method using three variables we mention 
here the S method, which basically uses (|2.13|) , but replaces 2Ee with the 
help of four-momentum conservation by the combination Eg -|- thereby 
making the determination of the variable y independent of collinear photon 
emission in the initial state, and replaces the hadronic transverse momentum 
PT,h by the leptonic pT,e, 


2 /s = 


Ql = 


Pl,e 


Xy. = 


Q\ 


(2.14) 


Eg + E;,’ 1-2/s 2 /sR 

For a thorough discussion of the essential features of these methods and some 
extensions we refer the reader to []37|, and references cited therein. 
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2.3 The Born cross section for ep eX 


Let us now turn to the model-independent description of the cross section for 
semi-inclusive deep inelastic scattering, where only the hnal state electron is 
measured (electroproduction). For the time being, we shall restrict to the 
case where the momentum transfer is far below the electroweak scale, 
i.e., M|. The Born contribution to the scattering amplitude taking 

into account only one-photon exchange reads: 

Sf, = z{2nYS\Px+p'-P-p)u{P)ru{p)^ {Px\j;^{0)\P) . (2.15) 

where Qe is the charge of the lepton in units of the proton charge {Qe = -l-l 
for the positron), and is the electromagnetic current operator. 

The differential cross section for the scattering of unpolarized particles 
reads {q = p — p')- 


da = 


with the lepton tensor: 


1 fP 


2\ 2 


2 a/A^ 




d^p' 


(27r)32E' 


= i;^[u{p')l^u{p)]*u{pX''u{p) 


spins 

= + Ap^^p\ 

ft u 


r = 9^^ 


TP ~n a aP-9 

^ , p^ =p^ - q^^T 


We normalize the hadron tensor as follows: 

= \Y.Y. (p|j'”'(o)|Px>{Pxvr('))i^’> 

pol X 



(2.16) 


(2.17) 


(2.18) 


Using current conservation, charge conjugation and parity symmetry, we de¬ 
compose this tensor as: 


H^u{P, q) = dTT ( -g^uFi{x, Q^) + P^P^ ^ ^2(0;, Q^) 

~ ~ 2x \ 

= 47r{-~g,,F,{x,Q^) + P^P, — F2ix,Q^)) , (2.19) 


Pu = Py- qi 


P ■ q 


X = 




2P-q 




, Q^ = -q^- 
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The functions Fi{x,Q‘^) and F 2 {x,Q'^) are the (electromagnetic) structure 
functions of the hadron. 

The quantity Xs appearing in the flux factor reads: 


X^ = . ( 2 . 20 ) 

Expressing the final electron momentum via the invariants {y, Q^) yields: 

d^p' 1 S 


{2nf2E', (47r)2yA^'^^'^^ ' 

Inserting into ( |2.1fi|) , we obtain for the differential cross section: 
d^cr S 


dydQ^^ 


Q" 2 As 
27ra^ S 




2(1 -y)- ) F 2 {x, Q^) 


2 1 


2 m^ 


xy^Fi{x,Q^) 


It is convenient to introduce the ratio R, 


R{x,Q^) = 


2, Fi(a;,Q2) _ ^ F2{x,Q^) 


= 1 + 4x 


Ft{x,Q^) V QV 2a;Fi(x,Q2) 

of the longitudinal and transverse structure functions, 

,M 2 

V W 
Ft{x,Q^) = 2xFi{x,Q^). 


Fl{x,Q‘^) = (1+ 4x^^ ) F2 (x, Q^) - 2xFi(x, Q^) 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


(2.24) 


Neglecting the small mass of the electron, and changing the kinematic vari¬ 
ables from (i/, Q^) to {x,y), we hnally obtain for the Born cross section: 


dV 

dxdy 


27ra^ 
Q'^xy 
27ia‘^ ^ 
Q‘^xy _ 


1 + (1 - 1/)^ + ) E2(x, Q'^) - y^Fiix, Q'^) 




2(1-2/)-2xV7^ 




-h 1 4x^-^ 


QV l + R{x,Q^)\ 


F2{x,Q^ 


(2.25) 


The r.h.s. of the cross section (|2.25|) depends on the two proton structure 
functions to be measured, F 2 {x, Q^) and Fl{x, Q^). For hxed center of mass 
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energy, only two of the three invariants x, y and can be varied inde¬ 
pendently. The separate determination of the structure functions at hxed x 
and requires the variation of y and therefore the variation of the center 
of mass energy. For the HERA collider, this means running with different 


beam energies [39 


2.4 The Born cross section at high energies 


The HERA ep collinear was designed to reach center of mass energies at 
the same order of magnitude as the masses of the electroweak gauge bosons, 
TR ~ 3 ■ 10 ^ GeV > Mz,Mw, so that both neutral and charged current 


reactions can be observed in the same experiments. For the neutral current 
process, ep —>■ eA, this means that we have to take into account not only 
photon but also Z-exchange. 

The generalization of the cross section (|2.25|) that takes into account both 
7 - and Z-exchange and their interference can be written as (see e.g.. 


d^ 


cr. 


NC 


dx dy 


27ra^ 

Q'^xy 


{Y+T2{x,Q^)±Y_xT^{x,Q^) 


y^TL{x,Q^)} , 


(2.26) 


with 

Y^{y) = \±{\-yf . (2.27) 

In ( p.26|) we have used the ultrarelativistic limit, S 3> M^, which is 

applicable at HERA energies. Upper and lower signs correspond to electron 
and positron scattering, respectively. 

The structure functions T 2 and Tl are generalizations of the structure 
functions F 2 and Fi^ of the previous section where only photon exchange was 
taken into account. The additional structure function JF 3 is due to the parity 
violating contributions from Z-exchange to the scattering amplitude. 

At HERA energies also the charged current reaction ep —>■ Vf.X has an 
appreciable cross section, which reads: 


dVgc ^ GIS f Ml, Y 
dxdy dvr \M^ + J 

X {Y+W^{x,Q^) ±Y.xW^{x,Q^) - y^W^{x,Q^)} , 


with being the muon decay constant, and the charged current structure 
functions denoted by >V 2 , TVs and Wl- 
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2.5 Structure functions 


In the naive quark parton model, the deep inelastic scattering cross section 
is calculated as an incoherent superposition of the electron scattering on free 
partons. The electroproduction structure functions are expressed in terms of 
probability distributions of massless quarks and antiquarks in the proton, 

= 2xF^^^{x) =x'^\Qi\^ [qoiix) + mix)] ■ (2.29) 

i 

Here the sum runs over all quark flavors, Qi is the charge of the quark qi in 
units of the elementary charge, and qoi{x) and qoi{x) denote the ‘bare’ quark 
and antiquark distributions. Bjorken scaling is implemented by assuming 
Q^-independence of the bare distributions. As the Callan-Gross relation, 
F 2 = 2xFi, holds in the naive quark parton model, the longitudinal structure 
function vanishes, = 0. 


2.5.1 Structure functions from perturbative QCD 

The calculation of radiative corrections to scattering processes involving 
hadrons in QCD requires the renormalization of the parton densities. This 
results in a logarithmic violation of Bjorken scaling^ and parton densities 
qi{x, Q^) that will weakly depend on the scale Q^. 

In the leading order of QCD, the structure functions Fi^ 2 {x, are given 
in terms of the renormalized parton distributions by 

F^^{x,Q^) = x'^lQil'^ [qi{x,Q^) + qi{x,Q^)] , (2.30) 

i 

Fl^{x,Q‘^) = [qiix.Q"^)+qi{x,Q‘^)\ , (2.31) 

i 


where the qi^x^Q"^) and qi{x,Q^) now denote the Q^-dependent quark and 
antiquark distributions. At this order there is no contribution from the gluon 
distribution g{x,Q‘^), and the longitudinal structure function still vanishes, 
F“(x,Q2)=o. 

Beyond the leading order, these simple relations between the quark and 
gluon distributions and the physical structure functions are modihed. From 


^For an introduction to QCD and hadron structure we refer the reader to the many 
excellent reviews and textbooks, e.g., H and 


46, 4 
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the operator product expansion it follows that they may be written as 

(2.32) 

where the sum runs over all quarks and antiquarks. The coefficient functions 
Cl and Cg, the quark and gluon distributions qi{x,Q‘^), g{x,Q‘^), and the 
QCD coupling as = asiQ"^) depend on the chosen renormalization scheme, 
and they may also become gauge dependent at higher orders. However, it is 
still useful to write structure functions at next-to-leading order as 


F^^^ix, g") = F^^{x, g") + AF{x, gg , 


;iLO 


(2.33) 


with AF{x, Q^) formally being O {as)- 

Common renormalization schemes used in deep inelastic scattering are 
the MS scheme and the so-called DIS scheme. In the DIS scheme the parton 
densities are redefined such that relation (|2.3U|) for F 2 also holds beyond 
the leading order The QCD corrections to F 2 are completely 

absorbed into the redehnition of the parton distributions, so that 


I^DIS 


{x,Q^) = x'^lQil"^ [qi{x,Q^)+ qi{x,Q^)] 


DIS 


AF^^^{x,Q^) ^ 0 . 


(2.34) 


The renormalized quark and gluon densities satisfy the Altarelli-Parisi evo¬ 
lution equations 


dqi{x,Q^) 

dlng2 

dg{x,Q^) 

dlng2 


as{Q^) 

271 

C^siQ^) 

271 


y 


p, 


<?<? 


X 


) (ii{yi Q^) + Pt 

yj 


qG 




P 


GG 


-)^(i/,g^ 

yJ 


with the splitting functions P being given e.g., in K 


X 

yj 


1 y{y.Pt 


(2.35) 


"‘Nowadays it has become customary to denote these integro-difFerential equations as 
Dokshitzer]^ - Gribov-Lipatov|p4 M - Altarelli-Parisi, or DGLAP equations. 
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The NLO correction to the structure function Fi in the DIS scheme reads: 
1 

?DIS/_ /n2^ _ ^ 1,0 |2„, ir^2\ \ f { \ ,n2^ , ,n2^ 


AFr"(a;,Q^) = 


\Qi?as{,Q‘^)^fq ( - ) [qiiy, Q^) + qi{y, Q^)] 

y \i \y. 


2 


gy^]9{y,Q^)\ , (2.36) 


where the sum runs over all active quark flavors. 

At next-to-leading order, the correction terms in (|2.36| ) are gauge inde¬ 
pendent, ultraviolet- and infrared-finite; they read 


(y. Q 4 

asAfq{z) = as [fq,l (z) - fq ,2 (^)] = - 2z , 

c 1 

asAfg{z) = as [fg,i (z) - {z)] = - 4z(l - z) . (2.37) 

The QCD corrections to the structure functions lead to a violation of the 
Callan-Gross relation and to a QCD prediction of the leading twist term of 
the longitudinal structure at next-to-leading order. 


TP r «s(Q^) f dy f ^ ^ 

n(x.Q) = y-y -F,(y,Q 



(2.38) 


X 


ygiy^Q^ 


showing that the ratio R = Fl/{F 2 — F^) ~ O {as) vanishes for —>■ oo 
because of asymptotic freedom. Expression (|2.38|) is also interesting as it 
predicts a strong variation of Fl at low x where the gluon distribution is 
expected to be strongly rising. For a discussion of the gluon distribution and 
further details regarding the DIS scheme we refer the reader to . 

The calculation of the QCD corrections in the MS scheme is extensively 
covered in the literature. The analogous expressions for the electroproduction 
structure functions in terms of renormalized parton densities will not be 
reproduced here, they can be found e.g., in jiSf . 

Since the Altarelli-Parisi equations predict only the evolution of the par- 
ton distributions in the initial conditions have to be fixed otherwise. 
Using certain assumptions about the shape of the parton distribution at 
some low initial scale Qq, which is typically of the order (2 ... 5) CeV^, sev¬ 
eral groups provide parameterizations of the evolved distributions that are 
selected by comparison with experimental data. A large collection of parton 
distributions is available in the program library PDFLIB . 


16 












2.5.2 Structure functions in the region of low 

Exact Bjorken scaling is a built-in property of the naive parton model. Be¬ 
cause of asymptotic freedom it also holds approximately true within QCD in 
the region of large momentum transfers, where it receives perturbatively cal¬ 
culable corrections. However, one should expect that a naive extrapolation 
of perturbative results valid for 3> 1 GeV^ down to the photoproduction 
region ~ 0 will fail sooner or later for several reasons. 

Physically, the spatial resolution power of the virtual photon on the pro¬ 
ton is given by ~ h/ \/^, and therefore a photon of virtuality < {h/rpY ~ 
(0.2 GeV)^, with Vp ~ 1 fm being a typical hadron radius, is not able to 
resolve the individual partons. Therefore, the probabilistic, partonic descrip¬ 
tion of the proton structure must break down at that scale. 

On the technical side, in a systematic analysis of the structure functions 
in the Bjorken limit using the operator product expansion (see e.g., [^), one 
obtains an asymptotic expansion of the structure functions in inverse powers 
of e.g., 

= (2-39) 

The coefficient functions Cnix^Q"^) depend only weakly, i.e., logarithmically 
on Q^. The terms with n = 0 and n > 1 are usually referred to as leading and 
higher twists, respectively. For approximate Bjorken scaling, the contribution 
of the higher twists in ( p.39|) must be negligible, which is the case only for 
of the order of or above a few GeV^. Furthermore, a perturbative evaluation 
of the coefficient functions becomes doubtful at small momentum transfers 
< 1 GeV^ where the running QGD coupling as{Q‘^) approaches unity. 
And hnally, as we shall see below, electromagnetic gauge invariance ac¬ 
tually requires that the structure functions F 2 and Fi vanish in the limit 
> 0. This clearly forbids Bjorken scaling at very low and calls for a 
formal resummation of the series (|2.39|) in that region. 


Constraints for 0 


To investigate the general behavior and to obtain constraints on the structure 
functions for small Q^, let us rewrite the hadronic tensor (|2.19|) in order to 
exhibit potential kinematic singularities as follows ||57||: 


Hi,u{P,q) = dvr 


V -rp Q/lQu 


-QiivFi + — F2 + 

P ■ q 
Pf^qu + Puq^l , 


Fi + ^F2 

qz 


(2.40) 
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As the total cross section for real photons on a proton is hnite, these singu¬ 
larities can only be artifacts of the tensor decomposition of into a sum of 
manifestly gauge-invariant terms. The limit —>■ 0 of this expression clearly 
must exist. This allows to impose the conditions 

F2 = 0{Q'’-) . 

n + = o [Q^) . 

on the structure functions in the limit 0. 

One can go further and relate the structure functions F 2 and to the 
total photoproduction cross section for transverse and longitudinal photons, 

The cross sections on the r.h.s. are functions crr,L = crT,L(Q^, where 
W'^ = -|- 2P ■ q — Q^ is the invariant mass of the hadronic hnal state. The 

Bjorken variable x appearing in the structure functions on the l.h.s. reads: 


(2.41) 

(2.42) 


X = 


1T2 _ + g2 • 


(2.44) 


Invoking conservation of the electromagnetic current, which guarantees the 
decoupling of the longitudinal degrees of freedom of the photon, the cross 
section ctl —0 for g2 q ^^d at fixed P ■ q. We therefore expect 0,1^ 


Fl = 0 (Q*) . 

SO that the ratio R should also vanish in this limit, 

R{x,Q^) = ^ = 0{Q^) . 


(2.45) 


(2.46) 


Obviously it is a crucial requirement for reasonable phenomenological pa- 
rameterizations of the structure functions in the region of low g2 they 
incorporate the above constraints. 


Models for low 

There exist several models for the description of high-energy photon-hadron 
interactions that automatically provide the proper behavior for g2 ^ q \ 
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complete overview is far beyond the scope of this report, so we will only 
name a few popnlar ones, and suggest the interested reader to consult the 
reviews |^, ^ for more detailed descriptions of several phenomenological 
approaches and dynamical models of the F 2 . 

In the Vector Meson Dominance (VMD) model (see e.g.. 


61 and refer¬ 


ences cited therein), the interaction of a photon with a hadron is described 
by assuming that the photon can fluctuate into a vector meson which then 
interacts with the hadron. 

A qualitative description of the low x behavior of cross sections is mo¬ 
tivated by the parameterization of the high energy behavior of total cross 
sections in hadron-hadron interactions in Regge theory |^. In this ap¬ 
proach, the amplitude for hadronic elastic scattering processes at center of 
mass energy s/s = W and momentum transfer t is written as a sum of terms 
of the form T{s,t) ~ with a{t) being a so-called Regge trajec¬ 

tory. Appealing to the optical theorem and the similarity of hadron-hadron 
and photon-hadron total cross sections, the high energy behavior of the total 
photoproduction cross section is written as 


f^tot - 




2\ai-l 


(2.47) 


where the a* = 0 ,( 0 ) and jSi = /9i(0) denote the intercept and the coupling 
of the corresponding Regge trajectory, respectively. Phenomenologically, the 
energy dependence of the total hadronic and photoproduction cross sections 
are fairly well described by taking into account two terms in (|2.47|) : the (soft) 
pomeron with intercept a-p ~ 1.08, which is identihed with an exchange of 
vacuum quantum numbers, and the Reggeon with intercept 0. 5 0 . 

corresponding to an exchange which is odd under charge conjugation. 

Using the relation between the total photoproduction cross section and 
the structure functions (p.43|) , the Regge parameterization (p.47|) then im¬ 
plies the following ansatz for a parameterization of the structure function 
F2(x,Q^), 

= (2.48) 


which is expected to be valid for high energies and for In this case, 

X = Q^/(IU^ — -|- Q^) — Q‘^/W‘^ -C 1. Therefore, the parameterization 

(|2.48|) implies a powerlike small x behavior of ^ 2 ( 0 ;, Q^), 


F2(x,Q'‘) = Y,h{Q'‘) (Q"-) 


Oi-l 


(2.49) 


Furthermore, from (p.48|) and (p.43|) we hnd that the functions Pi{Q‘^) 
to satisfy the conditions /3i{Q'^) = O (Q^) for 0. 


have 
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The HI collaboration has extracted the derivative 


|6l 


A(x,Q2) 


91nF2(x, 

d\n.x 


(2.50) 


and found A(a:, to be independent of a: for x < 0.01 within the experi¬ 
mental accuracy. This implies that the leading x-dependence of F 2 at low x 
is consistent with the power-law behavior F 2 oc x~^ for fixed Q^. 

However, experimental data from HERA on J/'^ photoproduction and the 
charm structure function of the proton show a steeper rise at low x, providing 
evidence for a deviation from the simple Pomeron behavior described above. 
A good description of all available data requires to take into account an 
additional contribution from a “hard” pomeron with intercept 1.4 |66| . 

A power-like behavior of structure functions for x —0 compatible with 
the hard pomeron can be obtained from perturbative QCD. Performing an 
all-order resummation of large a 5 ln(l/x) terms leads to the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) [0] evolution equation, Q whose solution in the lead¬ 
ing logarithmic approximation is interpreted as a pomeron with intercept 


«BFKL — 




TT 


(2.51) 


For this reason the hard pomeron is identified with the perturbative BFKL 
pomeron, while the soft pomeron that dominates in the total cross section is 
associated with non-perturbative effects. 

Recently the color dipole picture has attracted a lot attention. This ap¬ 
proach considers the fluctuation of the virtual photon into a quark-antiquark 
system as an effective color dipole, which then interacts with the proton 
through pomeron exchange. The proper low behavior is implemented 
by suitably parameterized wave functions for transversely and longitudinally 
polarized photons. Some model variants also consider the proton as a quark- 
diquark system, i.e., as a color dipole. For an overview see e.g., ^ and 

references cited therein. 

Parameterizations of the structure function F 2 for phenomenological ap¬ 
plications based on the Regge behavior described above and fit to the HERA 
data have been provided by Abramowicz-Levin-Levy-Maor (ALLM) [[70| , |7T| , 

|, Adel et ah m 


Donnachie and Landshoff 
Desgrolard et ah |]74 . 


Capella et al. 


and 


®See e.g., ref. Q for a recent review on BFKL, or the textbook by Forshaw and Ross 
[13 for a modern introduction to QCD and Regge theory. 
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2.5.3 Structure functions at high 


For neutral current scattering processes at high one has to take into 
account both 7- and Z-exchange and their interference. The generalized 
structure functions Q"^) describing the interaction of leptons with charge 
Qi and lepton beam polarization ^ may be decomposed (using the notation 
of 1^) as: 


xJ=’3{x,Q‘^) 


with Qf, = —1 and 


and 


— Q^) + 2\Qe\{Ve + Xtte) k{Q‘^)Gi^2{x , Q'^) 

+ 4 (Ug + Og + 2XVeCle) (Q^)hfl,2(2^, Q^) , 


-2sign(Qi)||Qg|(ae + Xve)K{Q‘^)xG3{x,Q'^) 


+ [2ng 

tte + A(ng + Og)] k,'^{Q‘^)xH3{x,Q‘^)^ , 

(2.52) 

A 

= ^sign(Qz) , 


Ve 

= 1-4 sin^ 9^ , 


Ojq 

= 1, 

(2.53) 

m = 

G,Ml Q2 

8v/27ra(Q2) Q'^ + ' 

(2.54) 


The quantity sin^ 9^ refers to an effective mixing angle for the renormalized 
fermion-Z boson couplings. 

The leading order expressions of the structure functions read: 


F2{x,Q^) 

= x'^lQql'^ [q{x,Q^)+ q{x,Q^)] , 


G2{x,Q^) 

= x'^\Qg\Vq[q{x,Q^)+ q{x,Q^)] , 


H2{x,Q^) 

Q 

2 2 

(2.55) 

xG3{x,Q‘^) 

= x'^\Qq\aq[q{x,Q^) - q{x,Q^)] , 



9 

xH 3 {x,Q^) = x^^^[q{x,Q^)-q{x,Q^)] , 
q 

where q{x, and q{x, Q"^) denote the quark and antiquark distributions. 
The neutral current couplings of the quarks are: 

Vq = 1 - 4\Qq\sin^ 9^ , 
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(2.56) 


ttg = 1 . 

The charged current generalized structure functions, 

xW 3 (x,Q^) = - sign(Qi) ^ ^ xW^‘{x, , (2.57) 

are obtained from the parton distributions by: 

W^{x,Q‘^) = 2x'^[di{x,Q‘^)+Ui{x,Q^)] , 

i 

W 2 {x,Q^) = 2x^[ui{x,Q‘^) + di{x,Q‘^)] , 

i 

xW^{x,Q‘^) = 2x'^[di{x,Q‘^) - Ui{x,Q‘^)\ , 

i 

xW^{x,Q^) = 2x'^[ui{x,Q^) - di{x,Q^)] , (2.58) 


with Ui and dj denoting the densities of the up-type (u, c, t) and down-type 
(d, s, b) quarks. 

At leading order, the longitudinal structure functions vanish. We will 
not discuss here the QCD corrections to the generalized structure functions 
above but refer the reader to the review 0. 
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Chapter 3 

QED Radiative Corrections 
and Radiative Processes 


The cross sections given in the previous chapter correspond to the Born 
approximation in perturbation theory and have to be improved to include 
higher order contributions, i.e., radiative corrections. At HERA energies, 
the radiative corrections in general need to be treated within the Standard 
Model of the electroweak interaction. 

The electroweak radiative corrections to deep inelastic scattering have 
already been extensively discussed in the literature (for a review see [^, 
and references cited therein). It is well known that the numerically most im¬ 
portant parts of the corrections arise from virtual photon corrections, photon 
emission and light fermion loops. Furthermore, in the kinematic region of 
low to moderate momentum transfers, only the QED part 

of the electroweak corrections is relevant. In this work we will be concerned 
only with QED radiative corrections to scattering processes. 

This chapter starts by briefly reviewing the most important features of the 
O (a) electroweak radiative corrections, but refer the reader who is interested 
in details to the literature. We will then focus on exclusive radiative processes 
at this order where the emitted photon is measured. As we shall argue, these 
radiative processes are of particular interest by themselves. Not only do 
they constitute an important part of the complete radiative corrections, but 
their measurement gives access to physical information that may be difficult 
to obtain otherwise. The treatment of QED corrections to these radiative 
processes, which are a subset of the higher order corrections to deep inelastic 
scattering, will follow in the next chapters. 


771 
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3.1 Classification of radiative corrections 


The basic process in deep inelastic scattering is lepton-quark scattering. It is 
well known that the O (a) radiative corrections to the neutral current cross 
section can be classihed in the following way [Q: 


(1) the leptonic corrections are described by Feynman diagrams containing 
an additional photon attached to the lepton line, i.e., the virtual cor¬ 
rections from the purely photonic correction to the lepton vertex, the 
photon loop contribution to the self-energies of the external leptons, and 
the photon emission from the lepton line, 

(2) the quarkonic corrections are represented by diagrams with an additional 
photon connected to the quark line, analogous to (1), 

(3) the lepton-quark interference part, consisting of box diagrams with at 
least one photon, and the interference part of the bremsstrahlung off 
leptons and off quarks, 

(4) and the purely weak corrections, that are given by all other diagrams that 
do not contain an additional photon. 


Each of the above classes can easily be seen to form a gauge-invariant (w.r.t. 
QED gauge transformations) subset of the full set of corrections.^]’^ The 
infrared-singularities of individual contributions therefore cancel within each 
of the classes (1-3), while each correction in subset (4) is infrared-hnite. 

The contribution of each loop correction diagram i to the cross section 
can be expressed by a correction factor hj, dehned via: 

dVi dVeorn ., X 

The total correction from virtual corrections can in principle be obtained by 
summing over the individual contributions ( |3.1| ). However, there are some 
contributions which are known to be universal and large and which should 
be treated in a special way. 

The vector boson self energies are dominated by fermion loops. Their 
contribution is large because the momentum transfers attained at HERA 

^For non-abelian gauge theories there exists an algorithm for the construction of the 
minimal gauge-invariant subsets of the tree amplitudes at a given order . 

^Class (4) is not irreducible. For example, one may treat the fermion-loop contributions 
to the vector-boson self energies separately. 
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are much larger that the masses of the light fermions, leading to large loga¬ 
rithms ~ In Q'^/mj. The fermion loops of the photon self energy (the vacuum 
polarization) can be taken into account by use of the running QED coupling 
a{Q‘^), thereby summing these leading logarithms to all orders: 


a{Q^) 


a(0) 

i-n](Q2) • 


(3.2) 


Here nJ(Q^) = denotes the fermionic part of the vacuum po¬ 

larization derived from the renormalized photon self-energy S'’', and q;(0) ~ 
1/137.036. Since the hadronic part is notoriously difficult to calculate for 
small Q^, it is preferably evaluated instead via a dispersion relation and pa¬ 
rameterized for practical applications, see e.g., Burkhardt and Pietrzyk |]79| . 


Another set of contributions to the corrections containing large logarithms 
arises from the infrared-divergent virtual corrections and from photon brems- 
strahlnng. This part in general does depend on details of the experimental 
setnp, i.e., whether the radiated photons can be resolved in the detector or 
not. We shall address this issue later. 

According to the above classihcation of the Feynman diagrams, one dis- 
tingnishes the bremsstrahlnng contribntions into the classes (1~3). The ra¬ 
diation from the lepton line, which belongs to class (1), is completely in¬ 
dependent of the modeling of the hadron side.^ The approach of radiative 
corrections that restricts to the leptonic corrections and to the vacnnm po¬ 
larization will be referred to as the model independent approach. The model 
independent framework has hrst been used in deep inelastic scattering in 
]80| . A very detailed acconnt with application to a large set of experimental 
determinations of kinematic variables can be fonnd in |Q. In this chapter 


we shall ontline the essentials. 

In the qnark-parton model, there are also contribntions from radiation off 
the qnark line and from lepton-qnark interference. While we shall comment 
on the former later on, the latter is in general rather small. For details and 
references we refer the reader again to . 


In the case of charged-current processes, there is no gange-invariant sepa¬ 
ration of the Feynman diagrams as in the NC case. A similar decomposition 
can only be performed in a given physical gauge. However, the attribution 
of the leading logarithmic terms of the corrections to the lepton and quark 
lines is still gange invariant [^, Q. Thus, while the leading logarithms of 
the leptonic corrections are model-independent according to the above dehni- 


^It nevertheless depends on the shape of the proton structure functions, i.e., their 
dependence on x and Q^. 
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Figure 3.1: Model-independent contributions to radiative deep inelastic 
electron-proton scattering 


tion, the complete set of radiative corrections to CC scattering is necessarily 
mo del-dep endent. 

The main target of this chapter will be bremsstrahlung in NC processes, 
which we shall treat in the model-independent approach. We will comment 
only briefly on radiation in charged-current processes which has a much 
smaller cross section, and close with a discussion of emission off the hadron. 

In calculations of cross sections below we always take into account the 
large fermionic contributions to the vector boson self energies in the form of 
running couplings. For the neutral current cross section ( p.25|) or ( |2.26| ), this 
implies replacing the prefactor by In the charged current case, 

however, we shall continue to use eq. ( p.28|) , which is accurate enough for our 
purposes. 


3.2 Radiative DIS 

In the following we shall be interested in the model-independent calculation 
of the cross section for the radiative deep inelastic scattering reaction 

e(p) + p(P) ^ e(p') + 7 (fc) + X{Px ), (3.3) 

taking into account only photon exchange^ and emission of the radiated pho¬ 
ton off the electron. This bremsstrahlung contribution is often called the 
Bethe-Heitler (BH) process. The contributing Feynman diagrams are shown 
in hgure |3.1| . The scattering amplitude for this reaction can be factorized 
into the Born amplitude for the process of scattering of a virtual photon on 

"^The cross section for this process taking into account both photon and Z-exchange is 
given in detail in 


26 









the lepton and the matrix element of the electromagnetic cnrrent between 
the proton and the hadronic hnal state, 


_^rad ^ ^ e {Px I Jr(0)| P) , (3.4) 

where q = p — p' — k. The differential cross section for the scattering of 
unpolarized particles can then be expressed as a contraction of the hadron 
tensor with a tensor describing the QED Compton process: 


da 




2y/)\s 2 2 

spins 


\M 


rad 12 


dp' dk 


X 




spins 


dp'dfc . (3.5) 


We shall hrst discuss the expression in square brackets and turn to the treat¬ 
ment of phase space later. 


3.2.1 The QED Compton tensor 


Consider the (sub-)process of Compton scattering of a virtual photon on an 
electron, 

e(pi)+7*(-g) ^ e(p 2 )+7(^) , (3.6) 


and let be the matrix element of this process with the index p describing 
the polarization state of the virtual photon. Adopting the notation and 
normalization of [Q], we define the unpolarized Compton tensor 


Kt^APuP2,k) 


1 

(^ 


^^e7‘-e'7(]Vf«7‘-e'7)* _ 
spins 


(3.7) 


With the help of current and momentum conservation, this tensor is conve¬ 
niently decomposed as follows: 


where 


Q^uPg T ^ ^ PigPjuBij ! 


~ _ (IgQu ~ _ Pi-q • _ 1 o 

Qgu Qgu ^2 ’ q2 ’ * 1 , 2 . 


We introduce the following invariants for this subprocess, 

3^2p2-k, t = -2p2-k, fl=(pi-p2)V 


(3.8) 


(3.9) 


(3.10) 
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so that, using pl= p\ = w?, = 0, we have 

s + i + u = . 


(3.11) 


The expressions for the quantities Bij in the Born approximation read: 


5ii 

Bi2 


= [(s + uY + {i + uY] — 2m^q^ 


1 1 

~ ~ 

5 t 


4g2 _ 
si 

= B 21 = 


Sw? 


B 22 — 


4g2 

si 


Srrr 


(3.12) 


Srrv 

si 


In the high-energy limit, when at least two of the variables s, f, u, q^ are large 
compared to m?, one can drop the terms proportional to m^/si, as their 
contribution is suppressed when integrating over any finite part of phase 
space for the real photon. Therefore in the high-energy limit of (p.l2|) we can 
set 


B, - 

Bi2 = 


J. [(s + uf + {i + uf] - 2m\^ (^ + ^) ’ 


B. 
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(3.13) 


3.2.2 The cross section for single photon emission 

With the expression for the lowest order Compton tensor of the previous 
section we can now write the unpolarized cross section ( p.5|) for the radiative 
process as: 


da 


(47ra)^ 


K^,{p,p',k) H>^'^{P,qh)dp'dk 


<ih=p-p'-k 


(3.14) 


Here we have introduced the hadronic momentum transfer 


Ql = -ql, qh = Px- P = p-p' -k , 

to distinguish it from the leptonic momentum transfer 

Ql = -if . m = p-p' ■ 


(3.15) 


(3.16) 


Note that hadronic and leptonic momentum transfer are related via 


Ql-Q] + 2p-k — 2p'-k. 


( 3 . 17 ) 










In the following we shall need kinematic invariants which are the “hadronic” 
and “leptonic” generalizations of the Bjorken variables, 


Vh 

yi 


P-q, P-{Px-P) Ql Ql 

P-p P-P ’ VhS 2P-qh' 

P-qi P-{p-p') Qf 

P-p P-p ’ yiS 2P -qi ' 


(3.18) 


The contraction of the Compton tensor with the hadron tensor can be de¬ 
composed as: 


= Svr [SiFiixh, Ql) + S 2 F 2 {xh, Ql)] • (3.19) 

Retaining only those terms that survive in the limit S S> but keeping all 
terms of order the coefficient functions Si ^2 can be expressed as follows: 


= 

^2 = 


+ 


where 


~ [(s + u)^ + {i + m)2] + 2771^“^ = 


-Bg , 


S 

XhM^ 

q2 


.B 


[xhS {iu)]^ + (1 - yi)[xhQ - yi)S - (s + m)] 


^1. 


R. 
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(3.20) 


s = 2p' -k , i= -2p ■ fc , u = qf , q^ = ql ■ 


The functions Si and S 2 agree with those given in P9 


Let us assume that the kinematics of the scattering process is determined 
by a measurement of the scattered lepton and the radiated photon. We can 
express the phase space integral for the scattered electron in terms of yi and 
Ql see (1^) : 


dp' = 


(dvr) 


■ ^yi dQ? 


(3.21) 


The differential cross section (p.l4|) then reads: 


327r2n'3 ~ 

da = —^[SiFi + S2F2]dyidQldk 

= [SiFi + ^2i"2] dyi dQl E, dE, dfl., 

'kQho 

= -[>^1^1 + S2F2] dxi dyi E, dE, dfi., . ( 3 . 22 ) 

TT Qh 
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In the last line we have change from the set (?//, Qf) to the set {xi, yi). Note 
that the cross section (|3.22|) may be generalized to also take into account Z- 
boson exchange by replacing [SiFi + 52 -^ 2 ] by with the Fi being 

the generalized structure functions, and S '3 given e. g., in We will return 
to it in section 13.41 


Parameterizations of the photon phase space 


In the radiative cross section (|3.22 ), we have chosen to express the phase 
space of the photon in terms of variables E^,'d,(j) that are measured in the 
HERA frame. Although Lorentz invariance is not obvious, this representation 
is useful for practical applications, as it allows an easy implementation of 
experimental cuts and constraints given by the detector. We will therefore 
only sketch the relation to the manifestly Lorentz covariant approach, which 
is described in great detail in |[5T[].p| 


Let us choose a coordinate system for the HERA frame such that the unit 
three-vectors in the direction of the incoming and outgoing lepton and the 
radiated photon are as follows: 


ei = n, = (0, 0,1) , 

62 = (sin^, 0 , cos 6 *) , 

e = (sin-d cos (/), slu'd sin 0 , cos'd) . (3.23) 

We then have 

c = ei ■ 6*2 = cos 6 , 

Cl = e • Cl = cos'd , 

C2 = e-62 = Goos'd-|-s slu'd cos 0 , (3.24) 

with s = sind = a/I — c^- The integral over the solid angle of the photon 
may be alternatively expressed in terms of ci, C2 as 


dr 2 ..j. = d(cos'd) d 0 = J(ci, C 2 ) dci dc 2 , 


(3.25) 


with the Jacobian 


T(Ci,C2) = 


|ssisin0| (p. - c\ - c\F 2cciC2 ’ a/P ’ 


(3.26) 


®For a classic textbook on particle kinematics and the treatment of multi-particle phase 
space see jS^. 
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Figure 3.2: Schematic view of the kinematic range for the integration over 
the angles Ti and T 2 for hxed r. The allowed range corresponding to "D > 0 
is given by the interior of the ellipse. 


where si = -^/l — cf. The allowed region for (ci,C 2 ) is given by the require¬ 
ment that V > 0. This representation is obviously symmetric w.r.t. incoming 
and scattered lepton, i.e., Ci <->■ C 2 . 

Furthermore, it is convenient to change variables to 


1 - Cl ,2 1 - Cl ,2 

"^li 2 ~ 2 ’ ^ ~ 2 ’ 

(3.27) 

so that 

A 


dfl^ = —= dri dr 2 . 

Vv 

(3.28) 

Expressing T> in terms of ri, 2 , we End: 


V = -ATT 1 T 2 - tI + tI - 2tti - 2tt2 - 2tiT2\ ■ 

(3.29) 

The domain V > 0 corresponds to the interior of the ellipse shown 
which can be parameterized as follows: 

in Eg. 13, 

0 < Ti < 1 , T^iTi) <T2< T^{Ti) , 

(3.30) 
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which is derived from the factorization 

V = (r+(ri) - Ta) (ra - r“(ri)) , (3.31) 

with 

= '^(1 - 'Ti) + ri(l - r) ± 2A/r(l - r)ri(l - n) . (3.32) 

Similar expressions are obtained under the exchange Ti r 2 . 

We assume now the leptonic determination of the kinematic variables, i.e., 
we keep xi,yi fixed. The upper limit on the energy of the radiated photon 
is calculated as follows. The true invariant mass of the hadronic final state 
reads: 

= Ml = {P + p-p' - kf 
= + (1 — xi)yiS — 2k ■ {P + p — p') 

= Wl,-2k-{P + p-p') , (3.33) 

where is the reconstructed (apparent) hadronic mass if only the scattered 
lepton is measured. The true hadronic mass in inelastic scattering must be 
higher than the threshold for pion production, 

Ml>M‘^ = {M + m^f . (3.34) 

This leads to 


2k-{P + p-p') < + {I - xi)yiS] - 

< + 2P ■ (p - p') - gf] - . (3.35) 

Expressing this kinematic limit in the HERA frame we obtain, assuming 
Ep M, Ee^ m: 




^[1 - Y{1 - r)] - - (M^ - M^) 

4 [Ep{l - Ti) + EeTi - YEeT2] 


Y 

T 


— = l-yi + xiyi— , 

Qf _ xiyiEp 

AE^Y (1 - yi)Ee + xiyiEp ' 


(3.36) 


(3.37) 


32 



Using the above parameterization for the photon angles, we can calculate 
the kinematic variables occurring in the Compton tensor (|3.12|) and in the 
coefficient functions 51^2 in the cross section (^.221): 


s = 2p' ■k = 2EeE^{l - /3'c2) ~ AE^E^ (^2 + ^ 2 ) ’ 

i = -2p ■ k = -2EeE^{l - f3eCi) ~ -AEeE^ ^ri + , 

u = -Q2, (3.38) 

Ql = -q'^ = - {s+ i + u) ^ Qf + AE^E^Ti - 4EeE^T2 , 

X, = ^ Ql 

yiS-2P-k yiS-AE^E^{l-Ti)- 


Note that we explicitly kept the leading lepton mass terms in s and t. It 
is necessary to retain them when they occur as denominator factors and we 
want to integrate over the regions ri —>■ 0 or r 2 —>■ 0, where the emitted 
photon is collinear to either the incoming or the outgoing lepton. 

Inspecting ( |3.38| ) and remembering (|3.28|) one easily sees with a little 
algebra that we can express the integration over photon phase space in terms 
of invariants, ^ 

dk = J’' ds d|f| dxh , (3.39) 


with the Jacobian J'' being given in m- With the help of relations like 
Q\ ~ one may switch to other equivalent sets of kinematic invariants 

as integration variables. 

For a fixed energy of the radiated photon, the cross section (|3.22|) ob¬ 
viously exhibits a strong dependence on the integration variables in three 
domains of phase space: (1) ri —>• 0, (2) T 2 —> 0, (3) —>• 0. One can 

decompose the cross section so that each contribution peaks in only one of 
these regions (the peaks were called “s-peak”, “p-peak” and “t-peak” in [^). 
We shall now investigate this peaking behavior in a simplihed fashion. 


Emission collinear to the incoming electron 

In the region ri —*■ 0, where the photon is almost collinear to the incom¬ 
ing electron (“initial state radiation”, ISR), the kinematic variable t of the 
Compton subprocess becomes very small; it approaches |f| ~ O (ttP) rather 
than being of order Qf. As a consequence, the expression for the differential 
cross section simplihes signihcantly. Defining the variable 

Ee — E^ 


(3.40) 
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which describes the relative energy fraction of the incoming electron after 
emission of the photon, we obtain the relations 


Qh ^Ql y Vh yi (1 

This implies 0 < yh < z and {I — z) < yi < 1. 
The coefficient functions 51^2 reduce to: 


xiyiz 


Ql 


yhS yi + z-I 


(3.41) 


Z 


1 + z^ 


1 - ^ (-t) 


2z^ 


(3.42) 


^2 ^ 


5 


h 


XhzS (1 


M^Q 




[zsy 


Q 


1 + z^ 


1 - ^ i-t) 


2z 


m 

~¥ 


Obviously, these functions factorize into the respective Born expressions S'f, 
Q) and a universal function given in parentheses. 


S 2 (see 


For initial state radiation and at high electron energy it is convenient to 
write the photon phase space as 




(3.43) 


The radiative cross section ( |3.22|) schematically simplihes to: 

da ~ da^°™(a:fe, yn/z, Ql) ■ ^ dz d|t| ^ . (3.44) 

When the radiated photon is not observed, we have to integrate over all 
available phase space as discussed above. The integral over t leads to large 
logarithm in the electron mass, of the order InlEl/m^) resp. ln{Ql/m^): 


d^ISR 


/ 1 + z^ In^ - Y 
27r\l — z m2 ^ — z J 




(3.45) 


The soft photon singularity encountered in the integration over z for z —>■ 1 
is canceled by matching infrared singular terms in the virtual corrections. 


Emission collinear to the outgoing electron 


For T 2 —>■ 0, the photon becomes collinear to the scattered electron (“final 
state radiation”, FSR). In this case we have s ~ 0{m?). Defining analo¬ 
gously 




E'e 


(3.46) 
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we obtain 


(3.47) 


Ql = z/Ql , yh = l- Zf{l - yi) , Xh = 


xiyiZf 


1 - Zf{l - yi) 
The coefficient functions again exhibit factorization, 


^2 ^ 


Ql 


^ + z'f I rn? 

Zf \Zf-l s 


S 


ZfQ 


h 


XhS U 


M^Ql 


S2 


Q 


1 + zi I m? 

7 ~ ~Z2 

Zf-ls 


(3.48) 


We write the photon phase for hnal state radiation as 


dk 


(4vr) 


■ dzf ds 


^ ’ 


(3.49) 


with the azimuthal angle 0' measured w.r.t. the outgoing electron. The 
radiative cross section (|3.22|) now reads: 


dcr ~ da^°'''^{xh,yh,Ql) 


a 


1 + ^/ 
Zf-1 



dzf _ d0' 

ds- 

-3 271 


z 


(3.50) 


/ 


The integration over the photon solid angle again leads to a large logarithm in 
the electron mass, In(i7e/n7,)^. The soft photon singularity for Zf ^ 1 cancels 
against the remaining infrared divergent parts of the virtual correction. The 
sum of virtual corrections at relative order O (a) and the contributions from 
initial and hnal state radiation at the same order is infrared hnite, as required 
by the Bloch-Nordsieck theorem. 

The above discussion assumes that the detector is able to distinguish 
the scattered electron from an almost collinear photon. This is the case 
if the detector determines the momentum of the electron directly, e.g., by 
the bending of its trajectory in a magnetic held. We will refer to such a 
measurement as an exclusive one. 

On the other hand, electromagnetic calorimeters essentially provide in¬ 
formation about the deposited electromagnetic energy and the location and 
prohle of the generated showers. When the showers of nearby hits overlap, 
it may no longer be possible to distinguish the particles that initiated these 
showers, and only the sum of their energies can be determined accurately. 
For the present work we will assume for simplicity that a minimum angle 
exists, below which electrons and photons cannot be separated anymore, and 
denote this type of measurement as a calorimetric one. 
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In the calorimetric case, collinear final state radiation therefore does not 
change the measured kinematic variables if the momenta of final electron and 
photon are combined, = p' + k, since 


QLi = -ip- pLiY = Qh ) a;cai = 


Q 


cal 


= Xh . 


2P ■ (p - 

and similarly pcai = Vh- Using p' ~ Pcai/^/) rewrite eq. (p.21|) as 

1 


(3.51) 


dp'^ 4 dp'c»i = 4 12 

Zf Zf (dvr)^ 


dpcal clQ(,a,i , 


(3.52) 


-/ 

and we identify the leptonic kinematic variables with the measured ones. 
After the appropriate variable transformation in the Born cross section, the 
radiative cross section becomes proportional to the non-radiative one: 


da ~ 


271 \Zf — 1 S ^ ) zj 271 


(3.53) 


Performing the integration over the photon solid angle, integrating over Zf 
and combining with the virtual corrections, one can show that - besides 
the cancellation of infrared singularities as in the exclusive case - the large 
logarithm in the electron mass cancels for final state radiation, and only 
a logarithm in the resolution parameter remains, in accordance with the 
Kinoshita-Lee-Nauenberg theorem . 


The QED Compton peak 

A moderately strong peaking behavior of the cross section is also found for 
—> 0. In this case the exchanged photon between the hadron and the 
electron becomes almost real; its momentum is almost collinear to the in¬ 
coming proton. The contributing region in phase space lies in the vicinity of 
the upper boundary in figure T 2 x^iri), corresponding to the outgo¬ 
ing electron and photon having opposite azimuthal angles and compensat¬ 
ing transverse momenta in the HERA frame. The hard scattering process 
with large momentum transfer occurs between this photon and the electron, 
leading to its denomination as the QED Compton process in deep inelastic 
scattering. Approximating the four-momentum of the exchanged photon as 

qh ~ -^P + q±,h , (3.54) 

such that 4 represents the momentum fraction of the proton taken by this 
photon, and q±^h being a small transverse momentum, we find for the invari¬ 
ants of the electron-photon scattering process: 

t ~-4(1 - 2//)S', (3.55) 
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Assuming that we determine the kinematic variables from a measurement of 
the scattered electron, xi ~ The coefficient functions factorize into one 
piece proportional to the differential cross section for Compton scattering, 


^2 ^ 


1 + (1 - yi? 
^ — y I 

1 + (1 - yif 

^-yi 


(3.56) 


Xh - xi XhM'^ 


xr 


Q 


and one can express the differential cross section in the form [^, 

27ra'^ 1 + (1 - yiY 


da ~ 


S xi{l - yi) 


'y{xi,Qf)dxidyi , 


(3.57) 


with 'j{xi,Q‘f) being the inelastic contribution to the photon distribution 
within the proton 





XI (Q^)™ 






(3.58) 


With the integration limits (^Q 2 ^min,max gjyg^ below in (|3.61|) , it is evident 
that the logarithm obtained from the Q^-integration, does not 

depend on a small mass. It therefore is much smaller than the corresponding 
logarithms in the electron mass for ISR and FSR. Nevertheless, the QED 
Compton contribution is signihcant in the region of large yi. 

An accurate calculation of the photon distribution, 'y{x,Q‘^), requires a 
detailed modeling of the contributions to the proton structure functions from 
low mass hadronic hnal states. However, it should be noted that also elastic 
proton scattering contributes to the QED Compton process at HERA, and 
it even dominates the inelastic contribution [Q . A separation of elastic and 


inelastic contributions may be difficult but possible via the residual trans¬ 
verse momentum of the hadronic system. On the other hand, the elastic 
contribution is essentially independent of Q^. Provided sufficient statistics, 
the QED Compton process could also be used as a test of QCD, as the 
evolution of the photon distribution within the proton differs from that of 
the colored partons 10 ■ No separation of elastic and inelastic contributions 
would be required for that purpose. 
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3.2.3 Kinematic effects of photon radiation 


It was already mentioned in the introduction to this chapter that the radiative 
corrections are sensitive to the shape of the structure functions. Following 
Krasny |^| we shall illuminate the connection of this shape dependence to 
the contributions with emission of hard photons. We assume that only the 
scattered electron is measured. 


Kinematic limits 

If we are only interested in the differential cross section in the lepton vari¬ 
ables, d^cr/dx/d?/;, we may integrate over the phase space of the photon. 
When no experimental cuts are to be applied it is convenient to express this 
phase space integral in terms of kinematic invariants, e.g., x^, and s or 
t, see (|3.39|) . Then the integration over s (resp. t) may be performed analyt¬ 
ically. The domain of integration for the other two variables is constrained 
by the inequality (ref. |Q, eq. B.4): 


[ylQl + ylQl - ywh^Qf + Ql)] - Qlf < o . (3.59) 

Taking (x/i, Q\) as remaining integration variables, the boundary of the do¬ 
main described by (|3.59|) reads: 


X. 


= Xi , 


m 


2\max(min) _ 


= Xh 


{xhyiS - Qf) {yiS ± a/\) + 2xhM‘^Qf 


where 


2[xhyiS - + xlM^] 

A, = yfS^ + AM^Q^. 


(3.60) 


In the high-energy limit, S 3> M^, the boundaries for following from 
( |3.6(]| ) are 


(Q 


2 Amin 


XhxfM'^ 


Xh -xiil- xlM‘^/Q^i 


^ XhyiS = 


XhQf 

Xl 


(3.61) 


The upper limit on Xh and a further restriction on is obtained from the 
inelastic threshold: 


Ml = M^ + ^——Ql > . 

Xh 


(3.62) 
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Obviously, = 1 corresponds to elastic scattering of the proton. For > 1, 
the upper limit (|3.61|) then implies 


~ 1 


Xl ( m 2 - m 2 ) 


Qi 


(3.63) 


and the actual lower limit on Q\ reads 

•^2 \min 


(Qi 


~ max 


x,-x,{l-xlMyQ])'l-x„'' > 


(3.64) 


Alternatively, one can switch to {xh, Uh) as integration variables. In that 
case, 


vr ^ 


xhiM^/S 


{xh - xi)yi + xlM'^/S ’ 


Vh = yi- 

Neglecting the proton mass, the above conditions simplify to: 

Xl <Xh<l , 0 <yh<yi . 


(3.65) 


(3.66) 


The radiatively corrected cross section 

As discussed above, the measured differential cross section depends on the 
structure functions in the whole domain Xh > xi and y^ < yi- At leading 
order this dependence can be written in the form of a convolution. 




y.max 


yi 


dxi dyi 


dxh / dyh K{xi,yi;xh,yh) 


da^°’^’^{xh,yh) 

dxh dyh 


(3.67) 




The analytic form of the kernel K{xi,yi] Xh,yh) can be found e.g., in 


We do not need its precise form here, but the relevant approximations have 


been given in section 3.2.2 


The implications of hard photon radiation are most easily understood 
by considering the following simple examples |^. The boundaries of the 
integration region contributing to the inclusive cross section at xi = 10“^, 
= 95 GeV and at a:/ = 9 ■ 10“^, Qf = 95 GeV are shown as dashed lines 
in fig. |T3| a and fig. |^b, respectively. The shaded areas roughly correspond 
to the kinematic domains accessible to fixed-target experiments (upper left 
corner) and the HERA experiments. The HERA band is essentially limited 


39 




















Figure 3.3: The topology of the (x, domains contributing to the radiative 
cross section at HERA, (a) xi = 10“^, = 95 GeV; (b) xi = 9 ■ 10“^, 

= 95 GeV. (Taken from ||]) 
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by the measurement of the scattered electron energy and the geometrical 
acceptance of the detectors |3.|^ The solid curve on the left side of each 


point {xi,yi) corresponds to the integration region contributing in the case 
when the photon is collinear to the incoming electron; it is described by the 
parameterization (see eq. (|3.41|) ): 


Xh = 


xiViz 


yi + z-I ’ 


Qh — zQI 1 


(3.68) 


with Hi = Qf/xiS and z = — E^)/Ef>. Each asterisk on the curve rep¬ 

resents a z{E.y), with E^j = 3j GeV. The asterisk closest to the full circle 
corresponds to j = 1, i.e., E^ = 3 GeV. The upper limit on the photon 
energy in this case follows from the condition (1 — yi)/{l — Xiyi) < z <1. 

Similarly, the dotted curve on the right side of the {xi,yi) point corre¬ 
sponds to the case where the photon is emitted collinear to the scattered 
electron. Here, 


_ xiyiZf 
^ l-Zf{l-yi) 


Ql = ^fQi , 


(3.69) 


with yi as above and Zf = {E'e + E^)/E'e. The kinematically allowed range 
is 1 < z/ < 1/(1 — yiil — xi)). The symbols on the curve correspond to the 
same choice of photon energies as above, E^j = 3j GeV. 

The dominant contributions to the measured cross section ( p.67| ) orig¬ 
inate from small bands centered around the curves discussed above; they 
correspond to the domain where the kernel K is large. The scale of width of 
these bands can be estimated using the approximations (|3.44|) and ( |3.53|) ; for 
the photon solid angle they correspond to cones of half opening angle ~ m/Eg 
and ~ ml Eg around initial and hnal electron, respectively. Another signih- 
cant contribution to the cross section is the region where the Born differential 
cross section in the integrand is large, i.e., where Q\ is close to its lower limit 
(|3.64|) . In this region the photon exchanged between the electron and the 
hadronic system is almost real, and the hard scattering process is interpreted 
as QED Gompton scattering, ey —>■ ey, between the electron and the photon 
emitted from a quark or proton. 

As can be seen from fig. |3.3| , moving the point {xi,yi) around within the 
HERA kinematic domain, one can easily find regions where the dominant 
contributions to the convolution integral (|3.67|) involve mostly “uncharted 
territory”. This is especially the case for small x < 10“^ and < 10 GeV. 


®The kinematic domain covered by the HERA experiments has been extended over 
time, mostly towards lower values of x and Q^, by successive upgrades of the detectors 
(e.g., small angle taggers) or using a longitudinally shifted position of the interaction 
vertex. However, an ‘unexplored gap’ between the shaded domains still remains. 
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At this point uncertainties enter the radiatively corrected cross section. 
While one might expect reliable predictions for the structure functions from 
perturbative QCD already for moderate momentum transfers, say, above 
4 GeV^, this appears questionable when the non-perturbative regime ;< 
1 GeV^ is approached. Therefore, extrapolations of the structure functions 
such as those mentioned in section |2.5.2| will be needed, leading to a model- 
dependence and uncertainty in the predictions of the radiative corrections. 

We would like to point out that in the case of HERA this strong model- 
dependence concerns mainly initial state radiation and the Gompton process. 
The Gompton contribution to the structure function measurement can how¬ 
ever be suppressed by imposing some minimum experimental cut on the 
transverse momentum or invariant mass of the outgoing hadronic system. 
For hnal state radiation, the calorimetric measurement of electrons and pho¬ 
tons leads to a cancellation of the large logarithmic contribution as discussed 
above. It also does not shift the kinematic variables. 

The uncertainties related to the model-dependence can be brought under 
control if one can experimentally detect events with emitted hard photons 
and either rejects them (see e.g., the discussion in for using tagged ISR 
at HERA), or if one turns necessity into a virtue and actually accepts these 
events and recalculates the kinematic variables of the hard subprocess taking 
into account the four-momentum of the emitted photon. This way one can 
actually extend the kinematic region in which the structure functions of the 
proton can be measured. Furthermore, we shall see another virtue of the 
measurement of radiative events in the next section in that it helps in the 
determination of the longitudinal structure function 

The above reasoning in principle also applies to other methods of de¬ 
termination of the kinematic variables, although the relations between the 
measured kinematic variables and those of the hard scattering process differ. 

Finally we would like to mention that the neglect of photon emission off 
the hadrons in the above discussion is not a problem. We shall argue in 
section that initial state radiation can be absorbed into the evolution of 
the parton distributions, while hnal state radiation off hadrons does not shift 
the kinematic variables at all. 


3.3 Collinear radiation and tagged photon 
processes 

The measurement of deep inelastic scattering with an exclusive photon can 
be expected to be a challenging task, as the radiative cross section (|3.22|) is 
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Figure 3.4: Schematic view of the HI luminosity system. (Taken from PB[]) 


formally suppressed by a factor o/tt relative to the non-radiative cross sec¬ 
tion. We mentioned three regions of photon phase space where the radiative 
cross section is enhanced and can partially compensate the large suppres¬ 
sion. However, only the emission collinear to the incoming electron appears 
to have been exploited yet. 

The main reaction used for the determination of the luminosity at the 
HI and ZEUS experiments at the HERA collider is the radiative elastic 
scattering process ep —*• ep-l- 7 . For the identihcation of the emitted photon, 
both experiments possess a small photon detector (PD) in the very forward 
direction, see £g. for HI and see fig. for ZEUS. The solid angle covered 
by these photon detectors is roughly described by 0 < 'd-,, < "do about the 
direction of the initial electron beam, with the maximum angle 'do being of 
the order of (0.45... 0.5) mrad for both experiments. 

This forward photon detector is an ideal device for tagging radiative deep 
inelastic scattering events. Since the maximum angle "do is small compared to 
the scattering angle of the electron 6 that is measured in the main detector, 
the differential cross section for these radiative events factorizes exactly. As¬ 
suming azimuthal symmetry, integrating the Compton tensor over the solid 
angle covered by the photon detector yields: 


71 


I dD., = 


PD 


1 + z'^ 
1-z 


In 1 -h 




)- 


2z 


m? J 1 — z 


1 + 


2 \ 
m \ 

'WoJ 


X - + ^zp^pj) + O (do) 


(3.70) 


The tensor structure in the last line is proportional to the non-radiative 
lepton tensor ( |2.17| ): 

—^—-L^,^{zp,p'). (3.71) 

Zj 
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Figure 3.5: General layout of the ZEUS luminosity monitor. (From p0[|) 


We therefore hud: 
d^cr 


dxi dyi dz 


X 



2 ; 


a 

1 + 

271 

l-z 


(3.72) 


J 1 — z 


1 + 


m 


2 \ -1 




Inserting the conditions at HERA, we hnd for Eg = 27.5 GeV and i9o = 
(0.45 ... 0.5) mrad: 


f;2 q2 

Co := ^ (6 - 7) ■ 10^ » 1 , 




(3.73) 


thus Co S> 1 even if i9o -C 1. Neglecting terms of order O (Co and O ('do), 
the cross section (|3.72|) simplihes to: 


d^CT 


dxi dpi dz 271 
where we have introduced 


— G, Lq) ■ - 


dxh d{yh/z) 


1 “|~ ‘^z 

Lq) = - —Lq — - -^ , Lq = In Co 


l-z 


1-z 


(3.74) 


(3.75) 


Evaluating the “large logarithm” Lq for the HERA forward photon detectors 
we obtain: 


Lo ~ 6.5 . 


(3.76) 
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The cross section (|3.74|) becomes even simpler when we interpret the emission 
of the photon as a reduction of the electron energy, Ef^ = zE^.. The energy 
fraction can be invariantly written as: 


z = 


2P ■ {p — k) 


(3.77) 


Furthermore we choose the following set of invariant kinematic variables that 
use the scattered electron and take into account the energy loss due to photon 
emission 1113: 


= -{p - p' - kf 


Obviously, 


X = 




2P ■ ij) — p' — k) ’ 


= xyzS 


y = 


P ■ {p — p' — k) 
P ■ {p — k) 


The physical range of the shifted Bjorken variables reads: 

0 < x,y < 1 . 


(3.78) 

(3.79) 

(3.80) 


The shifted variables (|3.78|) are related to the standard Bjorken variables via: 


= zQ 


i 1 


X = 


xiyiz 


y = 


yi + z-1 


yi + z-l ' " 2 ; 

and with the Jacobian for the variable transformation: 


d{x,y) 


d{xi,yi) 


yi 


yi + z-I 


zy 


(3.81) 


(3.82) 


In terms of shifted variables, the radiative cross section (p.74|) assumes a very 
simple, suggestive form. 


d^cr _ ^ ^ ^ dVs orn 

dxdydz 2 ti ° dxdy 


(3.83) 


The differential Born cross section on the r.h.s. is a function of the shifted 
variables (|3.81| ) only. For S', x and fixed, the variable y = Q‘^/{xzS) can 
still be varied through its dependence (|3.77|) on the tagged photon energy. 
Thus a measurement of cross section (|3.83| ) for different z can be used to 
separate P 2 and El at fixed center of mass energy |^, see also (p.25|) . Fur¬ 
thermore, from ( 3.81 ) it is obvious that smaller values of Q'^ can be attained 
for fixed energy and angle of the scattered electron than for the non-radiative 
measurement. 


45 



































Figure 3.6: The lowest order radiative cross section (|3.83| ) for = 10 GeV. 


ure 


To give an impression of the radiative cross section, we show in £g- 
3^ the differential cross section (|3.83|) for a tagged photon energy of 


= 10 GeV. As parameters we used 


E, = 27.5 GeV 


Ep = 820 GeV 


"do = 0.5 mrad , 


(3.84) 


the ALLM97 parameterization of the proton structure function E 2 , and 
for simplicity a hxed value R = 0.3. Similar results are obtained using 
modern structure functions from PDFLIB 
range allowed by experimental data (see e.g 

(see e.g., |^) can lead to a sizable variation of the cross section for large y. 

For later convenience we shall also introduce here the quantity S which 
is related to the lowest order cross section as follows: 


However, varying R in the 
6 |) or suggested by models 


S(£, y,Q^) = 


1 dVBorn 
y dx dy 


(3.85) 
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so that, together with 


“ h-4V = ^ = T ‘2") • 

Hi dxi ay I az y da; ay az ztt 

3.4 The general radiative DIS process 


(3.86) 


In the kinematic region of high momentnm transfers > M|, which are at¬ 
tainable at HERA, not only photon but also Z exchange have to be taken into 
account. The calculations of section |3.2| can be generalized by reconsidering 
the Compton tensor (|3.7|) and replacing the virtual photon by a left-handed 
vector current and contracting it with the generalized hadron tensor. The 


result is simply a replacement of (|3.19|) : 

^ Stt Ql) + S2J^2{xh, Ql) + S^J^^{xh, QD] . (3.87) 


The jFj are the generalized structure functions (|2.52|) , 51,2 are identical with 
(CT) , and H: 


5.1 = 


XhS ( 2 (2 “ A 1 


m 


TTl 

+ 7 ^ - 2(1 - yi)xhS) - — {q^ + 2xhS) . (3.88) 

L S 


Taking over the derivation of the differential cross section (^.14| ), we obtain: 

2 a^ yi 


dV = 


TT Ql 


SiJ^i dxi dyi dE^ dH.^ 


(3.89) 


This cross section exhibits the same peaking behavior as discussed above. 
For example, in the region of collinear initial state radiation we use the 


approximations ( 3.42 ) for 5i,2, along with 


5.1 ~ 


—I 


1 + z^ 


-t) 


9 

m 

2z^ 

E 


Consequently, the tagged photon cross section at high momentum transfers 
factorizes: 




d 2 ^Born 


NC 


dx d^ dz 271 dx d^ 

with the neutral current Born cross section 


x,^,Q 2 j , (3.90) 

appearing on the r.h.s. 
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For Q"^ <C M|, eq. ( |3.9CI|) agrees well with ( |3.83| ). At large momentum 
transfers > M|, where also Z exchange plays a signihcant role, the cross 
section is rather small, as it is roughly suppressed by a factor a/n relative to 
the non-radiative cross section, leading to a rather small number of events. 
Besides, in this region the QCD prediction of should be reliable, and it 
is more important to measure JF 3 . This, however, is better achieved by using 
different lepton charges rather than photon tagging. 

The radiative neutral current process can in principle also be used to test 
the validity of the Standard Model. In [^] the possibility of a measurement 
of possible anomalous triple gauge boson couplings (Z 77 , ZZ 7 ) was studied. 
However it was found that the sensitivity achieved at HERA is too small. 


3.5 Tagged photons in charged current reac¬ 
tions 


Radiative processes with hard photon emission do also occur in charged 
current reactions. However, the cross section for the non-radiative process 
ep —>■ uX is already quite small, and “paying” another factor a/n further 
reduces it, making it difficult if not impossible to gather sufficient statis¬ 
tics even after the HERA luminosity upgrade. Nevertheless, we shall briefly 
discuss it here for the sake of completeness. It is worth to mention that 
in contrast to neutral current processes the doubly differential cross section 
d^ucc/dxdj/ is still a rather flat distribution in y at HERA energies, as the 
denominator of the W propagator is: 

= xyS + ^ 0{S) ^ O {M^) . 


In the charged current case there are no leptonic variables at our disposal. 
The only option to determine the kinematics is thus the measurement of 
the hadron variables and PT,h-i and the energy of the tagged photon. In 
analogy to the non-radiative Jacquet-Blondel method ( |2.13D we dehne shifted 
hadron variables: 


2/jb — 


OjB — 


Pph 


2zEe ' 1 — ^jB 

Their relation to the unshifted variables is: 


a^jB = 


2/JB A2 ^2 1 ~ yjB 

9JB- — , - OjB • 


a^jB = a:jB 


ymzS 


1 - yjB 
1 - ym/z 


(3.91) 


(3.92) 


This apparently implies Qjb > Qjb^ contrast to = zQf < Qj for the 
lepton method. However, it should be mentioned that this does not mean 
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that one can really access events at higher than in the non-radiative case. 
It only tells us that the kinematic effect of collinear radiation makes radiative 
events (with untagged photons) appear at even lower measured values of . 

As we have argued at the beginning of the chapter, there is no gauge 
invariant subset of the Feynman diagrams for single photon emission in the 
case of charged current processes, and one cannot define leptonic or hadronic 
corrections. Nevertheless, the leading contribution to the radiative cross 
section for emission in the forward direction has the same universal structure 
as in the neutral current case. This may easily be verified by performing an 
explicit calculation, using a physical gauge for the emitted photon. As we 
intend to drop all terms at O ('dg) from the integrated cross section, the only 
relevant Feynman diagram that needs to be taken into account is the one 
with emission from the incoming lepton. 

Let us choose an axial gauge, defined with the help of an arbitrary light¬ 
like four vector c (c 7 ^ fc, = 0). The sum over photon polarizations then 
reads 

5 ^ tyA., AMfc, A) =+ LAthfii . (3.93) 

A=± ^ 

The choice of the vector c is not important as long as we make sure that its 
direction is far outside the cone around the incoming electron. In the present 
case it is convenient to take it proportional to the momentum of the outgoing 
(massless) neutrino, c = p^. 

An explicit calculation confirms the tagged photon cross section, inte¬ 
grated over the photon polar angle in the range 0 < d < do and expressed in 
terms of the shifted variables (|3.91|): 


decree 

dxjB d^jB dz 


a 


P{zM) 


dfjB d^jB 


djB, Q 



(3.94) 


with Qjb = ^JB^JB^'S'. At HERA energies, we may assume that the longi¬ 
tudinal structure function Wl is small and calculable, being proportional to 
Q!5 (Qjb)- a possible beneht of using radiative charged current events, more 
important than extending the accessible kinematic region, may be an im¬ 
provement in the separation of the structure functions W 2 and xWs, which 
also depend on the lepton charge, see (|2.58|) . However, considering the small 
variation of the relative size of the respective coefficients and F_ and the 
projected integrated luminosity of HERA this appears very difficult. 

Non-collinear emission of photons in charged current reactions can also 
be used to check the WW 7 vertex of the Standard Model. Helbig and Spies- 
berger have compared the sensitivity of the reaction ep —a u'-fX to the 


W production process at HERA and concluded that it is too small to be 
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competitive to other determinations, 
see e.g., but also the comment in 


(For constraints from e^e colliders 


3.6 Radiation from the hadron 

In the above discussions we have put aside the QED corrections from the 
hadron side. These correction are in practice quite small, despite the small 
masses of the light quarks, and may be treated perturbatively. We will follow 
here the reasoning of Kripfganz and Perlt 


100 


In calculating the radiative corrections to high energy processes we gen¬ 
erally encounter mass singular terms of the type 




21^ 


2n 




where is a typical large invariant and m a mass of any of the external 
particles. In the case of QCD it has been shown that the mass singularities 
of this type can be factorized and absorbed into the hadronic structure or 
fragmentation functions ||101|| . Performing the factorization at some reference 
scale Qo then shifts all logarithms in the small masses into the distributions, 
and one may expect that the resulting net correction will be typically of 
order a!/(27r) ■ ln(Q^/Qg). Therefore, the hnal result of a perturbative QCD 
calculation does not involve any large logarithm of a small quark mass. 

This procedure can be extended to take into account the inclusive QED 
corrections to structure functions. Starting from the evolution equations for 
the parton distributions (|2.35|), we replace the quark splitting function 


JD ( \ 






aag 


m 




2 <?<? 


27r 

+ 0{a‘^) . (3.95) 

The hrst line in (|3.95 ) represents the usual expansion of the splitting func¬ 
tion in QCD. The results for the can be obtained from the Pqq’^^ by 

replacing the gluons by photons, i.e., substitution of the group theoretical 
factors of SU{3) by the corresponding ones of the U{1), 


Cf 


1, Ca ^ 0, 


and taking into account the quark electric charges. Note that the leading 
logarithmic contributions Pql’^ and Pqq’^'^ are proportional. 


p(o,i) _ 

<?<? 


'1 p(i>o) 
gg 
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Figure 3.7: dependence of the QED corrections to the structure function 

F 2 for deep inelastic lepton-proton scattering at a; = 0.001, x = 0.1 and 
X = 0.505 in percent. (Taken from [102 ). 


As long as we are only interested in the leading order QED correction to the 
quark distributions, we expand them in powers of of a. 


q{x, Q^) = g® (x, Q^) + QQ + . .. , (3.96) 

g{x, Q^) = g^^\x, Q^) + -^g^^\x, Q^) + ... , (3.97) 


where the upper index (0) denotes the pure QCD evolved parton distribu¬ 
tions, and (1) the O {a) correction. To this order we may take a hxed; any 
Q^-dependence is then shifted to order Kripfganz and Perlt explicitly 
demonstrate in ng that the O (a) correction g^^^ obtained from the 
solution to the generalized DGLAP equations with (|3.95| ) do not contain 
logarithms in the fermion mass. 

The QED corrections to parton distributions have also been studied in 
|102|| . Although depending on the input parton distributions, the corrections 
to structure functions in general stay at the per mille level, except for x —> 1 
where they may reach of the order of 1 percent. A typical result for the 


structure function F 2 in shown in fig. 

Photon radiation from hadrons at large angles w.r.t. the hadron is small 
and calculable within perturbative QCD. For the tagged photon process dis¬ 
cussed above, its contribution is O ('dp) and thus negligible. 
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Chapter 4 

QED Corrections to Radiative 
Scattering 


From the viewpoint of perturbation theory, the radiative cross sections dis¬ 
cussed in the previous chapter are interesting hybrids. On one hand, these 
radiative processes are already contained in the radiative corrections to the 
deep inelastic scattering process as part of the hard photon corrections. On 
the other hand, they correspond to the Born approximation for the scattering 
process with the given number of detected particles which is “DIS plus one 
exclusive photon”. For a precise description of the cross sections, radiative 
corrections to these processes need to be considered. Due to the abovemen- 
tioned relation the corrections to radiative processes are necessarily part of 
the higher order corrections to non-radiative processes. 

This chapter is devoted to a calculation of the QED corrections to the 
tagged photon process at the next order in perturbation theory. Note that 
they are of order O {a) relative to the radiative process, but of order O (a^) 
relative to the non-radiative DIS process. Since the QED corrections are 
dominated by large logarithms in the electron mass, we start with a discus¬ 
sion of these corrections in a simplihed framework by focusing just on the 
leading logarithmic terms |^. We then illustrate the influence of the choice 
of reconstruction method of the kinematic variables on the size of the cor¬ 
rections. The qualitative features of the corrections turn out to be similar to 
those in the case of non-radiative DIS. 

We will not repeat here the calculation of the subset of next-to-leading 
logarithms which were presented in [^, for leptonic variables and in 
for the S variables, but explain in some detail the full calculation of the 


model-independent O (a) corrections by the present author |^. A system¬ 


atic approach to QED corrections taking into account the leading logarithms 
at all orders follows in the next chapter. Our restriction to the corrections 
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to the tagged photon process excludes the consideration of e’''e pair pro¬ 


duction corrections in DIS ( |[103|| , but see also ||104|| and references) which at 


this order do not contribute. Furthermore, we will not discuss here the QED 
corrections to the QED Compton process but refer the reader to 


4.1 Radiative corrections in the leading log¬ 
arithmic approximation 

As discussed in the previous chapter, there are two potential sources of large 
corrections to the cross section that contribute terms of the order of a/vr ■ 
ln(Q^/mj), with rrif being the mass of a light fermion. The hrst one are the 
vacuum polarization corrections to the propagator of the exchanged boson 
between the electron and the hadronic system, which we assume to be taken 
into account and resummed using renormalized (running) couplings. The 
other source of large corrections is the radiation of photons off light fermions. 

Following [^], we shall suppose a calorimetric experimental setup, where 
a hard photon radiated collinearly to the outgoing electron line cannot be 
distinguished from a bare outgoing electron, so that final state radiation 
can be neglected to the desired accuracy. We also assume some minimum 
lower cut on the transverse momentum of the outgoing hadronic system, in 
order to suppress the contribution from QED Compton events for a leptonic 
measurement of the kinematic variables. Therefore, the only large corrections 
we explicitly need to consider below originate from the radiation of photons 
off the incoming electron. We will treat these in the so-called structure 
function formalism. 


4.1.1 Structure function formalism 


In the structure function formalism for QED ||105| , p.06| , |107|| , a radiatively 
corrected inclusive cross section is written as a convolution of the electron 
non-singlet structure function with the Born cross section (Tq for 

the hard scattering process taken at reduced center of mass energy. Schemat¬ 
ically, 

o'Rc(s) = j QQ ao{zs; QQ , (4.1) 

with being some typical large scale of the hard process. The electron 
non-singlet structure function. 


= + + pl=>W + 0((aLf) , (4.2) 
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which is a QED analog of the parton distributions in QCD, depends on the 
large scale only via the large logarithm L = ln{Q‘^/ml). It is known to 
properly sum the leading logarithmic contributions (aL)” to all orders in 
perturbation theory ||107| . 


For our purposes we need here the hrst two coefficients of the power series 
expansion ([4.21) of Introducing a small auxiliary parameter e that serves 

as an infrared (IR) regulator to separate virtual and soft from hard photon 
contributions, these two coefficients read: 


- z) + ■ e{l - e - z) 


with 


Pr. 


(1), 


e 


z = 


1 + z^ 


pP = 21ne + 


(4.3) 


and similarly (see e.g. |p.07|| ), 


P. 


( 2 ), 


»l-e 


Z = 


= 2 


1 + Z^ 
1-Z 




( 1 ) 


0 




(4.4) 


21n(l 


, 3 

ln. + - 


I 1 + ^1 

H- Inz 


1 + ^ 


Inspecting the lowest order tagged photon cross section (|3.83|) , 

d^aio _ p, ^ . dVs orn 

dxdydz 2n ° dxdy ' 

and comparing with the O (a) contribution in the integrand to the inclusive 


cross section (hffi) , we conclude that the logarithmic piece of P{z,Lo) is 
indeed contained in the first order correction; the difference in the logarithms 
{L — Lq) is accounted for by the cross section for photons emitted at angles 
larger than 'do- 


4.1.2 Leading logarithms at the next order 

Let us now turn to the contributions to the cross section at order O (a^) |l30| . 
The typical maximum emission angle 'do is about (0.45... 0.5) mrad for the 
HERA photon detectors. The logarithm Lq appearing in (|3.83|) is moderately 
large, Lq ~ 6.5, see (|3.76|). The complement 


Ti = L — Lr 


6 ... 16 S> 1 (for e.g., = 0.1 ... 1000 GeV^) 


is of similar magnitude or even larger than Lq. Starting at this order we thus 
have two competing large logarithms Lq,Li. Asa consequence, in the leading 
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logarithmic approach at O (a^) we need to consider all donble logarithms of 
the type o^LqLi. There cannot be terms a^Ll since we reqnire at least 

one tagged photon. 

The corrections receive contribntions from virtual corrections to collinear 
emission, one collinear plus one soft photon, two collinear photons, and one 
collinear plus one non-collinear photon. Since the logarithmic contribution 
from the non-collinear photon arises mainly from emission at small angles 
outside the photon detector, we shall also denote it as a semi-collinear one. 

The sum of the contributions of virtual and soft corrections to collinear 
photon emission, with the emission angle of the soft photon being integrated 
over the full solid angle, but the hard photon only over the solid angle of the 
PD, can be obtained from the expression for the one-loop Compton tensor. 


The logarithmic part reads [84 


2ti) 


(Lo + f^i)Pf-(f^o + 2Li)ln^ P^^\z)-ao{x,y,z) . (4.5) 




For the double collinear contribution we assume that only the sum of the 
energies of photons in the PD is measured. In this case z = 1 — 

The contribution from two hard photons in the PD, with the energy fraction 
of each being larger than e, then yields ||108||: 


a \ 2 Lq 


2 r 


2tt) 


2 L 


- 2P'‘’(z) 


m, 


-In^ 




(4.6) 


see also ([4.4|) . 

The hnal contribution is associated with one collinear photon hitting the 
PD and the other one being emitted at an angle larger than "do. Let us try 
to understand the origin of the large logarithms in this case. For the hard 
photon that hits the PD, the major contribution comes from the region of 
polar angles~ rrie/Pe ^ 'do- Similarly, for the other photon the biggest 
contribution comes from polar angles close to the lower limit, > -yQ. 
Thus, the leading contributions come essentially from the region where we 
have strong ordering in angles, ~ Co ^ 1- One can easily see by 

an explicit calculation that one obtains the double logarithms entirely from 
the contribution where the photon hitting the PD is emitted hrst and the 
“lost” photon is emitted second, while the reversed case is suppressed as long 
as the following condition is satisfied: 


2^2 ■ Co > 1 , where X 2 = 


E. 


( 2 ) 




(4.7) 


The contributions from regions where condition ( |4.7|) is not met are sublead¬ 
ing. Thus we shall assume in the following e -C 1, but e ■ Co S> 1. 
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Taking the ordering of photon emissions into account, we calculate the 
contribution from one tagged plus one undetected photon using the quasireal 
electron method ||1U9||, 






dxs (1) 
~ -^0 


Mx,y-,z - X 2 ), (4.8) 


where cxo is understood to be expressed in terms of the “true” kinematic 
variables Xt, yt of the hard subprocess, 


a(i{x,y\ z -X2) = aoixtyPfyZ - X2) - J(x,y;l- xzfz) 


(4.9) 


Note that this contribution explicitly depends on the experimental determi¬ 
nation of the kinematic variables x and since the almost collinear emission 
of the second photon shifts the “true” kinematic variables {xt, yt) with respect 
to the measured ones {x,y). The Jacobian 


JT” = det 


/ d{xt,yt) 

V d{x,y) 


(4.10) 


accounts for scaling properties of the chosen kinematic variables under radi¬ 
ation of the second photon. The upper limit of the a: 2 -integration in (|4.8|) is 
given by either some experimental cut on the maximum energy of the second 
photon, or by the kinematic limit, which also depends on the choice of the 
reconstruction method. 


After a change of variables, X 2 = zu, Uq = the integral in (^^) 

may be conveniently decomposed into IR divergent (e dependent) and IR 
convergent contributions as (suppressing the arguments f, y)\ 


dX2 (1) 

pUQ 

Je/z 

Mz) 


(^1 - 0-0(2: - X2) = 


du Pe\l - u) [(o-o(2:(l - m)) - ao{z)) + cro(2;)] 




duP^\l-u) 


uo 


+2\nz-P, 


( 1 ) 


Cro(2:(l - u)) _ 
doiz) 

duP^\l-u) 


' UQ 


(4.11) 


where in the last step we have extended the w-integration of the IR convergent 
piece to 0 due to the smallness of e, and we have used the property 


duP^^\u) =0 
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Method 

Ql 

Xt 

yt 

■ 

^min 

J{x,y;z') 


z'CP 

xyz' 

y+z'-l 

i-y 

y 



y+z'-l 

z' 

l-xy 

y+z'-l 

JB 

i-fijB /n2 

1-§JB y 

ym 

ym 

1-i/jB 

i-yje/z' 


z' 

1-^jb{1-i/jb) 

z'-ym 

S 

Ql 

xs 

z' 


XT, 

1 

2' 


Table 4.1: Scaling properties of the kinematic variables nnder initial state 
radiation for different experimental methods of their determination. 


in the simplification of the IR divergent piece. 

Adding np the contribntions ([4.5|), ([4.6|), and ( |4.8| ), we see that the de¬ 
pendence on the anxiliary parameter e cancels, as it shonld. 


4.1.3 The radiatively corrected cross section 

We are now able to write down the resnlt for the leading logarithmic radiative 
corrections. The correction is dehned via: 


dx d^ dz 


d^CTl 


lo 


dx d^ dz 


(1 +'5rc( 4, a, z)) . 


(4.12) 


Since we restricted onrselves to the leading logarithms, we mnst retain in 
the correction factor (1 -|- 5 rc) only the logarithmic terms Lq, Li from ([4.5| ), 
(|T^ ), and (|4^ ). We thus obtain:^] 


aLo I-z 
Orc — ——^ 


( 2 ), 


ItT 1 -1- x2 ® 




+ 


aLi 


2ti 

f 

J UQ 




du P^q\1-u) 


IJQ 

du Pd)(i_„) 


cro{z{l - u)) 
do^z) 


(4.13) 


Remember that the contribution from the undetected hard photon depends 
on the choice of kinematic variables and on the upper limit uq. 

^In contrast to we do not count here the contribution from vacuum polarization 
as part of the radiative corrections, as it is contained in the leading order tagged photon 
cross section dd a\o /dx dy dz. 
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Figure 4.1: The leading logarithmic radiative corrections to the tagged pho¬ 
ton cross section for = 10 GeV. Lepton method. 


4.1.4 Numerical results 


To give an impression of the signihcance of the radiative corrections at leading 
logarithmic accuracy as given by eq. (|4.10| ), we shall present some typical nu¬ 
merical results. The relations between the measured and the true kinematic 
variables in ( |4.9|) for some reconstruction methods are given in table 


they agree with those for the calculation of the leading logarithms from ini¬ 
tial state radiation to non-radiative DIS (see e.g., |^)- The variable z' 


corresponds to 1 — u in ([4.13|) , so that uq = 1 — z\ 


We use the same set of HERA parameters as in section]^, eq. (|3.84|) , the 
ALLM97 parameterization and a hxed ratio R = F^/Ft = 0.3, unless stated 
otherwise. As a representative intermediate value for the tagged energy we 
take E^ = 10 GeV. No cut will be applied to the energy of the second 
(“lost”) photon. 

Figure ^ shows the radiative correction 5 r,c (EiH) reconstruc¬ 

tion of the kinematic variables using the electron method. In the region of 
very small yi, the radiative corrections may become strongly negative. The 
contributions from virtual and soft photon corrections dominate, because the 
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Figure 4.2: i?-dependence of the leading logarithmic radiative corrections to 
the tagged photon cross section for = 10 GeV. Lepton method. 
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kinematic limit on the energy of the undetected photon tends to zero as 


Uq — 1 — 


yi{^ - xi) 


xwi 


(4.14) 


For yi ^ 1, the phase space for photon emission increases, and we hnd the 
typical strong rise of the corrections, which is due to the large shift between 
the “true” and the measured the kinematic variables under hard photon 
radiation as discussed in section |3.2.3| , and due to the steep dependence of 
the lowest order neutral current cross section for small y, resp. Q^. 

The radiative corrections for the electron method exhibit a very strong 
dependence on the ratio R = Fl/Ft^ see £g. Again, this effect is largest 
for large y. Increasing R increases the correction, mainly because it reduces 
the hard cross section cxo for y —^ 1, while hard collinear radiation probes 
Xt > X, yt < y. This strong dependence of the radiative corrections on a 
quantity (here: R) poorly known in other regions of phase space is a nice 
example of the kinematic effects we discussed in section |3.2.3| . 

Figure |4.3| a shows the leading logarithmic radiative corrections for the 
Jacquet-Blondel method, again for a tagged photon energy of 10 GeV. In 
this case the corrections are negative for ^ 1, because in this limit the 

phase space for the undetected photon tends to zero. 


^0 1 ^min 


(1 — 3ijB)(l — ^Jb) 

1 —^jb(1 —^jb) 


(4.15) 
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Figure 4.3: Leading logarithmic radiative corrections to the tagged photon 
cross section for = 10 GeV: (a) Jacquet-Blondel method, (b) S method. 
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whereas for small yjB the corrections remain moderate. Since the Jacquet- 
Blondel variables correspond to a “more inclusive” measurement than the 
leptonic variables, the corrections due to radiation of an additional photon 
are generally smaller. 

Finally, hgure |4.3| b displays the corresponding results for the S method. 
Here the corrections are surprisingly small. This apparent suppression is 
easily traced back to the weak dependence on initial state radiation of the 
reconstruction of the kinematic variables using the S method, see table ^ 
At the same time the /^-dependence is negligible 


4.2 Complete leptonic corrections 

As explained in the previous chapter, the subset of leptonic QED corrections 
to deep inelastic scattering is gauge invariant. Assuming one-photon ex¬ 
change it also factorizes, thus allowing a discussion isolated from the hadronic 
part. In the present section, we shall therefore consider the complete set of 
QED corrections to the Compton subprocess 

e(pi)+7*(-?) ^ e(p 2 )+7(fc) , (4-16) 


with the emission angle of the photon being integrated over the PD, while 
taking into account the corrections from virtual and real QED corrections. 
While performing this integration, we require that the remaining part of the 
amplitude of the full process, Ad (7* +p —> X), depends only weakly on 
the transverse recoil momentum of the forward photon. This implies the 
condition pi ■ fc -C pi ■ P 2 , i-e., 'do -C 9. 

The presentation below follows the outline given in pi 
more details. 


but provides 


4.2.1 Compton tensor 


We start our calculation with the Compton tensor as dehned in (|3.7|) . 
For the discussion of the radiative corrections to this tensor we write its 
decomposition up to and including 1-loop contributions as follows 


7^1-loop 

1 


2 

p 

^ flip 

~ PfU 

pfj.l' 

~ 9 iJ.i 




+ PvPj 

*7=1,2 

Pi ■ Q 


a 

Bii H- Ti^ 

^ 27r ^ 


(4.17) 


q 


2 ) PifJ, PifJ. QfJ,' 




i = 1,2 . 
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The expressions for the quantities Bij corresponding to the Born approxima¬ 
tion are given in (|3.12|) . The general results for the one-loop QED contribu¬ 
tions Tg, Tij in the high energy limit can be found in [Q; they are lengthy 
and will not be reproduced here. 

The next step will be the integration over the solid angle of the photon 
that is tagged in the forward detector. Since we are interested in the region of 
almost collinear emission, k ~ {l—z)pi, we neglect the transverse momentum 
of the emitted photon in the tensor decomposition ( [4.17] ), as it is of order 
Oi^Ee'do), and use momentum conservation to set p 2 = zpi. To the same 
accuracy, the kinematic variables of the Compton subprocess are related to 
those of the radiative DIS process via: 


u = -- 




= {pi-k- P 2 Y ~ , s ~ 




4.2.2 Virtual and soft corrections 

The virtual corrections to the Compton tensor are furthermore conveniently 
decomposed into a piece containing the universal infrared singular contri¬ 
butions, which are proportional to the Born contributions, and an infrared 
hnite remainder: 


Tg = pBg 


T' 

9 


Tij = pBij + Tlj , i,j = 1,2. 


(4.18) 


The precise form of the infrared singular expression p depends on the reg¬ 
ularization procedure. The calculation of the one-loop contributions in K 


uses a fictitious photon mass A in the calculation of the virtual corrections 
and obtains: 


P — 4(-^q — 1) In--h31n2;-F—— - , 

m ^ 3 2 


with 


Tq = In 


A=ln2: 


(4.19) 


(4.20) 


The integration of the Compton tensor over the solid angle of the collinear 
photon is described in some detail in appendix jA.lj . We quote here only the 
hnal result: 


TT 


= {-Qfgg^ + Azpigpi^) X 
1 


(4.21) 


1-z 


l + ^p)p(z,Lo)-^T +0{PlC^^) , 


27r' 


a 


2ti 
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where 


T = {A\nz + B)P{z,Lo) + CLo + D , 

A = 2Lq — Lq — 2 ln(l — z) , 

B = - 2Li2(l - 2 ;) - ^ , 

C = - Inz — z, (4.22) 

1 — z 

1 — Gz -4- 

D = --(Li2(l-^) + ln^ln(l-z)) 

1 — z 

- 2zln^(l-z) + -^lnz-—z+l. 

1 — z 3 


Here 


Liofa;) = - 


y 


ln(l - y) 


(4.23) 


is the Spence function or dilogarithm [110|. The single and double logarith¬ 


mic terms in Lq and Lq of expression ( [4.22|) agree with [^. 

The dependence of the virtual corrections on the unphysical parameter A 
is canceled by the contribution from emission of an additional soft photon, 
calculated using the same regularization method. Since soft photon emission 
factorizes, its contribution is proportional to the Born piece B^^ and can 
be combined with the p-piece of the virtual corrections. Requiring that the 
energy fraction of the second (soft) photon in units of the energy of the 
incoming electron does not exceed e, with e -C 1, and adding the contribution 
from soft photon emission to the virtual correction then amounts to the 
replacement of the quantity p in ([4.211) by p, see [Q: 


p = 2 (Lq-1) \n — + ?,LQ + ?,hiz-\r?Y- 


TT 


9 


■ + 2 Li2 


1 + c 


(4.24) 


with 


E' 

Y =-^ and c = cos 9 = cos Y{p, p') (4-25) 

Ee 

being the relative energy of the scattered electron and the cosine of the 
scattering angle in the lab system, respectively. 


4.2.3 Double hard bremsstrahlung: ISR 

Kinematics of the double Compton process 

The other contributions to the radiative corrections involve the emission of 
two hard photons. This part is quite elaborate, as we need to keep the elec- 
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Figure 4.4: Three of the six Feynman diagrams contributing to the double 
Compton process (^4.26|) . The remaining three diagrams are obtained by 
exchanging fci ^ 2 - 


tron mass different from zero even at high energies. Techniques for the effi¬ 
cient calculation of processes with double photon emission have been devel¬ 
oped by the CALKUL collaboration (see |111| and references cited therein). 


It is however not necessary to deal with the full expressions for our purposes, 
since we are only interested in the hard photon corrections to the tagged pho¬ 
ton cross section. We shall therefore consider the double Compton process 
of £g. 1^, 

e(pi) + 7 *(-?) ^ e(p 2 ) + 7 (^ 1 ) + 7 (^ 2 ) • (4.26) 


with the implicit assumption that one of the hnal photons is emitted almost 
collinearly to the incoming electron. Furthermore, we will keep only those 
terms that contribute to integrated cross sections at high energy, i.e., for 
\{pi -P2)\ \q^\ > ruT 

For explicit calculations in the HERA frame the following parameteriza¬ 
tion of the momenta is convenient: 


Pi = (EelPeC^) = Ee ■ {!] Pe^z) , (4.27) 

P 2 = (E';Pgn) = Ee • (l;/3'u) , n = (sin^, 0, cos^) , 

ki = XiEe-{l;ni), Uj = (sin"dj cossin"dj sin 0j, cos'd*) , i = l,2. 


The Xi = E^^i/Ef, are the energy fractions of the photons in units of the 
electron beam energy. 

Instead of the polar angles the photons, "d*, we shall make frequent use of 
the substitution 


(4.28) 


Collinear emission of photon i obviously corresponds to r* —0. 

For the kinematic invariants of the Compton subprocess we introduce the 
notation: 


= 2pi - ki , z[ = 2 p 2 ■ ki , 
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a = 2ki- k2 = (ki + ^ 2 )^ , 

^ = -[{Pi - ki - k2Y - rn^] = zi + Z 2 - a , 

~ [(p2 + + ^ 2 )^ — Tn^] = + ^2 + CT , 

Ql = -iPi-P2)\ 

Ql = -q^ = Q‘i + Zi + Z 2 - z[ - z'^ - a . 


In the kinematic region of interest at least one of the invariants Zi ^2 will be 
small, i.e., O {m?)\ 

Zi = 2p- ki = 2xiEl{l - j3e cos{}i) 

= 2xiEl{l - /3e + 2/3eri) ~ Xi {w? + 4Egri) . (4.29) 

In the last line we neglected terms of relative order 0{m?/E‘l). Note that 
also the invariant A may become small in the double collinear region, as 


a 


2ki ■ k2 = 2 xiX 2EI{1 - ni ■ 122 ) 

2xiX2E‘1 [1 — cos -di cos ^2 — sin di sin ^2 cos(0i — 02)] (4.30) 


AX1X2EI 


T 1 +T 2 - 2tiT2 - 2^/tiT2{1 - ri)(l - T 2 ) cos(0i - ^ 2 ) 


The kinematic restriction to at least one almost collinear photon allows cer¬ 
tain simplihcations of the expressions we encounter. To exploit these sys¬ 
tematically we used the symbolic manipulation program FORM ||112| . 


Double Compton tensor 


For the treatment of two photon emission off unpolarized electrons, we nor¬ 
malize the double Compton tensor analogously to (|3.7|) as p4| p| 

E , (4.31) 

^ ^ spins 


where now Mf:'^ 


is the matrix element of the double Compton process 
(4.2t), with the index /i describing the polarization of the virtual photon. 


Double collinear emission 

When both photons are emitted almost collinearly to the incoming electron, it 
is convenient to express the photon four-momenta occurring in the numerator 
of the expression for the double Compton tensor schematically as: 

ki = XiPi + ki^± . (4.32) 

^Note that the normalization chosen here differs from [Q by a factor of 1/4 in order 
to conveniently absorb factors of 4 in other places below. 


65 









The components of the transverse parts ki^± are typically of the order O (m) 
and limited to be smaller than XiEgi^o; they can be safely neglected. Similar 
simplihcations can be applied for 2 ;', A' and Q\. We will therefore set: 

z- ~ 2xiPi ■p 2 ^XiQ‘^ , 

A' ~ (xi + X2)Q‘f , 

Qh - {^-Xi-X2)QJ, 

ki ~ XiPi , (4.33) 

and by momentum conservation: 

P2-Pi-ki-k2 = {l-xi- X2)pi . 

The expression for the double Compton tensor simplifies tremendously when 
we note that in the double collinear region: 

^ ^ ~ ^ (^ 0 ) • 

Introducing the abbreviations 

ri = 1 - Xi , r2 = 1 - X 2 , z = 1 - Xi - X 2 , (4.34) 


the double Compton tensor takes a rather compact form: 


7^2 —coll 


X 


[-9i^uQ^ + 4 ^ iPii^Piu)] 

1 + 2 :^ 1 ^ ^2 A 1 /rf + zr2^ 

_ X1X2 2:12:2 A 2 ^2:2 2:1/ X1X2 V 

Pr^ + z"^ ^2 _|_ ^2 ('x _ 2;)(r^r2 + 2:) \ 

A \ X2ZI Xizl X 1 X 2 Z 1 Z 2 ) 


r2 + 2:ri \ 

2:2A ) 
(4.35) 




+ 42 : 


m 

A2 



This expression is consistent with 
logarithms were calculated. 


108||, where leading and 


next-to-leading 


The integration over the photon phase space 

We assume that the photon detector measures only the sum of the energies 
even if the two photons hit in different positions. Taking into account the 
symmetry factor 1/2! for the emitted photons, we need to calculate: 

4 1 /* 

^ ^ / dfci dk2 0(xi - e) 0(x2 - e) 6 {xi + X2 - {I - z)) . (4.36) 

PD 
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The infrared cutoff e is necessary to regularize the singularity when either 
photon becomes soft. 

For a general shape of the PD the integrations in ( |4.36|) can only be 
performed numerically. To proceed analytically we require, in addition to 
"do ^ 1, azimuthal symmetry of the PD. The result obtained then also serves 
as a useful cross check of a numerical implementation. The integrations over 


the photons in ( |4.36| ) are therefore understood as 


dki = 


PD 


(dvr) 


Xi dxi ^ f dDj 0('do - '&i) , 


(4.37) 


with no restriction on the azimuthal integration. 

Under these conditions the integrals over the solid angles of the photons 
can be performed completely. A list of the relevant integrals is given in 
appendix |A-^ ; along with a description of their calculation. 

The remaining integration over the relative photon energies however is 
very tedious. Nevertheless, one can decompose the result of ([4.36|) in the 
following form: 


2 r 


+ 42: iPi^Piu)] X — [PiS2(2:) + 

(4.38) 

The previously known leading terms containing double and single logarithms 
Lq are contained in 



1 + ^2 


In 


1-z 1-z 

3 + z^ 


+ (1 + z) lnz-2(l-z) 
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-^0 


6(1 — z) + " z + 4 


{1 + zf 


1-z 


In ■ 


PfW + 2 


6(1 - z) + 


1-z 

3 + z^ 


Inz -2 In e 

2 
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-^0 


(4.39) 


In^ z + 4 


(1 + Z )2 


In ■ 


Lq , 


with the second-order leading-log radiator 



1 - 1 - 2 ;' 
1 - 





— \nz + 



+ -(1 -h 2 ;) In 2 ; - l-f 2 ; 


The nonleading terms are split into an infrared divergent piece that depends 
on Ine, 


dIR— div 
nonlog 


(^) 


8z 

1-z 


In 


e 

1-z ’ 


(4.40) 
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and an infrared finite piece -Pnoniog”^" (^)- ontline of the calcnlation and 
more information on the last piece are given in appendix |A.2| . 

Inspecting ( |4.39| ), ( [4.40| ) we find that the terms depending on the soft- 
photon cntoff parameter e in (|4.38|) do factor nicely, as expected from the 
nsnal soft-photon factorization: 

[-QiiuQf + 4:Z ipif^piu)] X (^—j P{z, Lo) ■ 2(Lo - 1) In - . (4.41) 

4.2.4 Final state collinear radiation 

Consider now the kinematic region where one photon is emitted almost 
collinearly to the incoming electron and the other one close to the ontgo- 
ing electron. We snppose that the angle between the scattered electron and 
the second photon is smaller than some angnlar separation parameter Pq, 
i.e., 

Z(k 2 ,p')<Po, with (4.42) 

Withont loss of generality we therefore set: 


ki = xipi + , k2 = 




1-e 


P2 + ^ 2 ,± , 


(4.43) 


where we again assnme that the indicated transverse parts of the photon 
momenta are small, being typically of the order O (m), c.f. the discnssion of 


Ql 


(4.32). These simplihcations 

lead to: 

A 


XiQ‘f , 

A 


(1 Xi)- 

A' 


Ql 



1-C' 

Ql 




1 ^ 

ki 


XlPl , 

and by momentnm conservation: 

P2 + k2= Pi 

— 

ki 


i-e 


ko ~ 


1-e 

QL 




ZP2 , 


1-r 

p2 ~ (1 -a;i)(l -Opi • 


(4.44) 


Fnrthermore: 


2:1 


q: 


2 :£ o 


Q 


i ZOW) , 
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and the double Compton tensor considerably simplifies to 


^ 27 ,FSR 

flU 


X 

X 


Q1 


1 + (1 - Xi^ 1 

Xi Zi 

i + (i-0^ 1 

^ 4 


+ 4(1 - Xi) ■ (pi^piu) 


-2{l-x, 


m 


J 


2 ■ 

m 

^C) 


(4.45) 


exhibiting a complete factorization of collinear initial and final state radia¬ 
tion, respectively. For the application to the tagged photon process, we have 
to identify xi = 1 — z m the above expressions. 


The integration over the photon phase space 


Integrating the expression ([4.45 ) for the double Compton tensor over the 
solid angle of photon 1 at fixed energy we obtain: 


^ j d!J, A'J.™ _ 

PD 




i-e 


-F4(1 - Xi) • (pi^pi,.) 


X 


Xi 


P(1 -Xi,Lo) X 


l + (l-0^ 1 




4 


0— 
J2 


Z. 


2 J 


(4.46) 


The corresponding integration over the final state photon depends on the de¬ 
tails of the experimental situation, i.e., whether the detector is able to resolve 
a photon collinear to the electron, or whether it just measures the sum of 
their energies. We denominated these cases in the discussion in section [3.2.2 
as exclusive and calorimetric measurement, respectively. 

Let us consider first the exclusive case, when only p 2 of the scattered 
electron is measured but the momentum k 2 of the almost collinear photon is 
missed. The phase space for the second photon can be rewritten as: 


dfco = 


1 E: 


'/ 2 


{Airy 


TT 


(1-0= 


dn' 


2 ) 


(4.47) 




with Eg being the measured energy of the scattered electron and the polar 
angle1^® photon being measured with respect to the scattered electron. 
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Performing the integration over the solid angle leads to: 




i + (i-0^ 1 






(dvr) 


with 


L' = In 


(1-0^ 


m 

*7^ 








+ 21nF, r = — 
Ep 


(L'-l)+? 
E' 


The lower limit of the ^ integration in ( 4.48|) is given by 

e 


'Cmin 


F ’ 


(4.48) 

(4.49) 

(4.50) 


while the upper limit depends on details of the full process and on the kine¬ 
matic reconstruction method. It needs to be discussed in the concrete appli¬ 
cation. 

In the calorimetric case we assume that only the four-momentum of a 
cluster (or electromagnetic jet) can be measured when condition (|4.42|) is 
met: 

VL\ = Poster =P2 + k2- (4.51) 

Expressing the occurrences of p 2 in the numerator of (|4.46| ) in terms of 
the tensor structure in the hrst line becomes independent of 

[-9^.u (QHcai + > (4-52) 

because the calorimetrically measured invariant momentum transfer is: 


(^')cal = -(Pl-Pcal)' = 




Ql 


(4.53) 


1 - ^ 1 - xi 

In the dehnition of measured cross sections only p'^^^ can be used. We there¬ 
fore replace the phase space of the outgoing electron by the cluster using the 
substitution: 


dp2 / dfcs 


dPcai / dfca (1 - 0^ • 


(4.54) 


Combining the Jacobian (1 — ^)^ with the integration over the solid angle of 
the photon yields: 

r ^ ri L n _ !r^2 1 ^2- 

/ dk 2 (1-0^ 


1 + 


(dvr) 


de 


^ A 

i + (i-0^ 

e 


12 


{L" -1) + e 


(4.55) 
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with 


/?2 0 / 2 

L'’= In 


+ 2 In y + 2 In 


EL 




■^cluster 


= L' + 21n(l-0, 


since in the inclusive case the ratio Y is given by: 

pf 

Y _ '^cluster 


(4.56) 


(4.57) 


Inserting the limits for the .^-integration which are now independent of the 
full process, 


Y 


<e<i, 


(4.58) 


and evaluating the r.h.s. of (^4.55|) we hnally obtain: 


1 

(dvr) 2 



I) + 3 


27r2- 

“ 3 “ 


(4.59) 


4.2.5 Semi-collinear emission 

The hnal situation covers the conhguration where one photon is emitted 
almost collinearly to the incoming electron, while the other photon is emitted 
at a large angle with respect to both the incoming and outgoing electron 
directions, i.e., for '&2 > 'do and '&2 > 'dj,. We denote this kinematic domain 
as the semi-collinear one. 

To be consistent with the above treatment of the double collinear emis¬ 
sion, we need to perform the solid angle integration over the collinear photon 
and drop all contributions of the order O(i9o) and 0 (^ 0 "^)• In doing so, a sub¬ 
tlety arises from the propagator denominators A = zi Y Z 2 — cr. For strongly 
ordered angles "di -C 'd 2 , corresponding to zi Z 2 , one may approximate: 

ki ~ Xipi , A ~ (1 — Xi)z 2 + O (zi) (for -C i92) ■ (4.60) 


This leads to a factorized expression not only for the differential cross sec¬ 
tion, but also after integration over the collinear photon, as long as i92 3> 'do. 
However, for i9i ^ '& 2 , which can happen if photon 2 is emitted close to 
the forward cone that is specihed by the solid angle of the PD, approxima¬ 
tion (|4.60|) no longer works well. One obtains contributions from further 
terms of an expansion around the limit 'di —0, which lead to a more com¬ 
plicated, steeper 'd 2 -dependence that spoils the factorization. Fortunately, 
these terms fall off rapidly and essentially contribute only in the small region 
■do ^ 'd 2 < few • "do. Since these extra-terms are important only for this par¬ 
ticular kinematic conhguration of almost double collinear emission we shall 
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denote them the quasi-collinear contribution. The details of this calculation 
are given in appendix 

Assuming that photons in this narrow region outside the PD will not 
be measured, we explicitly integrate these additional terms over the photon 
solid angle and split their contribution schematically as follows: 


TT 


4 / 


dDo 


flU 




e 

TT 


^ dDa 


' 1^2 >' 1^0 


/^2 




{ri{Q^ + Z2)r + {r,Q^ - z’) 




riZ2Z2 


riZ2Z2 


T+^P2P2 


Q 


riZ2Z2 


XiTi 


P{n,Lo) 


1 


+ [-gtiyQl + 4(ri - X2) {pl^,Plu)\ - H{xi, X2) , 

X 1 X 2 


(4.61) 


with ri = 1 — xi- The expression for the function H{xi,X 2 ), which collects 
the non-factorizing, quasi-collinear terms, is given in appendix P, eq. (|B.25|) . 
It is infrared-finite, 'do-independent for small 'do, and does not contain any 
logarithm of a large scale. 

The integrand in the first part on the r.h.s. of ( |4.61| ) can be rewritten as 


- —Piri,Lo) ■ —K^°’^^{ripi,p2,k2) , 


Xi 


n 


(4.62) 


where for the sake of consistency one should drop terms of order on the 
r.h.s. in the expression for the lowest order Compton tensor (|3.7| ). 

Since in our decomposition of phase space only photon 1 reaches the PD, 
we have to identify ri by 2 ; and Xi by 1 — 2 ; in the above expressions. However, 
we still need to integrate over the phase space of the other photon that is 
emitted at large angles. This calculation depends on the complete scattering 
process and in general requires a numerical integration. 


4.3 Corrections to radiative DIS 

Let us turn to the description of the radiative deep inelastic scattering process 
with photon tagging, 

e{p)+p{P) ^ e{p') + X{Px)+jiki) ( + 7(/^2)) , 

including the complete leptonic QED radiative corrections calculated above. 
It is straightforward to contract the expressions for the radiatively corrected 
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Compton tensor with the hadron tensor (p.l9|) . We assume here the kine¬ 
matic reconstruction by measurement of the scattered lepton and later outline 
the changes for other methods. 

Applying the results from the previous section, we hnd for the contribu¬ 
tion from virtual and soft corrections to the cross section: 


1 


a 


y dxdydz 


[P(;^,Lo)p-T]S(i,y,g^), 


(4.63) 


where p is taken from (|4.24|) with 


y = = z{l-y)+xy — 

H/fp il/fp 


= cos 9 = 


z{l - y)Ee - xyEp 
z{l - y)Ee + xijEp 


(4.64) 


In the calculation of the contributions from the emission of two hard 
photons, we decompose the phase space into three regions discussed in the 

4) both hard photons hit the forward photon 


(see also 32 


previous section 

detector, i.e., both are emitted within a narrow cone around the electron 
beam ( 191^2 < 'do, 'do -C 1); (ii) one photon is tagged in the PD, while the 
other is collinear to the outgoing electron {d '2 = Z(fc2,p0 — 'do)j hnally 
(hi) the second photon is emitted at large angles (i.e., outside the dehned 
narrow cones) with respect to both incoming and outgoing electron momenta. 
For the sake of simplicity, we assume that m/E^. '9^ 9. 

The contribution from the kinematic region (i), with both hard photons 
being tagged in the PD, but only the sum of their energies measured, reads: 


1 


dVg? 

y dxdydz 


a 


2 r 




0(2) r,( 2 ),IR—div. I 73,(2),IR—fin. 

^ log ' nonlog ' nonlog 


(4.65) 


see eq. ([4.38|) . 

In region (ii) we need to distinguish between the cases of whether the 
outgoing electron in the presence of a collinear photon can be measured 
separately (exclusively), or whether its energy and momentum are detected 
together with the electron (inclusively), as this affects the reconstructed kine¬ 
matic variables. 

For the exclusive event selection, when only the scattered electron is de¬ 


tected, we obtain from ([4.46|) and ([4.48|) : 


1 a 


y dh dy dz 




de 


( 1 - 0 ^ 








(4.66) 
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where 


S/ = L' = \n^^ + 2\nY , 




Xf 




xy 

y-i ’ 

e 

Y 


yf = ^- 

'Cmax 


y 


Q^f = 




i-c i-e’ 

(1 - x)y - (m2 - m2) /{zS) 
1 - (M 2 _ ]Vf 2 ) /(zS) 


(4.67) 


- (1 • 


In the case of a calorimetric event selection, where only the sum of the 
energies of the outgoing electron and collinear photon is measured and taken 
into account in the determination of the kinematic variables, the correspond¬ 


ing contribution from ([4.46|) and ([4.591) is proportional to the lowest order 
contribution and reads 


1 dVS 


(ii),cal 


y dxdydz 


a 

47r2 

47r2 


P{z,Lo)Jd^ 

^min 




P(z,Lo) 


21n|-^)(P-l) + 3 


(P-l + 21n(l-0)+e 

27r2 


S, (4.68) 


see also [31, 32, 33 


For the calculation of the contribution from the semi-collinear region (iii) 
we apply the decomposition ([4.611). Due to the factorization property ([4.62|) 


of the leading term, it is useful to introduce the “radiation kernel” 


P{zp,p',k2) = —Kp^{zp,p',k2)HP''{P,zp-p'-k2) 

on 


(4.69) 


1 


ZZ 2 Z 2 _ 


G ■ Fi{xh, Ql) + Xh{zSf[l + (1 - yf] 


XhM^^ 

Ql 


G 


+ (z5)[(i - mQl - 4) - (Ql + ^^ 2 )] )f2{xh,QI) 


which can be obtained from eqs. ( |3.19|) , ( |3.20| ). Here: 


Ql — Ql + ^^2 


Xh = 

G = 


Ql 


4 , Z2 = 2p-k2 , 4 = 27»' ■ ^2 , 
Ql 


2P ■ {zp — p' — /C 2 ) yzS — 2P ■ k 2 

Ql + ZZ2^ + (^QI — ^2 j • 


(4.70) 
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The semi-collinear contribution then reads 

>,2 


1 dV- ^ 


y df d^dz TT^ 


l^l Qh 


+ -^ I dx2 —^— H{1-z,X2)T.{xt,yuQt) ■ (4.71) 


47r^ 


Z - X2 


In the integration over the solid angle the cones corresponding to double 
collinear ISR (region (i), half opening angle 790 ) cind FSR (region (ii), half 
opening angle 'dg) have to be excepted. The upper limit on the photon energy 
is obtained from (p.36|) by substituting Ef, —>• zE^, S —>• zS, Y Y/z: 


77 max 
-^2 — 


R[z - R(1 - r)] - AE^zYt - {M^ - M^) 
4 [Ep{l - Ti) + zEeTi - YEeT2] 


(4.72) 


For the second, quasi-collinear part, we have used the abbreviations 


z 


Xt = 


X2)xy 


zy - X 2 


yt = 


zy -X 2 
- X2 


Ql = 


X2 


(4.73) 


and under HERA conditions and for small x the upper integration limit 
simplifies to 

{l-x)y- /{zS) 


x\ = z- 


^ ^ — zy . (4.74) 

1 — xy 

The total contribution from QED radiative corrections is hnally obtained 


by adding up ( |4.63|) , (H:.65|) , ( [4.71| ), and, depending on the chosen event 
selection, ( ^4.66| ) or ( [4.68| ). As has been demonstrated explicitly in the 
unphysical IR regularization parameter e cancels in the sum. 

In the case of a bare electron measurement, the angle 'dg plays the role 
of an (unphysical) regularization parameter that serves to define the phase 
space region of collinear hnal state radiation. It can be shown to cancel 
between the contributions from regions (ii) and (iii) and to drop out of the 
hnal result if chosen small enough, see |]3^ . 

For a calorimetric measurement of the scattered electron, ■dg corresponds 
to a resolution of the detector and is actually physical. Therefore the cross 
section will depend on it. However, in this case the mass singularity that 
is connected with hnal state collinear emission cancels in the sum of the 
contributions from virtual and soft corrections and from regions (ii) and 
(iii), in accordance with the Kinoshita-Lee-Nauenberg theorem |8^. For 


sufficiently small "dg the resulting cross section depends logarithmically on 
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Figure 4.5: Radiative corrections to the tagged photon cross section with 
leading logarithmic and with full leptonic accuracy for the bare electron 
selection and for = 10 GeV. 


i9q, while for a coarse detector, i.e., for ■dg ~ C>(1), the result agrees at the 
leading logarithmic level with the correction obtained in section E3 

Let us now compare the radiative corrections calculated with full leptonic 
accuracy to those obtained by considering only the leading logarithms. For 
the same set of parameters as in section H. hgure displays the radiative 
correction 


f^RC = 


d^cr 

d^ClBorn 


(4.75) 


at leading logarithmic and with full leptonic accuracy for the bare electron 
measurement for x = 0.1 and x = 10“^ and for a tagged energy Ej = 10 GeV. 
We hnd a large difference of the order of 5 to 10 percent between both 
calculations, especially in the region of small x and for large y. 

In the case of the HERA main detectors, the measurement of the energy 
of the scattered electron in the electromagnetic calorimeter is performed with 
a typical separation from an almost collinear photon of the order of (several 
cm)/(l meter) ~ O {50 mrad). Figure [4.6| displays the corresponding com- 
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Figure 4.6: Radiative corrections to the tagged photon cross section with 
leading logarithmic and with full leptonic accuracy for calorimetric event 
selection with "dg = 50 mrad and for = 10 GeV. 


parison of the radiative corrections for the calorimetric selection, assuming 

■dg = 50 mrad . (4.76) 

Apparently the difference between the leading logarithmic and the full lep¬ 
tonic result becomes signihcantly smaller, especially at smaller x and for 
small y, reaching a level of the order of 5 percent. This suggests that the 
leading logarithmic approximation gives a reasonable estimate for the ra¬ 
diative corrections for a sufficiently inclusive measurement, at least for the 
reconstruction of the kinematic variables using the lepton method. 


In the references [^, ^ a calculation of the radiative corrections was 
performed that takes into account the leading and next-to-leading logarithms 
of the full corrections. We can now compare this approximation with the full 
result. It turns out that it is rather close to the full calculation 133, and 


that the non-logarithmic terms, typically contributing at the per mille level, 
are potentially signihcant only at small x and for a large tagged photon 
energy. To demonstrate this we compare in hgure |4.7| the corrections in the 
calorimetric case taking into account the next-to-leading logarithms vs. the 
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Figure 4.7: Comparison of the radiative corrections with next-to-leading log¬ 
arithmic and with full leptonic accuracy for a calorimetric event selection 
with dg = 50 mrad and for = 20 GeV. 


full leptonic corrections for a higher tagged energy of = 20 GeV. Similar 
results are obtained for the bare electron measurement as the non-logarithmic 
terms do not change the contribution from region (ii) and the large-angle part 
of region (iii). 

The size of the contribution from hard photon emission to the radiative 
corrections can be reduced if one applies experimental cuts. As an example 
which is applicable at HERA we dehne a “missing longitudinal energy” using 
longitudinal momentum conservation: 


F/miss Ee — E-pD — 


S/j + E'f,{l + COS 6) P ■ k2 


2Er, 


(£f + pW) , (4,77) 


where Epp, represents the actually measured energy in the photon detector. 
Assuming 100% efficiency of the PD and the second photon being lost in the 
forward beam pipe outside the PD, the quantity E^^ss is roughly equal to the 
energy of this lost photon. We then apply the following cut: 


A F'miss „ _ 
miss • ^ U.O , 


El 


PD 


(4.78) 


78 






Figure 4.8: Comparison of the radiative corrections with leading logarithmic 
and with full leptonic accuracy for a calorimetric event selection without cuts 
and with cut (|4.78|) for = 10 GeV. 


and compare the radiative corrections with and without this cut for the 
calorimetric selection. For the same set of parameters and for x = and 
= 10 GeV the result is displayed in hgure K 


We observe that the size of the corrections is much reduced, and that the 
difference between the leading logarithmic and the full result is also lessened 
slightly. A main reason is that this cut is mostly efficient for the contribution 
from very hard lost photons which would lead to the largest shifts between 
reconstructed and “true” kinematic variables. Therefore this cut is effective 
only for sufficiently large y, see also ( [4.14|) . It also affects the QED Gompton 
contribution which plays a role only at large y. 

For other reconstructions of the kinematic variables the above formulae 
have to be adapted in the contributions from hnal state radiation (ii) and 
for the semi-collinear contribution (iii). Ref. discusses the modihcation 
of the calculation of the next-to-leading logarithms for the S method. The 
derivation of the boundaries of the photon phase space is more elaborate 
than for the electron method, so that we will not treat it here and only quote 
the results. 
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Figure 4.9: Radiative corrections with leading and with next-to-leading log¬ 
arithmic accuracy for the S method for a calorimetric event selection with 
■dg = 50 mrad and for = 10 GeV. 


We found in section [4.1| that the leading logarithmic corrections for the S 
method were quite small, see £g. Hb. The reason for this behavior can be at¬ 
tributed to the scaling properties of the kinematic variables under collinear 
photon emission, see table However, as can be seen from £g. the 
more accurate next-to-leading logarithmic calculation adds terms that are 
the order of a few percent and thus of similar size as the leading logarith¬ 
mic estimate, even for the calorimetric selection. As can be expected the 
corrections become signihcantly larger for a bare electron measurement but 
still stay well below those for the kinematic reconstruction using the lepton 
method. Furthermore, the R dependence of the corrections remains small 
even after taking into account the next-to-leading logarithmic terms. For 
more details we refer the reader to 


33 
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Chapter 5 

Estimates of Higher Order 
QED Corrections 


In this chapter we shall attempt to obtain an all orders estimate of the 
leading logarithmic contributions to the radiative corrections to the deep 
inelastic scattering process with tagged initial state radiation. To this end, 
we will hrst rederive the hnal result of section |4.1| in a simpler framework. 
We then generalize the arguments to all orders to hnd a closed expression 
after resummation. 

Our starting point will be the factorization theorem for the QED ra- 
diatively corrected, photon inclusive cross section for high energy scattering 
processes. 


(^Rc(s) ^ / d2 D(2,^p)a/,(zs;/i]7) . 


(5.1) 


Here s stands for a typical large scale of the process. All large logarithms in 
the small electron mass m are contained in the radiator function 
while the “hard cross section” ah is hnite for m —> 0. Also, the cancellation 
of infrared singularities from QED virtual corrections and from soft photon 
emission takes place within this radiator function at all orders. The factor¬ 
ization theorem (|5.1|) is a corollary of the corresponding theorems for the 
factorization of infrared and mass singularities in QCD |p.01|] , but has also 
been proven explicitly for QED in ||107||. Note that we have suppressed in 


(5.1) the dependence of the hard cross section on the renormalization scale 
which is formally unrelated to /i|^. The purpose of factorization in QCD 
in the above example is to extract all soft (long distance) contributions into 
the distribution D, so that the hard process receives only perturbatively 
calculable short distance corrections. 
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Eq. ( |5.1| ) introduces on the r.h.s. an arbitrary unphysical scale, /i^. Physi¬ 
cal results, and thus the l.h.s. should in principle be independent of the choice 
of this so-called factorization scale. This formal independence can be used to 
derive renormalization group type equations. In practical applications, how¬ 
ever, things are more complicated. Calculating the radiator function using 
a partial resummation of large terms at all orders in the coupling constant 
while treating the hard part at hxed order in perturbation theory leads to an 
unavoidable residual dependence of the r.h.s. on /i|, in higher orders that are 
not completely calculated. On the other hand, one can show that the leading 
dependence of the uncalculated corrections to ah is logarithmic in If all 
small scales can be extracted into the radiator function D, one may hope 
that the missing logarithmic terms in ah will be small if this factorization 
scale is chosen close to a typical large scale. 

In the simplest case, the leading logarithmic approximation, the hard 
cross section is approximated by the so-called “improved Born cross sec¬ 
tion” for the hard scattering process by replacing coupling constants by the 
running couplings taken at the scale For the photon inclusive cross sec¬ 
tion with no explicit registration of the radiated photons there essentially 
is a single large scale, e.g., the total energy s or a momentum transfer . 
The ambiguity in selecting one of these scales reflects the uncertainty of this 
approximation; it can only be resolved by a calculation with higher (next- 
to-leading) accuracy, i.e., taking into account terms beyond the Born cross 
section. Nevertheless, the leading logarithmic result often gives a reasonable 
estimate. The structure function method in the form as applied in section [4.1 
corresponds to this simple approximation. 


5.1 Exclusive photons in the leading logarith¬ 
mic approximation 


The radiator function appearing in ( f).l| ) satishes an Altarelli-Parisi 

type evolution equation. As we shall be interested in the leading logarithmic 
contribution only, we restrict ourselves to the electron non-singlet structure 
function D^^{z, /r^), which is the solution of the integro-differential equation,[| 


1 



2 : 


^For the sake of simplicity we shall ignore the running of the QED coupling. 
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with the electron splitting function 


Pee{x) = \im (Q{1 - e - - 5(1-x) [ ^ d|/| (5.3) 

6^0 I 1-x J 1-y J 

and with initial condition 

m^) = 5(1 — z) . (5.4) 

The normalization of the non-singlet structure function, 

1 

J dz = 1 , (5.5) 

0 

valid for any /i^, is a consequence of 

1 

j dx Pee{x) = 0 , (5.6) 

0 

which expresses the order-by-order cancellation of the infrared singularities 
between virtual corrections (x = 1) and soft photon emission (x —>■ 1). 


The formal solution of the integro-differential equation (^.2|) with initial 
condition (^.4|) reads: 

^^)=Hl-z) + -Xj^ j . (5.7) 


A perturbative solution is obtained from ( ^.7|) by iteration. It may be written 
in different equivalent forms, e.g.. 


D^^{z,fi^) = 6{l-z) + 


^ f dQl 

271 J Ql 


P^Jz) 


(5.8) 


a \ 2 r c1Q2 r ^q 2 j- 


27r/ J Ql J Ql 

m? 


Pee{y)Pee (“)+••• 


= 5(1 - ^) -h — 


a f dQl 


271 J Ql 


Pee{z) 


(5.9) 


1 


/ 1 

a f dQ^ 


271 J Q2 


V 


dy 


Pee{y)Pe. 




y 


+ ... . 
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As long as we are interested in photon inclusive radiative corrections only, 
both forms (|0|) and (|5.9|) are equally useful. However, the hrst variant 
turns out to be powerful for the treatment of exclusive photon emission in 
QED. This can be seen as follows. 

For QED radiative corrections, the origin of the leading large logarithms, 
L = ln(/i^/m^), is identihed by analyzing the peaking behavior of the squared 
matrix elements for photon radiation in the region of phase space where the 
photons are widely separated in polar angle (measured w.r.t. the radiating 
fermion). 

Let us start by considering a typical expression corresponding to single 
photon emission. In that case, the large logarithm results from the integra¬ 
tion over the photon polar angle, which is of the type (r = (1 — cos'd)/2):[| 

f d{p ■ k) C d(cos'd) f dr f dr 

J p-k J 1-(3 cos'd J T + {I - f3)/2 J T + im?/{AEl) 

(5.10) 

In the high-energy limit, one may drop the term m? j (4Eg) in the denominator 
if one cuts off the r-integration at 


.j-min 


2 

m 

4 ^’ 


(5.11) 


corresponding to a hctitious minimum angle 'dmin — m/Ee. The large log¬ 
arithm L of the non-singlet structure function is then mimicked by 
assuming a hctitious maximum^ 


max _ k' 


(5.12) 


It must be emphasized that this clearly is an effective but not a physical 
cutoff. First, at large angles the approximations leading to (|5.10|) will no 
longer be valid. Second, in (|5.12|) obviously is not restricted to be 

smaller than unity for arbitrary /r^. 

The n photon contributions (n > 2) in eq. (^) can be interpreted anal¬ 
ogously. However, we need to take into account the fact that the leading 
contributions are obtained from the integration over conhgurations in phase 


^In the approximation of small angles, one may equally well work with the photon 
rapidity y = ^ ln[(l -|- cos'd)/!! — cosd)] ~ ^ Inr. 

^For a general discussion of inclusive radiative corrections, the choice of r™'" and 
is completely arbitrary, as long as = yd/rrk. We do know, however, the full 

O (a) result, and can match the leading logarithmic result at this order to the full answer, 
confirming the above choice for r™™. 
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space where the photon polar angles are strongly ordered, 


■di -C '^2 -C • ■ ■ -C , (5.13) 

which is eqnivalent to saying that the (negative) virtnalities of the electron 
after radiation, ronghly given by 2p ■ ki, are strongly ordered for comparable 
photon energies. This snggests to replace in (|5.8|) the nested integrations 


by 



(5.14) 


(5.15) 


The method described above corresponds to the introduction of “uninte¬ 
grated splitting functions” for the generation of exclusive photons from QED 
photon showers in Monte Carlo event generators like KRONOS |^, which 
are chosen in such a way that they reproduce the angular distribution for 
single photon emission in the region of small angles as given by the peaking 
behavior of the electron propagator ,0 


e 


2 



l + (l-x,)^ 

Xi 


.^max 



with 


Pe^{x^ 


1 -I- (1 — x)"^ 

X 


0(x 




d(a;) , 


(5.16) 


(5.17) 


and the small parameter e serving as infrared-regulator. 

In the following we shall assume the photon angle distributions ( 5.15 ) for 
an estimate of the radiative corrections to the tagged photon process. 

Although the rewriting of the scale-ordered integrations (|5.14|) in the 
perturbative expansion ( b.8| ) by the angular-ordered integrations (|5.15|) looks 
very suggestive, a remark on the proper cancellation of infrared singularities 
is in order. 


^Although the Monte Carlo event generator KRONOS generates the proper large in¬ 
clusive logarithms, the actual implementation of the minimum and maximum values for 
the polar angles generally leads to emission of exclusive photons at too large polar angles. 
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In fixed order perturbation theory, the Bloch-Nordsieck theorem guaran¬ 
tees that the emission of soft photons is a Poisson process. The probability 
for emitting n photons in the energy interval oj G [tUmin, cUmax] is proportional 
to 


n! [vr \m^J V'^min 


n 


(5.18) 


The obvious divergence for uj^am —> 0 is known to be canceled by taking into 
account all appropriate virtual corrections at the same order of perturbation 
theory. 

From the reasoning further above, we know that in the regime of hard 
photon radiation the emission of photons is ordered in virtualities resp. an¬ 
gles. Applying this ordering condition to the phase space of the photon leads 
to a cancellation of the combinatorial factor 1/n!; the same will happen in 
eq. (|5.18|) with the combination 1/n! ■ ln"'(/i^/m^). 

We shall interpret this outcome in the following way: emission of a single 
(hard) photon results in an off-shell electron. The finite virtuality effectively 
restricts the phase space accessible to further photon emissions to the range of 
larger angles. At the same time, the virtual corrections to the amplitude after 
the Ath photon emission are regularized or “cut off” by the virtuality of the 
electron, Qf, instead of m^. The ordered perturbative expansion (|5.8|) then 
expresses this order-by-order cancellation due to the normalization property 
( |5.ti| ) of the electron splitting function even for our interpretation of ordered 
photon emission. 


5.2 Leading radiative corrections to DIS with 
tagged ISR 

Let US now turn to an alternative derivation of the leading radiative correc¬ 
tions to DIS with tagged initial state radiation using the above reasoning. 


We shall hnd that not only can the result of section |4.1| be reproduced with 
little effort but also generalized. 

As a start, the lowest order cross section for the tagged photon process 
(|3.83|) follows from the first term of the expansion of the electron non-singlet 
structure function. Integrating over the solid angle of the forward photon 
detector, ?? < 'do, we obtain: 


TO 


^ p f ') dVBorn(a^, 

dxdydz 271 J t dxdy 


aLo 


P^^\z)ao{x,y;z) , 


(5.19) 
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with To = (1 — cos'(9o)/2, Lq = ln(ro/r™™), and the abbreviation (c.f. section 

O): 

, . d'^aBorn{x,y;z) 

At order O (a^), corresponding to donble photon emission, we decompose 
the contribntions into two cases. Case (a) corresponds to both photons hit¬ 
ting the photon detector. Assnming that only the snm of their energies is 
measnred, we hnd: 


(^3^(2,a) 

dx dy dz 


^min n Z 

^ . (5.21) 


In case (b), one photon enters the photon detector, while the other is emit¬ 
ted at a larger angle. Again, to logarithmic accuracy we may neglect the 
transverse momentum of both photons, as we did in section 0- Taking the 
angle ordering into account yields: 


(£) / ^ I ^ I (^o{z,y]z,z' 

^min TO 

1 

(£) ^oLiP^^\z) j dz' Pee{z') ao{x,y,z,z') , (5.22) 

z' 

^min 

with Li = ln(r™®'^/ro) = /m?) — Lq. Remember that the semi-collinear 

contribution does depend on the method of reconstruction of the kinematic 
variables. Therefore, do in the above expression is understood to be expressed 
in terms of the “true” kinematic variables of the hard subprocess. 


d3^(2,b) 

dx dy dz 


Cro(x, y] Z, z') = Cro(Xtrue, 2/true; Zz') ■ J (x, y\ Xtrue, 2/true; z!) , (5.23) 


with the Jacobian J depending on the chosen set of kinematic variables, see 
table |0| . The integration limit z' ^^^ also depends on kinematic cuts. 

Inserting the regularized form of the electron splitting function ( ^.3|) in 
the last line of ( |5.22|) leads to: 


dV^ 

dx dy dz 


(^)T„AP<y.) 
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X 


dz' P^^\z') {ao{x,y]z,z') - ao{x,y;z,z' = 1)) 


- ao{x,y;z,z' = 1) / dz'P^^\z') 


(5.24) 


Adding contributions ( b.21| ) and ( b.24| ), identifying z' = 1 — u, and noting 
that 

do(x, y, z, z' = 1) = ao{x, y, z) , 
we recover result (|4.13|) . 

The above procedure may be generalized to higher orders in a; we only 
need to pay attention to the angle ordering to hnd the proper combinatorial 
factors. As a further explicit example we sketch the application at order a^. 

The contribution from three forward (collinear) photons reads: 


(^) ^')Mx,y]z,z') . (5.25) 


1 / aL 


The semi-collinear part receives contributions from two collinear and one 
large-angle photon, 


1 / uLq 


^ Pi\z) j d^' {^^^^P^^\z')aQ{x,y-z,z'), (5.26) 


2! V 27r 

and from one collinear and two large-angle photons: 




27r 


P^^’iz) dz' 


2! V 27r 


P^^\z') 


ao[x,y]z,z 


(5.27) 


It is now straightforward to perform a summation of the leading contributions 
to all orders. Those parts corresponding to collinear emission with no photon 
at large angles obviously add up to 


D {z;Lo)- / dz' 6{1- z')ao{x,y,z,z') , 


(5.28) 


with the suggestive notation D^^{z; Lq) indicating that it coincides with the 
inclusive non-singlet structure function, 

n=0 ^ '' 


(5.29) 













and thus satisfies an evolution equation, 




NS, 

y 


D^^{z;0) =6{l-z) . 


(5.30) 

The slight change in notation here is intended in order to emphasize that the 
present version of has no exclusive meaning in the sense of the discussion 
in the previous subsection; it just represents the formal sum on the r.h.s. of 
(|5.29|) . Furthermore, and even more importantly, while D^^{z;L) appears 
to be well defined, the partial cross section (|5.28|) with only forward photons 


is by itself not a meaningful, physical expression, as it lacks contributions 
from large angle photons. This is tacitly expressed by the ^-function in the 
variable z' that is related to the energy loss outside the solid angle covered 
by the forward photon detector, indicating that it is a virtual correction 


part. Strictly speaking, (|5.29|) has to be multiplied by the probability for no 
emission at large angles. 

The full result for the cross section is therefore obtained by summing 
the contributions from any number of collinear and any number of large- 
angle photon emissions while taking into account the angle ordering. It is 
easily seen that this procedure finally leads to the following expression for 
the all-order leading logarithmic result: 


73 LL(oo) « 

= D^^{z;Lo)- j dz' D^^{z'-L,)d,{x,y-z,z') 


= D^\z-Lo) 


(To{x, y,z) 


\ 


1 - j D^^{z';Li) dz' 
0 


(5.31) 


+ / Li) ((To(a:, !/;2, i') - !/; 2, l))dj' 


In the last step we have used the normalization (|5.5|) of the non-singlet 
structure function and exhibited the cancellation of the singular behavior 
for z' 1. 


Expanding the all-order result (|5.31|) in powers of a confirms that the 
terms up to order O (a^) terms agree with the result of section |0. 
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5.3 The non-singlet structure function 


In the previous section we found a very compact expression for the tagged 
photon cross section including the leading QED radiative corrections re¬ 
summed to all orders. It involves the electron non-singlet structure function 
L), where L stands for any of the large logarithms Lq = ln(Eg'dQ/m^), 
Li = ln(Q^/m^) — Lq. During the calculation of the one-loop radiative cor¬ 
rections to this process in chapter we already encountered the hrst two 
orders of a perturbative expansion in a of We shall now sketch methods 
aimed at practical determinations of this function for numerical applications, 
mostly following Jezabek Przybycieh [p.l4|| , 

these papers. 

Introducing the parameter^ 


and references quoted in 


f3 = —L, 

TT 


(5.32) 


the Gribov-Lipatov evolution equation for the electron non-singlet structure 
function reads: 


dD^^{z,f3) 




(5.33) 


with boundary condition 


D^^{z,0) = 5(1 -^) . 


(5.34) 


The integro-differential equation (|5.33|) can be solved in various ways. 
A perturbative solution is obtained by a power series ansatz, 


P<”I(j) + 0(/3~«) , (5.35) 

n =0 V / 

®In the usual treatment of the Gribov-Lipatov equation, the running of the QED cou¬ 
pling a is incorporated by replacing j3 by 

m2 

This modification serves to take into account the additional corrections from pair radiation. 
However, the forward taggers at HERA clearly detect only real (on-shell) photons. We 
shall therefore disregard these contributions at least for the tagged collinear initial state 
radiation and work with constant a. 
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yielding 


P^°\z) = 6{l-z) , 


(5.36) 


and the recurrence relation 

1 

P,,{y) , n > 1 . 

Z 

From (|5.6|) we find: 


1 

jdzP^^\z)=0, n>l. 
0 


(5.37) 


(5.38) 


The coefficient functions are clearly divergent for z ^ 1 and need regu¬ 
larization. Using a small auxiliary parameter e —> 0 as in (|5.3|) , the first few 
terms obtained from recurrence relation ( |5.37| ) read (see also ||115| , |114| , [Uf 
and references): 


pW| 


z = 


feh) 

feh) 


pfi ^ 




Pi‘>5(l - j) + P|’(z)e(l - z - e) 


= 2 


lim 

e^O L 

1 + Z^ 


1-z ’ 


(5.39) 


Pf^ = 21ne + 


2 ’ 


3 \ 1 

2 ln(l — z) —Inz + - \ + -(1 + z)hiz — 1 + z 
1 z \ 2/2 


21„e + - 


27r^ 


24 ^ ^ ln^(l — 2 ;) -I- 2 ln(l — z) — -\nz ln(l — z) 

1 z y 2 4 2 

1,2 3, 9 TT^ 

-I- \n z -In z H- 

12 8 32 12 


(5.40) 


6(1 -|- z) lnzln(l — z) — 12(1 — z) ln(l — z) + -(5 — 32:) Inz 
3(1 - z) - ^(1 -h z) In^ z + 6(1 -h z) 112(1 - z) , 


pP = 


21„e + - 


27r^ 21n€ + - + 16C(3) . 


As can already be seen from these expressions, the power series expansion 
(|5.35| ) does not converge uniformly. As z —> 1, the coefficients behave as 



n-T 


ln’"-^(l-z) 

1 -z 


(5.41) 
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This is an obvious hint towards a closed form of the solution to the evolution 
equation in the limit z —>■ 1 which was found by Gribov: 


n exp [P/2 ■ (3/4 - 7 )] p 

- T{l + p/2) -• 


(5.42) 


Here 7 = 0.57721... is Euler’s constant. The softening of the perturbative 
behavior of the non-integrable singularity, (1 — z)~^, to an integrable one 
is a consequence of the resummation of multiple soft photon emission to all 
orders (“exponentiation”). 

The knowledge of the solution ( b.42|) has lead to several proposals to 
improve the accuracy of finite order perturbative approximations to 
This was pioneered by Kuraev and Fadin ||105|| who suggested the prescription 


N 


0) = Dc(z, 0) + Y, • ( 5 -«) 


n=l 


where the functions ^n(^) are derived from the requirement that (|5.35|) and 
(|5.43|) are identical up to order P^. 


An alternative prescription was formulated by Jadach and Ward ||117 
They suggested a factorized ansatz, 


D^^(z,0) = Da(z,0)Hz,0) 


— Daiz 


.!})'£{//) Uz) +0(0^*') . ( 5 . 44 ) 


Again, the coefficient functions pn{z) can in principle be determined by com¬ 
paring the expansions (|5.35|) and (|5.44|) . It is however advantageous to first 
derive an evolution equation for the function <l>(z,/3), which can in turn be 
used to obtain a recurrence formula for the pn{z) |p.l3|| . The first three coef¬ 
ficients read: 


= 


1 + z^ 


Pi(^) = -g [2(l-z)2 + (l + 3z^)lnz] , 


(5.45) 


p2(z) = - 


2 I - 4:Z + 3z‘^ l + 7z‘^ 2 

(1 - zY H-In 2: H-— Iffi 2; 


12 


{1-zYUYl-z) 


Further coefficients of the Jadach-Ward series were calculated analytically 
and numerically by Przybycieh ||114|, where it was also found that the higher 
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order terms approach a very regular pattern. It is worthwhile to mention 
that the known coefficients of this series are much more compact than those 
of the perturbative series ( ^.391 ). 

Numerical studies of the solutions to the Gribov-Lipatov equation ob¬ 
tained by Monte Carlo calculation and by inverse Mellin transforms show 
that, if truncated at the same order, the Jadach-Ward series converges sig- 
nihcantly better than the Kuraev-Fadin series and provides a particularly 
good approximation to the full result already with the first three terms for 
typical LEP and HERA energies ||118|, |119|, |115|| except for very small z. How¬ 
ever, even this can be remedied. Inspecting the coefficient functions (|5.45|) 
for z ^ 0 it is evident that the poor convergence in this limit is due to 
0n(^) ~ In” ;2. Studying the z ^ 0 behavior of the Gribov-Lipatov equation, 
this has lead Jezabek to suggest a further improvement of the Jadach-Ward 
series that we will not pursue here. For details we refer the reader to ||113 . 


5.4 Numerical results for higher order cor¬ 
rections 


In order to assess the importance of the higher order QED corrections to the 
tagged photon process, we have calculated the radiative correction factor 


< 5 h.o. — 


/ d^CTh.o. / dVi.o, \ ^ 

\ dx dy dz / dx dy dz J 


(5.46) 


with (Ji.o. being the purely leading order, logarithmic part of the lowest order 
radiative cross section (p.l9|) , and cxh.o. being the radiatively corrected cross 
section taking into account corrections at the level of leading logarithms as 
discussed above. 

We have used the same sets of parameters as in the previous cases and 
the ALLM97 parameterization of the structure function of the proton, with 
hxed R = 0.3. As an example, we have performed the numerical calculation 
for the electron method. Tables |5.1| and compare the numerical results 
for the radiative correction factor Jh.o. for the leading logarithmic correction 
as obtained in section with the all-order resummed result (|5.31|), evalu¬ 


ated using the Jadach-Ward approximation (|5.44|) of the electron non-singlet 
structure function, taking into account the hrst three coefficients (|5.45|) of 
the expansion. 

As can be seen from these tables, the difference between the first-order 
corrections and the all-order resummed result is typically at the level of a 
few per mille. The difference slightly increases for small x and 0, where 
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^h.o. 

a; = 10-^ 

a; = 10-1 


O (a) LL 

all-order LL 

O (a) LL 

all-order LL 

y = 0.05 

0.074 

0.068 

-0.040 

-0.043 

l/ = 0.2 

0.093 

0.089 

0.024 

0.019 

l/ = 0.5 

0.159 

0.156 

0.120 

0.116 

o 

bo 

0.344 

0.344 

0.351 

0.350 

y = 0.95 

0.842 

0.849 

1.412 

1.438 


Table 5.1: Comparison of higher order leading logarithmic radiative correc¬ 
tions at O (a) versus the all-order LL resummed calculation for the electron 
method for a tagged energy of 5 GeV 


(^h.O. 

a; = 10-^ 

a; = 10-1 


O (a) LL 

all-order LL 

O (a) LL 

all-order LL 

y = 0.05 

0.118 

0.113 

-0.016 

-0.019 

y = 0.2 

0.148 

0.146 

0.049 

0.045 

y = 0.5 

0.212 

0.213 

0.147 

0.145 

o 

bo 

0.385 

0.393 

0.379 

0.387 

y = 0.95 

0.681 

0.702 

1.428 

1.486 


Table 5.2: Comparison of higher order leading logarithmic radiative correc¬ 
tions at O (a) versus the all-order LL resummed calculation for the electron 
method for a tagged energy of 20 GeV 
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the phase space for photon emission is strongly constrained and effects of 
multiple soft photon emission become important. 

For ?/ —>■ 1, there is also a signihcant difference between the O {a) and 
the all-order corrections. To understand and to check this hnding, we have 
also calculated the O (a^) leading logarithmic corrections using eqs. (|5.25| ff). 
It turns out that the difference between the O (a^) and the all-order result 
is at the level of only a fewxlO”^ for 0.01 < y < 0.99, establishing that the 
differences in these tables also for large y are truly an indication of the size 
of higher order corrections beyond the first order, being essentially saturated 
by the contributions at O (a^). Note, however, that the effects from higher- 
order corrections are small compared to the uncertainties of the radiative 
corrections due the present errors in the ratio R even for rather small values 
of y. 

Finally, there is still the question of the appropriate scale to be used in 
the leading logarithmic approximation. For the lowest order contribution to 
the radiative cross section, we can determine the “optimal” scale by “match¬ 
ing” the cross section (|5.19|) to the cross section obtained by a fixed order 
calculation, by comparing 


1 + z^ 
1-z 


-LLA 


VS. 


1 + Z^ 

1-Z 


Ln- 


2z 


1 + z^ 


1-z 1-z 


(Lo-l) + l-z. (5.47) 


Obviously, the complete lowest order cross section is better approximated by 
^LLA ^ ^ corresponding to emission of a very hard collinear pho¬ 

ton, while it appears more reasonable to set ~ Lg — 1 for z ^ 1. Taking 
this variation as an indication of the intrinsic uncertainty of the leading loga¬ 
rithmic approximation, we have checked the contributions from higher order 
corrections by varying the collinear logarithm from Lg to Lg — 1. We found a 
change in the contributions from higher order corrections being typically of 
the order of a fewxlO”^ except in the region of large x and y —>■ 1, where it 
can reach the per mille level. However, this remaining scale ambiguity, which 
is resolvable by a calculation of the next-to-leading logarithmic contributions 
at order O(a^), appears to be well below the anticipated statistical accu¬ 
racy of the corresponding experiments even after the upgrade of the HERA 
collider. 
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Chapter 6 

Concluding Remarks and 
Outlook 


In this report we have considered QED radiative processes at HERA, focus¬ 
ing on deep inelastic scattering with an exclusive photon being tagged in a 
forward photon detector. These processes can be used to extend the effec¬ 
tive kinematic range for structure function measurements accessible with the 
HERA experiments down to lower values of the invariant momentum transfer 
Q^. The information obtained this way is important for accurate calculations 
of the radiative corrections to non-radiative deep inelastic scattering. The ex¬ 
tended range allows testing the domain of applicability of perturbative QCD, 
as well as studying the transition into the non-perturbative regime which is 
described by phenomenological models. Furthermore, a measurement of the 
triple differential cross section (]?a/dxdydz, eq. (|3.83|) , enables the separa¬ 
tion of the structure functions of the proton F 2 {x,Q‘^) and Fl{x,Q^) and 
thus a direct measurement of the longitudinal structure function without the 
need to run the HERA collider at different c.m.s. energies. 

The understanding and control of the radiative corrections is crucial for 
precise theoretical predictions of cross sections. We discussed the calculation 
of the most important QED corrections to the tagged photon cross section 
in various approximations. The leading logarithmic approximation provides 
a rough estimate of the radiative corrections. It is quite useful for inclusive 
measurements, when hnal state radiation off the electron is unimportant 
(“calorimetric selection”), and describes the full leptonic corrections with an 
accuracy of several percent. Its compact expressions make it suitable for a 
qualitative understanding of the leading contributions to the corrections for 
different kinematic reconstruction methods. 

The most important result presented in this work is the calculation of the 
full gauge-invariant set of leptonic QED corrections. We find that the differ- 
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ence between the full corrections and the leading logarithmic approximation 
is very signihcant. Depending on the reconstruction method and whether the 
measurement is calorimetric or exclusive, this difference can easily reach 5 or 
even 10 percent. Nevertheless, an approximation that keeps only the leading 
and next-to-leading logarithms appears to work well within an accuracy of 
typically a few per mille. 

Finally, we estimated the size of the QED corrections at higher orders. 
Motivating an exclusive interpretation of the Gribov-Lipatov equation similar 
to QED photon shower algorithms in some Monte Carlo event generators, we 
obtained a compact, closed expression for the all-order leading logarithmic 
contributions to the tagged photon cross section. For this particular process, 
these higher order corrections turn out to be less important than the full 
corrections at relative order 0{a). Furthermore, they are typically smaller 
than those induced by the present uncertainty in the longitudinal structure 
function, F^, at low Q^. 

The apparently large dependence of the radiative corrections on Fl for 
the lepton method is not a real problem. Provided sufficient data is available, 
it can be solved similar to the case of non-radiative DIS by iteration, where 
the result of the analysis is fed back into the calculation of the radiative 
corrections. We thus conclude that the QED radiative corrections to the 
tagged photon process at HERA are well under control with the calculations 
described in this report. What remains to be done is an implementation of 
the QED corrections in a Monte Carlo event generator that is able to perform 
the calculations also for less symmetric and more complicated setups than 
assumed here. 

For an overview of the experimental situation, hgure |6.1| displays the kine¬ 
matic coverage of the {x, Q^) plane as presented at the ICHEP98 conference 
in Vancouver ||120|| . The contributions from HERA to the determination of 
the structure function F 2 essentially he to the left of and above the dashed line 
that roughly describes the acceptance limit of the HERA experiments, while 
those from the hxed-target experiments are to the right and below this line. 
The HERA experiments cover a wide range in y = Q‘^/{xS) at sufficiently 
large x and by combining several measurements, but in the interesting 
region of low x and for < 1 GeV^ the kinematically domain accessible to 
non-radiative DIS becomes very narrow. The large, dark blue area in that 
region, indicated “HI 96 ISR prel.” corresponds to the kinematic domain 
that was added by analyzing events with tagged initial state radiation 


23 


it overlaps with other measurements only at the left and upper boundaries. 
The analyzed radiative events obviously populate only the region of rather 
low Q^, justifying a posteriori our calculation that took into account only 
photon exchange and neglected contributions from the Z boson. 
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Figure 6.1: Measured regions of F 2 in the (x, kinematic plane, as shown 
at the ICHEP98 conference. The nominal acceptance region of the HERA 
measurements corresponds to |/hera > 0.005. The fixed-target experimental 
data occupy the region of high x at low Q^. (Taken from ||120||). 
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Figure 6.2: The structure function F 2 as determined from tagged ISR events 
as a function of x for fixed values of Q^. The data are compared to the 
ALLM97 parameterization (solid curve) and a NLO QCD £t to non-radiative 
F 2 data for > 3.5 GeV^ (dashed curve). (Taken from P5| ). 
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An improved determination of the structure function F 2 from tagged ISR 
events including the HI 97 ISR data has been presented at the DISOl con- 

P3| and is shown in fig. p.‘2\ . The results 
There is good agreement with results 


ference by the HI collaboration 


of this analysis look very promising, 
on F 2 from non-radiative events in the kinematic regions where both meth¬ 
ods overlap. Furthermore, it should be emphasized that the ISR data on 
F 2 describe the important domain where the transition from perturbative 
QCD S> 1 GeV^) to photoproduction 0 GeV^) occurs, within a 

single, coherent data set. This eliminates systematic errors stemming from 
uncertainties in the normalization when combining different data sets. 

In principle it should be possible to simultaneously obtain the ratio R = 
Fl / Ft during this analysis from the ^-dependence of the tagged photon cross 
section (|3.83|) at fixed x and Q^. The extraction method described in |2^ 
has been applied to the HI 1996-97 data set |Q, resulting so far only in 
upper limits on F^ and R which are not yet competitive with more model 
dependent indirect determinations. Improvements in the detectors and more 
data may change the situation. 

In the present paper we have concentrated on DIS with unpolarized 
beams. For polarized protons there are additional spin-dependent struc¬ 
ture functions gi, g 2 - These are being measured at the HERMES experiment 
at HERA and at low using the polarization of the lepton beam and a 
polarized fixed target. 

At the upgraded HERA machine, large longitudinal polarization of the 
lepton beams can also be expected for the HI and ZEUS experiments. A 
primary motivation for lepton polarization is the spin dependence of weak 
interactions which can be exploited in searches for physics beyond the Stan¬ 
dard Model, e.g., extensions to the electroweak sector with heavy vector 
bosons which couple to right-handed currents, or supersymmetry. 

If acceleration of polarized protons were possible in the HERA ring, 
the kinematic domain accessible to spin-dependent lepton-proton scattering 
could be dramatically extended. Gakh et ah ||121|| have considered DIS with 
tagged photons for polarized beams and performed calculations for the corre¬ 
sponding leptonic QED radiative corrections with logarithmic accuracy. The 
leading logarithmic corrections coincide with with those discussed in this 
work, but the next-to-leading logarithms are spin-dependent. The model- 
dependent lepton-quark interference terms which contribute to asymmetries 
such as spin-dependence have not yet been calculated. However, the initial 
plans for proton polarization in the HERA ring were set back with the recent 
HERA upgrade. 

In the description of radiative processes in this report we have neglected 
the contributions which are due to emission off the hadron side; we argued 
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that they are not relevant for the determination of the leptoproduction struc¬ 
ture functions. On the other hand, the case of radiative (quasi-)elastic scat¬ 
tering, 

e + p^e + '^ + p, 

with the photon being well separated from the incident and scattered lep¬ 
ton, has recently received much attention. This process receives contribu¬ 
tions from bremsstrahlung off from the lepton (Bethe-Heitler process, BH) 
as discussed in this report, but also from Deeply Virtual Compton Scattering 
(DVCS) on the proton, 

7*+P^7 + P, 


which gives access to the new class of so-called generalized or non-forward or 
skewed parton distributions (GPDs) ||122|| . Briefly put, GPDs are probability 
amplitudes to knock out a parton from a hadron and to put it back with 
a different longitudinal momentum. DVCS has a reliable theoretical basis 
within perturbative QCD. Its measurement provides a further opportunity 
to gain more insight into nucleon properties such as spin content [|123|| . The 


extraction of DVCS from the HERA data also requires to properly take into 
account the irreducible background from the BH process including the QED 
radiative corrections. 

Radiative processes also play an important role at e+e“ colliders. Most 
prominent are the “radiative return events” at LEP when running at energies 
far beyond the Z resonance: an incoming electron or positron may lose just 
as much energy as needed so that the effective collision energy matches the 
Z mass. Similar radiative return processes are also visible at other colliders, 
e.g., B and (j) factories, and can be used to measure hadron production cross 
sections and resonances below the nominal machine c.m.s. energy ||124|| . The 
high luminosity of these factories can easily compensate the large suppression 
of the radiative process by a factor a/vr despite the absence of a dedicated 
photon detector at small angles. 

Another important application of radiative return events is the improve¬ 
ment of measurements of the ratio 


.^had ("S) 


a{e'^e hadrons) 
cr(e+e“ —>• ’ 


which is needed as input in the determination of the hadronic contribution to 
the running QED coupling a{Q'^). Measurements performed at the DA<hNE 
collider have provided encouraging results ||125|| and have also been used in 
the determination of the hadronic contributions to the muon g — 2 


m 
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Appendix A 

Auxiliary Calculations 


In this appendix we present some details on how to efficiently perform the 
integration over the solid angle of the radiated photons in the forward region 
for the virtual corrections and the double collinear contributions. 

Denoting the energy fraction of the i-th photon in units of the incident 
electron energy, Xi = E^^i/Ee, we write: 


dkj = 


(Att) 


Xi axi — 
TT 


dD,- 


(A.l) 


When a photon is almost collinear to a fermion it is convenient to rewrite 
the scalar product between the fermion and photon momenta as follows: 

= 2p-ki = 2EeE^^i{l - Pe cosdj) = 2EeE^^i{l - /3e + 


1 + (3e 


Xi { m + 


2/5e(l + Pe 


■ AEtTi 


~ Xi (m^ + AEjri) =: Xim'^{l + Ci) 


(A.2) 


In the last line we neglected terms of order 0{im?/E‘l), and we introduced 
the variable 

AE'i 


0 = 




Ti 


(A.3) 


which varies between 0 and AE^/m3. 

Inserting this relation between in the expression for the photon phase 
space we obtain the parameterization: 


dki = 


(47r) 


Xi dxi 


d(j)i 

2tt 


dO 


(A.4) 


103 



A.l Integrals for virtual corrections 

For the calculation of the virtual corrections to the tagged photon process 
we need to integrate the virtually corrected Compton tensor (|4.17|) over the 
phase space of the (single) radiated photon, which is the solid angle of the 
photon detector PD: 



(A.5) 


PD 


We shall assume throughout that the photon detector be azimuthally sym¬ 
metric with respect to the direction of the incoming electron direction and 
covering the range of polar angles 0 << 'do; with -do 9, where 9 rsj 0{1) 
represents the scattering angle of the outgoing electron. 

Introducing 



(A.6) 


m- 


we can treat Lq as a large logarithm for the conditions of the HERA photon 
detectors. It is large in the same sense as Lq = IniQf/m?) S> 1. We therefore 
denote as large logarithmic terms those ones that contain at least one factor 
of Lq or Lq. 

Consistently neglecting of terms of order O ('dg), we may drop terms pro¬ 
portional to the transverse momentum of the radiated photon with respect 
to the incoming lepton. Therefore we take fc ~ (1 — z)p\ and set 


P2 - zpi 


(A.7) 


in the integrand. Among the kinematic invariants of the Compton subpro¬ 
cess, s, t, it and only t is small and may become of order m^. We 
parameterize it as 


i = -2pi ■ k = -m^(l - z){l + C) ; 0 < C < Co • (A.8) 


Furthermore, we set 

s = q‘^ -i-u , u = -Q^ , = -{1- z)Q] . (A.9) 


Inserting the decomposition (4.18) into (4.17), we obtain: 



PD 



PD 


+ PliiPlu {T[i -I- z‘^T22 + ^ {T[2 +T 21 )) 
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The straightforward integration of Tg is essentially elementary and leads to 
transcendent functions including logarithms and dilogarithms. In the con¬ 
vention of Lewin ||11CI||, the dilogarithm (Spence function) reads: 


Lio(x) = 


j — \n{l-y) . 

J y 


(A.ll) 


Neglecting terms of O (Co we hnd: 




TT 


[{A In z + B)P{z, Lo) + CLo + D] =: 


1-z 

A = 2LQ-Lo-2ln(l-z) , 
B = In^ z — 2 Li2(l — z)-, 


Q 1 

l-z 


T 


C = 
D = 


2z 


\nz — z , 


l-z 
1 — 62; -|- 4 : Z ‘^ 
l-z 
2 / 


(A.12) 


(Li 2 (l — z) + ln; 2 ln(l — z)) 


- 2z ln^(l -z) + 


2n^ 


+ 


l-z 


In z -I- 1 . 


The integration over solid angle for the other coefficients of the tensor de¬ 
composition, T'j, i,j = 1,2, is more involved. In that case the evaluation 
of the non-logarithmic terms (in the sense of being free of Lq and Lq in the 
hnal result) unavoidably leads to trilogarithms, dehned as (see ||110|| ): 


Li3(() = Li2((') . 


(A.13) 


The trilogarithms originate from a common integral that can be expressed 
in the following way: 


df 


tQ-\-S 


t t — tr 


Li 2 (l —f) — ln(5) Li 2 (l — fo) + Li3(fo) — Li3(l — fo) 

TT^ 1 

+ Li3(l)—— Into + - In^ to ln(l — to) 

0 2 

+ 0{5) . (A.14) 


This expression, regularized by a small parameter, 5 -C 1, is valid for 0 < 
to < 1- As the integral converges at the upper limit, we have extended the 
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upper limit to infinity. We shall not reproduce here the somewhat lengthy 
calculation which can be done with the help of ||110|| and only requires some 
care in hnding the proper analytic continuations of intermediate expressions. 
The results for the integrations over the photon solid angle read: 


TT 


71 


TT 


E\ 


e 

71 


- (i_^)3L 


-2 {A \nz + B) 


1 + (1-^)' 


1-z 


Lq — (3 — z)ALq 


Az{2 - z) 
1-z 


E + 


3-8z 


z 


C + 


5 - llz + 5z^ 
1-z 


Lq + 2C 


2fl — Az *4- 

+ 2z - 4 - ^ ^ (Li2(l -z) + Inz ln(l - z)) 

+ {3-z) (\n\l -z) + ~ ln(l - ;^) 




-2z^ {A\nz + B) 1±12— 
1 z 


+ (1 — 3z)ALq + 8z(l — z){A\n.z + B) — 4(1 — z^D 

1-4^2+ 8^3 (i-2^)(3-2z)(l-^-^2y 

C H-^- I Lq 


2(1+ 4z-15^2+ 8^3) 


1-z 


1-z 


(Li2(l — z) + \nz ln(l — z)) 


{l-llz + Sz‘^) ln2 (l-z) + 


27r^\ Az{4 — 3z) 


1-z 


E 


2{l-z){l + 2z) 


dr2.v T/n = 


2z{2 - z) 


^ {l-z)^[ 1-z 

+ (3 — z)ALq 


ln(l — z) + 2C + 4 — 6^ 
[(A In 2; + i?) Lq + 2£'] 


(A.15) 


1 — z 


2(3 — 62: + 2;^) .. 

+ 2z -\ -;-(Li2(l — 2;) + In ;2 ln(l — z)) 


1-z 
2 / 


(3 - z) ( ln^(l - ;2) + ^ V ln(l - z) 


dn^T', - 


(Aln^ + B) 


2z^ 

1-z 


Lo - (1 - 3z)ALo -2C + 2z 


l-Az‘^ + 8z^ (l-2z){3-2z){l-z-z^) 

-C H- 


2(1-6z +15^2-8^3) 
1-z 


1-z 


(112(1 — ^) + In ^ ln(l — z)) 
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+ (1-11z + 8z 2) ln2(l-z) + 


27r^ 


+ 


4z(4 - 3z) 


1-z 


E 


2{l-z){l + 2z) 


\n{l - z) - 4{l - zfCL^ 


with the abbreviations A, B, C and D as given above and 


E = Li3(l - z) - Li3(z) + Li3(l) - ^ln(z) 

0 

— ^Li2(l — 2;) + - In 2; ln(l — 2;)^ ln(l — 2;) . (A. 16 ) 

The single and donble logarithmic terms in Lq and Lq of the above expres¬ 
sions (|A-12|) and ( |A-15|) agree with ref. 

It is easy to verify that the above integrals satisfy the relation 


dfl, {4 zT; + Q] [T[, + z^T', + z {T[, + T'^)] } = 0 , (A.17) 


PD 


which, after insertion into ( |A.1(]|) , leads to the factorization of the virtnal 
corrections as in eq. (|4.22|) . 

As a farther check we test the proper behavior in the soft photon limit, 
i.e., for z —>■ 1. The coefficients A, B, C and D are at most logarithmically 
divergent, and we hnd np to terms of order O ((1 — z)^): 


E'i r 
PD 




z^l 


Z 

Ln — 1 


a [ 1 

^ + 2^ (''+2 


Here E^\—Q‘j) is the Dirac form factor of the electron in the one-loop ap¬ 
proximation. The resnlt ([4.22|) thns fnlhlls soft photon factorization as re¬ 


paired, see also eq. (36). 


A.2 Integrals for double collinear emission 

The calcnlation of the contribntion from donble collinear emission in sec¬ 
tion 4.2.3 assnmes that only the snm of the photon energies, {1 — z)Ee, will 
be measnred in the forward photon detector. 
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With the help of (|A.1|) we split the integration of expression ([4.38|) over 
the restricted two-photon phase space in the following way: 


dki d/c2 0(^9o - '^i) - ^ 2 ) 5 ((1 - z) - (xi -f 0 : 2 )) 

1 


(dvr) 


Xi dxi / X 2 dx 2 6 ((1 - z) - (xi 0 : 2 )) 


In the last line we adopted the notation of Arbuzov et ah |p.28 


E 


TT 


dffi df22 ©("do — "di) ©("do — ^ 2 ] 


(A.19) 


(A.20) 


for the angular part of the integrals. 

Equation (|A.2|) shows that in the collinear region the Zi are typically of the 
order of or larger than even in the limit of high energies. This knowledge 
allows us to perform the proper counting of powers in selecting those terms 
in the differential cross section for single and double photon emission that 
contribute in the collinear region. For example, 

a = 2ki ■ k2 = 2 EIxiX 2{1 — ni ■ ^ 2 ) 

= 2 xiX 2 E ‘1 [1 — cos 'di cos '§2 — sin t?! sin 'd 2 cos(0i — 02 )] 

= AX1X2EI T1 + T2- 2tiT2 - 2y/rir2(l - ri)(l - T2) cos( 0 i - ^2) 


= XiX2m 


~ xiX2m 


Ci(l-r2) +C2(l-ri) 


- (1 - ri)(l - T 2 ) cos(0i - 02 ) 

Cl + C2 - 2 V'^cos(0i - 02) 


(A.21) 


In the last line it is assumed that both angles ri _2 <C 1 in the collinear region. 
We then End: 

A = zi + Z 2 - a 


~ m 


xi+X 2 +XiC,ir2-\r X2C2ri + 2xiX2\/C^2Cos{5(t)) , (A.22) 


with ri ^2 = 1 — a^i ,2 and 50 = 02 — 0i. 

In the double collinear limit, the only non-trivial dependence on the az¬ 
imuthal angles of the photons appears in the expression A. Performing the 
integration over the relative azimuthal angles of the two photons yields: 


d(50) 1 
27r A 




{xiCir2 + a;2C2ni Xi -f- X2) - 4a;ia;2CiC2 


-1/2 
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d((5(/)) 1 
2n ^ 


— 

m? 

1 a;iCi ’^2 + 3^2(2’^! + X1 + X2 

£)3/2 


(A.23) 


At this point it turns out to be useful to introduce yet another substitution 
for the polar angle variables of the photons: 

T]i = Xir 2 {l + Cl) , h 2 = a: 2 ?^i(l + C 2 ) • (A.24) 

Inserting these dehnitions into the above expression for D, we obtain: 


D{vi,V 2) = rjl + Tjl + 2r]iri2 


1 — Xi — X 2 — X 1 X 2 4:XiX2 


rir2 


+ 


rir2 


-{rirji + r2T]2) 


= Vi + V 2 + ‘^ViV 2 Cosi> + 2{riT]i + r 2 T] 2 ){l - cosi/j) , (A.25) 


where 


, 1 - xi - X 2 - X1X2 , 2a;ia;2 

cosip := -= 1- 


(A.26) 


rir2 rir2 

Obviously, —1 < cosip < 1 in the physically allowed range 0 < Xj < 1, 
0<Xi+X2<l. Therefore we hnd: 

} d{6(j)) 111 


271 A m2 ^D{r]i,r]2) ’ 

r 

d((50) 1 1 r]i + t]2 + 2xiX2 


271 A2 




A.2.1 Integrals over photon angles 

We shall now provide the relevant angular integrals needed for the contri¬ 
bution of two photons emitted almost collinearly to the incoming electron. 
With the replacement (|A.24|) dehnition (|A.20|) reads: 


Co 

= J dCidC2 

0 


d0i d02 
27r 27r 


m 


XiX2rir2 


a:ir2(l+Co) a;2ri(l+Co) 

/ d77i / d772 


d0i d(j)2 
271 271 


(A.27) 
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and we have: 


(A.28) 


rn^rji w?rj2 

= - , Z2 = - 

r2 ri 

The calculation is performed under the assumption that Co ^ 1 and dg 1 . 

We begin with the integrals that give rise to double and single large 
logarithms Lg. 

l/(ziZ2): 

This one is really trivial: 


1 

2 : 12:2 


1 

X 1 X 2 



(A.29) 


l/(ziA): 

We start with the integration over azimuthal angles: 


/ dCidC2 


d0i d02 1 


27r 27r ziA xiX2ri 


dr]i 

Vi 


dr]2 




(A.30) 


Next, the integral over 772 is easily evaluated: 

/ .d*-' =ln£to.rb). (A.31) 

yD{rii,ri2) 

with 


E{vi,V2) = Vd + 


IdD 

2drj2 


? 7 i cost/’ + 772 + ^"2(1 


— cos Ip) + \/15 . 


(A.32) 


Let us also dehne 

D { r ] i , V 2 ) ■= 77^ + 77^+ 2771772 cost /’, 

Eirii:V2) ■= = T/icosC’+ 772 + , (A.33) 

2 (7772 

which reproduce the leading asymptotic behavior of D and E for 771^2 3> 1. 
Some special values that will be needed below are: 


a/-D(?7i, 772 = a;2ri) 
-^(771,772 = a; 2 ri) 


771 + 0:2(1 + Xi) , 

- ixVi + 2 : 2 ^ 2 ) , 

rir2 


(A.34) 
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where 


D{Vi,V2 = 0) 


E{r]i,r]2 = 0) 


Vi , 

(1 + cos' 0)771 


rir2 


hi , 


\/D{r]i = 0 , 7 / 2 ) = 7/2 , 

^(hi = 0 , 772 ) = 2772 , 


X = I- Xi- X 2 . 


(A.35) 


Before proceeding with the 771 integration, let us note that we are hnally 
interested in dehnite integrals over 772 , with lower limit 77 ™™ = X 2 ri. Choosing 
the integration constant in (|A.31|) appropriately so that the antiderivative 
vanishes if 772 is taken at the lower limit, we can decompose the logarithm in 
the integrand in the following way: 


^(771,772) ^ ^(0,772) _ ^(771,3727-1) ^ ^(771,772) ^ ^(771,772) 

^(771,3:27-1) ^(771,0) ^(771,0) ^(0,772) ^(771,772) 

(A.36) 

Since we shall be interested in the case t/™*™ 3> 1, there will be certain 
simplihcations possible when considering the integral over 771 . 

The hrst integral is readily evaluated. 


r ^ ^( 0 , 772 ) 

J 771 E( 77 i, 0 ) 



2772 

(1 + cos 0)771 


In 771 In 


2772 

1 + cos 0 



(A.37) 


We shall later see that this term contributes to the leading double logarithms. 
The second term in the decomposition ( A.36 ) gives: 


d77i E(77i,a;2ri) 

771 ^( 771 , 0 ) 


d 77 i xVi + 3:27-2 


Vi 


XVi 


d77i 

m 


In 1 + 


3:2 r2 
XVi 


— Li2 

The third pieces gives: 

f ^ ^( 771 , 772 ) 

J 771 ^(0,772) 


V XVi J 


f dr]i 771 cos 0 + 772 + VId 
J 771 (1 +cos 0)771 



(A.38) 


(A.39) 
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where we introduced the auxiliary function 
E{t;x) := 


d'C v^l + 2t^ + + 1 

T "" 2 


(A.40) 


Evaluating the integral for — 1 < f < 1 and using various identities for 
dilogarithms (see e.g., || 11 U||) yields: 




= -In^ 

2 


Li2 

Li2 


Vl + + tx + 1 

2 

(1 + t)x 


a/1 + 2tx + x^ + fx + 1 
(1 — t)x 


(A.41) 


\/l + 2tx + x^ + tx + 1 
The asymptotic behavior of this expression for x —> +cxd is found to be: 

1 


l{t;x) = - In 


2 f {^ + t)x\ vr 


+ — + Li2 

D 


d (A.42) 


We are hnally left with the fourth contribution, 

f dT]i E(? 7 i, 772) _ f dT]i 771 COS'/+ 772 + ^^2(1 - COS'/) + 

J Vi -£( 771 , 772 ) J hi 771 cos / + 772 + vd 

(A.43) 

which appears quite elaborate but can nevertheless be reduced to diloga¬ 
rithms. Fortunately, this effort is not necessary as we need this expression 
for large 772 = X 27 ’i(l -|- ( 0 ) S> 1. As a consequence the ratio ^ = 771/772 does 
not become very large, since we integrate symmetrically over the polar angles 
of the photons. Substituting 771 = ^772 and expanding the integrand for large 
772 , we obtain: 




V2 


(r 2 (l - cos/) O ( 0 ) -h O 


hi/J 


h2 h2 


(A.44) 


Therefore, this expression is (9 (I/ 772 ) and therefore always negligible for sym¬ 
metric integration. 

Collecting the contributions from the first three terms of the decompo¬ 
sition, inserting the integration limits and keeping only terms that are not 
suppressed by factors 1/Co) we find: 


1 

1 

ziA 

XiX2ri 


1 r2 r , ^2rf 

—Lq -|- Lq In-h L 12 

2 xiz 


X2 


XiZ 


+ 


COS'/; 


Xir2 

X2ri 


(A.45) 
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m^/(z^A): 

With the same steps as in the previous case we hrst hnd: 


/ dCidCs 


d0i d(j)2 rn^ _ r2 f diji C dT]2 

271 271 zfA XiX2ri J vf J a//^(?7i,772) 

— TT 

= [ ^ 111 ^( 771 , 772 ). (A.46) 

3 ^ 10 : 2^1 J rjl 

Again we apply the decomposition (A.36|) to the logarithm in the integrand. 
Integrating the sum of the hrst and second term, we obtain: 


d77i 

772 


In 


E{d,7]2) 

E{rii,X2ri) 


d77i rir2772 

^ 5 - In- 

m XVi + X2r2 

X (XVi +X2r2 


(A.47) 


X2r2 


Vi 


1 , 

-In- 

771 xVi + X2r2 


The third term can be neglected for symmetric integration over the photon 
polar angles. This can be seen by substituting 771 = 772 /.^', leading to 


1 

V2 


de' In 


E{r]2/^',r]2) 

E{0,7]2) 


1 

V2 


, ^^ + cos-0+ v/l + 2^Aos7/> + f2 

1 + cos t/j 



(A.48) 


The contribution from the fourth term of decomposition 
be negligible for the same reason. We thus obtain: 


(A.36|) turns out to 



1 

a;ia;27"i 


Lq + In 


a:27"i 

1-z 


z 


In 


1 ^ 1(1 - z) 

XiZ 


(A.49) 


Z2/(ziA2): 

After integration over azimuthal angles we have: 


/ dCidC2 


d0i d(j)2 Z2 


271 271 xiX2ri 


dr]i 


, h 2 (hi+ h2+ 20:1X2) 

77i[D(77i,772)]3/2 


(A.50) 

It is interesting to note that the integral over 772 can be split into one piece 
that is precisely the same as occurred in the discussion of l/( 2 ;iA), plus 
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another one which at this stage is free from logarithms and involves only 
algebraic fnnctions. 



^2(^1 +V 2 + ‘2.X1X2) 

T]i [D (771,772)] 3/2 


A f d772 
Vl J a /£>( 77 i , 772 ) 

1 1 

l + COS 7 /> (771,7/2) 

( 2 -n) ^ (A. 51 ) 

7 ]i ^£>(771,7/2) 

f 1 XV2+X2iXi - x) 

\1 +COS^|J / {xVi+X2r2)VD{rii,ri2) 


The integration of the r.h.s. over 771 now is quite simple. The first term 
leads to an expression that has been discussed further above. The remaining 
antiderivatives read: 


with 

F{Vi,V2) 

G{Vi,V2) 


K{Vi,V2) 


dr]i 


V2 


J Vi VD{Vi,V2) 
XV2 + X2ixi - x) 
ixm + a; 2 r 2 ) a/£>(771,772) 


d77i 

d77i 


In £(771,772) , 
^(^ 1 ,^ 2 ) , 

In £(771,772) , 


(A.52) 


1 d 

\fD + --— = 772 cosip + 771 + ri(l - cosil)) + \/1d , 

2 OTji 


V2 


772 + 2r2(l - cosijj) 


In 771 


— In 


(A.53) 




77 I + 771772 cos ^jJ + {1 - cos Ip) (77iri + 2r]2r2) 
1/772(772 + 2r2(l - cos'ip)] 

Vd — 771 cosip — 772 — r2(l — cosip) (1 — x)(l — cosip) 

XVl + X2r2 XV2 + X2{xi - x) 


Collecting the above expressions and inserting limits, we obtain for the con¬ 
tribution in the doubly collinear region: 


^2 

2:1 A2 


1 

XiX2rf 




XiX 


- 1 + 


a:ia;2 
X . 
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+ 


Li2 

xix 
rir2 - 2x 


\ / , Xir2 

+ ^ cosip; 

X2ri 


(A.54) 


X 


,n(2.,) + !l(ri^ln 


X 


1 - X 


ri(3r2-4x), + 4a;iX , , 

— In X--In ri-— In r 2 


2X 

4x - rir2 

2^ 


2x 2x 

In ij] + X 2 ri + Xir 2 cosip) 


In (r/ + a;ir 2 + 0 : 2^1 cos ip) > , 


rir2 

with the abbreviation 

V = \JD{xir2, X2ri) = \/{xi + X2){xi + X 2 - 4a;ia;2) • (A.55) 


Comparing so far the above results with ||108| , |128|| , we hnd that we reproduce 
the double and single logarithms in Lq. Furthermore, our expressions also 
contain the ‘hnite’ terms for (^o ^ 1 - 

The double Compton tensor in the double collinear region ( |4.35| ) con¬ 
tains also further terms that contribute only non-(logarithmically-)enhanced 
terms. The corresponding integrals are convergent if we remove the upper 
limit by taking C) 0 . Inspecting the expressions for these additional terms 
reveals that the corresponding integrands contribute only when rji < 0{1), 
i.e., the angle of both photons is of the order of '§ 1^2 ^ O {m/EfP). Taking 
one photon at a large angle immediately leads to a strong suppression. We 
therefore take the upper limit of both integration variables rji as inhnite. 


m^/(ziZ2A): 


rn^ / dCidC2 


d^i d(p 2 rnp 


dT]i dr]2 


27r 27r Z 1 Z 2 A X 1 X 2 J r]i r]2 ^/D{T]i,r]2) 

— 71 

One integration, e.g., over 772 , is elementary: 


dr72 


In 772 — In 


rp+r]iri2 cos i/)+(l-cos ■i/))(2?7iri+<72r2) 
a /m [m + 2 r-i ( 1 -cos p)] 


(A.56) 


+ Cl 


t 'jD{’qi,T]2) + 2ri(l - cos^ 

(A.57) 

The integral is convergent for 7/2 —>■ cxd, so we are free to choose the constant 
Cl (w.r.t. 7 / 2 ) in such a way that the r.h.s. of ( A-57|) vanishes in that limit: 


Cl = In 


1 + 


771 cosip + r 2 (l — cosip) 

V^hbhT2r^Tw^cos^ 


(A.58) 
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We therefore obtain for the definite integral: 


dr]2 


In 


/I I 'I I vl+rim [X 2 cos •i/>+ 2 (l-cos i/>)]+ 2 a:ix| 
^ ^ ^ a/)?! [m + 2 »'i (1 -cos i/i)] 


V2 ^/D{rii,r]2) 


^/rii[r]i + 2ri(l - cos 
ln(a;2ri) + ci 


The integration over t]i is attacked with the help of the substitution 
2ri(l —cos-^) drii dpi 


(A.59) 


Vi = 


Pi{pi + 2 ) ’ + 2n(l -cos^/>) ri(l-cos^/>) 


, (A.60) 


leading to: 

00 0 

f dr^i / 

0 

dp2 

r2 


Pi J 

P2JD{PI,P2) ‘2XiX2 

3^ir2 

X2r 

1 


Its evaluation yields: 



m? 

1 1 



_ZiZ2A 


(1 - x) ln(l 


2x1jri 


dpi In 


{pi + l)(pi + 2xi/r2) 


Pi(Pi + 2 x/^2) 


(A.61) 


— Xi Inxi — X2 \n.X2 


(A.62) 


m^/(ziA^): 


/ dCidC2 


d0i d02 


2ti 2ti ziA‘^ XiX2ri 


dpi d?75 


V1+V2 + 2x1X2 

Vi [D{Vi,V2W/^ 


Performing the integration over p 2 , one finds for the r.h.s.: 


(A.63) 


1 

f dr/i 

P 2 - 2 X 2 

XP2 + X2{xi - xY 

2xixl j 

^D{pi,p2) 

. Vi 

XPi + X2r2 


(A.64) 


Obviously, this can be expressed in terms of the already known expressions 
G and IniP given in ( |A.53|) . However, this is not really necessary, as we only 
need the definite integral which is absolutely convergent. Evaluating (|A.64|) 
in the limits p2 = X2ri.. .00, we arrive at a much simpler expression: 


1 

xixl 


dpi 


X2r2 

Pi [xpi + X2r2\ 


(A.65) 
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This immediately leads to 


1 j ri(l — x) 

xixl Xix 


(A. 66 ) 


m^/(z^A^): 

This integral can be dealt with in a similar way. Starting from the expression 


rr/ / dCidC2 


d0i d(l)2 


r2 


2n 2n xiX 2 ri 


dr]i dri2 


V1+V2 + 2x1X2 

ril [D{r]i,r]2)fl‘^ 


1 


(A.67) 

and performing the integration over 772 , we arrive at the following expression 
for the r.h.s.: 


1 

1' drii 

■772 - 2X2 

XV 2 + X 2 {xi - xY 

2xixl j 

7/1 A* ( 771 , 772 ) 

_ hi 

XVi + X2r2 


Again, we evaluate this in the limits 772 = a; 27 "i .. . 00 , to obtain: 


r2 

xixl 



X2r2 

Vibcii + ^^ 2 ^ 2 ] ’ 


The remaining integral is trivial, yielding: 



1 

'1 aiix ri(l - x)' 

_zfAY 

^ 2^2 

X2 XiX 


(A. 68 ) 


(A.69) 


(A.70) 


m^/(ziZ2A^): 

The last and most involved integral among the non-logarithmic ones that is 
needed for the double collinear region is the following: 


/ dCidC2 


d(j)i d(j)2 


2n 2 t\ Z\Z2^ X\X2 


d77i d772 


We start with the integration over 772 to obtain: 


hi + ^2 + 2.XXX2 

771772 [£>( 771 , 772 )] 3/2 ■ 

(A.71) 


1 


a:ia;2 


d77i r r2(ri772 + 43 ;i3;2) ^ rix[xV 2 + X 2 {xi - x)] 
VA)( 771 , 7 / 2 ) L 4a;ia;277f 2xl{xi - x)[xVi + ^^ 2 ^ 2 ] 

^rf + xj-xl riri2[x^ + Xi{Ax - Xi)] 

AxiX2rii Wxlxlrji 


(A.72) 
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rir‘^T]2 


r|(l + X 2 ) 


+ 


dr] 


IGxjxl {xi - x) [r 2 Vi + 4:XiX2] 4:XiX 2 [r 2 r]i + 4:XiX2] _ 

T]i + rir 2 (l — cos'ip) 


^1^2 J + 2 ri(l — cos'ipy 


\nrj2 


In 


rFu -\ , + V 1 V 2 cos'ip + {1- cos 'Ip){2riiri + r] 2 r 2 ) 

VD{vi,V2) +- / i /^i - 

^yr]l[nl + 2ri(l - 008 - 0 )] 


Again, these expressions confirm that the integrals are absolutely conver¬ 
gent. We therefore immediately insert the limits for r] 2 . The remaining rji 
integration is straightforward for the hrst part of ( |A.72| ); the second is part 
is tackled with the help of the substitution (|A.6CI|) : 


riX 


X1X2I 2x1X2 xKx-^i) 

2x21ri 


^^ ri(l-x) ^ rl{l + X2) 


In 


r2 


Xix dxixlpx-Xi) I + X2 


+ 


/v>3 /^3 


dpi 


Pi (2 -h pi)[r2(pi + IY + 2- r2] (pi + 2)(pi 2xiIt2) 


16(pi + 1)2 


6x^x1 |_ 


xl{3 — 2xi) In — ^pyppd. _j_ ^ 2 ^^ 


xix 


2 x 2 ) In 


Pi(Pi + 2x/r2) 

(1 - x)r2 


X 2 X 


In 


A 

rir2 


- 2 x 1 X 2 


(A.73) 


Fortunately, an intermediate spurious singularity in (y —xi) cancels, and we 
hnd: 



1 

_ZiZ2^‘^ _ 

Qx\x\ 


a;2(3 — 2a:i) In 


(1 - x)ri 


xix 


+ xl{3 — 2 x 2 ) In 


In 


A 

rir2 


- 2a;ia;2 


(1 - X)r2 

X 2 X 

(A. 74) 


The remaining integrals can be obtained from those given above by exchang¬ 
ing the labels 1 and 2, i.e., Xi •«-*• X 2 and ri <-»• r 2 - 

Once again we stress that the coefficients of the double and single logarith¬ 
mic terms (Lg and Lq) in eqs. (|A.29|) , ( |A.45|) , ( |A.49|) and ( |A.54|) agree with 
refs. [p.08| , |128|] . Furthermore we have calculated also the non-logarithmic 


contributions. 


A.2.2 Integrals over relative photon energy 

Having performed the angular integration, we still need to integrate over the 
relative photon energy, see (|A.19|). As this integral will be infrared-divergent. 
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we introduce a soft-photon cutoff e on the minimum energy fraction of each 
photon, which is identical to the one used in the soft-photon contribution. 
We thus dehne: 

1 1 


:= / Xi dxi / X 2 dx 2 5 ((1 - z) - {xi -f X 2 )) 


(A.75) 


With the substitution xi ^ {1 — z) ■ u and the abbreviation e = e/(l — z), 
we have, after elimination of the trivial ^-function: 

i-f 

= {1 — z) / du Xi X 2 ■ ■ ■ 

J L J Xl={l-z)u,X 2 ={l-z){l-u),... 

€ 

For the logarithmically (Lq) enhanced leading terms, we find the familiar 
result 


The analytic calculation of the remaining, non-enhanced terms is quite te¬ 
dious, leading to lengthy expressions involving many dilogarithms and trilog¬ 
arithms (see e.g., H). As these terms also contain IR-divergent contribu¬ 
tions, we shall pursue here the following approach. We analytically extract 
those terms in the integrand that either contribute to the infrared-divergence 
as e —0 or survive in the limit z —> 1, before performing the integral over 
the remaining expression numerically. Besides, this separation improves the 
stability of the numerical integration. 

For the IR-divergent pieces of the non-enhanced terms we find 

i-f 


plK—div 

nonlog 


4z 




4z 


du 


8z 


u{l — u) 1 


Ine . (A.76) 


We then decompose the infrared-finite pieces as follows: 

= -Pn^lo, + . (A.77) 

so that P^l^iogiz = 1) = 0. The hrst term on the r.h.s., obtained in the limit 
z —>■ 1 reads: 


pz^l 

nonlog 


(— 

\3a;ia;2 


xi + 3x2 , 3x1 + X 2 , 

In xi -5-In X2 


X 2 


^2 

Xi 




2^2 


X -] X‘ 


1-^2 


X 1 X 2 


1-e 


3 U(l-u) ' (l-u)2 


3 — 2u l + 2u 

+ 7 ^- 7 T In M H-^— ln(l - u) 




16 


(3 + vr^) , 


(A.78) 
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Figure A.l: The infrared-finite, ‘normalized’ part of the double-collinear con¬ 
tribution, P^Zlogiz)- 


where we neglected terms of order e. 


Finally, we plot in figure |A.1| the ‘remainder’ 


ated numerically. 


-Pnoniog(^)) which we evalu- 


120 









Appendix B 

The Semi-Collinear 
Contribution 


In this appendix we shall describe the steps leading to the almost factoriza¬ 
tion (h-61|) of the contributions from double photon emission in the region of 
semi-collinear kinematics. 

We assume that one photon, say photon 1, is almost collinear to the 
incident electron and thus detected in the forward PD, i.e., "di < "do, while 
the other one, photon 2, is radiated at larger angles, {^2 > "do- For the 
computation of this contribution we formally need to perform the phase 
space integration for the two photons. The solid angle part of this integral 
schematically readsiQ 

^ J dn, J (B.l) 

The hadron tensor depends on the photon angles only through the trans¬ 
ferred momentum q = p — p' — ki — k 2 - It is thus natural to expect that the 
tensor is a smooth function of the photon angles, as long as the 

invariant mass of the final hadronic system, = (P + q)^, stays well above 
the inelastic threshold of pion production, S> (M -f . Within the 
desired accuracy, which implies dropping terms of O (^90), we set 7?! to 0 in 
the argument of 

In principle, given this allows us to perform the integrations 

over 01 and ^91 immediately. Before doing this, however, we must hrst de¬ 
termine those terms of the full expression for the double Compton tensor 

^For the discussion below it is useful to temporarily consider the contraction of the 
Compton tensor with the hadron tensor. 
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that lead to nonvanishing contributions in the semi-collinear region. Second, 
we have to hnd a procedure to systematically expand the result in the limit 
of large Cqj which means to hnd the analogue of the logarithmic and non- 
logarithmic terms in the double collinear region. This last step is actually 
non-trivial, which will become clear in the discussion below. 

Let us hrst assume that photon 2 is sufficiently hard, X 2 ~ O {xi) ~ 
0(1), and that it is emitted at a large angle, '§2 ^ "do, so that Z 2 3> zi. 
We hnd that under these conditions the double Compton tensor simplihes 
considerably. Approximating ki ~ xip, one has 1/A ~ l/(riZ 2 ), with ri = 
1 — Xi, etc. Terms of order O ((-2i)°) can be dropped, as they are expected 
not to contribute when integrated over the solid angle of photon 1. It is 
now straightforward to verify that the double Compton tensor in this limit 
factorizes into a collinear radiation factor and the Compton tensor for single 
photon emission, as in the quasi-real electron approximation 


K 


semi—coll 


{p,p',ki,k2)\ 


223>^1 


1 -1- (1 - Xif 

ri L xi 

X K]Z{rip,p ,k2) . 


1 , , 

--2{l-xr) — 
zi zi 


(B.2) 


Strictly speaking, this expression is correct for emission at large angles and 
for radiation collinear to the hnal state electron, while one should drop those 
terms from the Compton tensor in the last line that are of order m? and 
contribute only for small angles, i.e., terms like m? jz\. 

However, if '^2 is small, eq. ([B.2|) not only fails to reproduce the indicated 
terms that we dropped, it completely misses the behavior of the angular 
distribution of photon 2 in that region. To understand this, let us consider 
the integral over the azimuthal angle 0i. Neglecting terms that are always 
suppressed by a factor m ?/relative to the leading ones, we have: 


d(/>i 
2ti a 


{r2Ziy‘ + {riZ2Y + 2rir2ZiZ2 cos^) 


(B.3) 


-1- 4xia;2m^(^i -f Z 2 )) + 


{Zi + Z2)ziZ2 


1 -1/2 


with cos'ijj defined in ( |A.26| ). 

Requiring one photon being almost collinear and the other one at a large 
angle, we see that the terms in the last line of ([B.3| ) can be neglected in 
comparison to those of the hrst line. Thus, for the semi-collinear case we 
obtain: 

/d0i 1 1 


27r A (^ 2 ^ 1)2 + (^ 12 : 2 )^ + 2rir2ZiZ2 cos -i/ 


(B.4) 
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Under the same conditions, we find: 


d0i 1 r2Zi + riZ2 

_ _ _ 

2tt is? [{r2Ziy + {riZ2Y + 2rir2ZiZ2COS'iljY/‘^ 


(B.5) 


It is easy to see that it is not snfficient to keep just the first terms of the 
expansions for large Z 2 when we are going to smaller angles for photon 2. 
Consider e.g., the r.h.s. of ( [B.4| ): 

1 1 r2ZiCOsV’ f z 


+ 0 i . (B.6) 


\/{r2ZiY + {riZ2y + 2rir2^i;22 cos^’ ^ 1^2 (^ 2 : 2 )' 

A similar expansion holds for 1/A^. Retaining only the first term of this 
expansion in the expression for K integrated over leads to the factorized 
form as discussed above. However, when integrating over the polar 

angle 'di, we find that the higher terms in (|B.6|) lead to expressions of the 
type: 

(B,7) 


Eldl 


^2 


^3 


Although these expressions appear to be of formal order C>('do), they do 
contribute when integrating over the polar angle of the second photon, ^ 2 - 
Since they fall off much faster as a function of '^2 than the leading terms (p.2|) , 
their contribution is essentially concentrated in a small region '§2 ^ 'do- We 
express this by formally splitting the double Compton tensor as follows: 


K 


semi—coll 


1 

ri 


1 + (1 - 1 




2:1 


-2(l-xi) 


m 


Kllinp^p', k2) 


+ Rf,„{p,p',ki,k2) , 


(B.8) 


thereby introducing the “remainder” which falls off rapidly as a function 
of '^2- 

Inserting (|B.8D into (|B.1|) yields 


1 

- P(ri,Lo)-^ 

Xiri TT 


dn-2K';i{np,p',k,)H'‘''(p,q) 


q=rip—p'—k 2 






2 


q=(r\—X 2 )p—p' "TT 


dOi d02 R^uip.p', ki, ^ 2 ) • (B.9) 


i 9 i<i 9 o<i 92 


In the second line we exploited again our assumption that the hadron tensor is 
a smooth function of its arguments and the above finding that R^^p, contributes 
to the integral only in a narrow region of the '^2 integration. 
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We have not yet discussed the actual form of From current conser¬ 
vation and the kinematic restrictions it is obvious that its tensor structure 
should be as follows: 


Rtiu{p,p, ki, ^ 2 ) ^ -g^JiuRg + p^lPuRll ■ (B.io) 


Furthermore, we can actually guess the coefficients Rg and -Rn without great 
effort just from the considerations above and from the knowledge of the 
double Compton tensor in the double collinear region, for the simple reason 
being that the “problematic” semi-collinear terms must fall off at least like 
1/^2, as well as contribute in the double collinear case, too. The “remainder” 
Rfj,i, is thus given by: 


Looking at the expression for the Compton tensor in the double collinear 
region (j4.35| ) and requiring that the resulting contribution not be suppressed 
by a factor l/(^o after integration over the solid angle of photon 1, we identify 
the following candidates with a denominator A for closer investigation: 


I I Z2 Zi 

ziA ’ Z2A ’ ziA"^ ’ Z2A‘^ 


(B.12) 


All other terms, e.g., m^/(^^A), will be suppressed and thus harmless. 

Analogous to the double collinear case it is convenient to dehne the fol¬ 
lowing abbreviation for the semi-collinear situation: 




dfli df22 ©('do 


A)e(i>2-A) 


(B.13) 


We shall also use the substitutions described in appendix [A.2| and exploit the 
fact that all integrals below will be convergent even if the upper limit on the 
variable ^2 is removed. 

By straightforward calculation we obtain: 


1 1 
ZiA Zi{riZ2) 




xir^Co 


XiX2ri 

1 

XiX2ri 

1 

XiX2ri 


f d „ 


Vi 


a:2?’lCo 


S (cos'll:, —"j 

xir2Co 

_L V V2j\ 

rii =xir2_ 


\/D{r]i,r]2) 


-I m=^2riCo 


Hi{xi,X2) , 


1 

V2 


(B.14) 
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with the function S being dehned in ( |A. 4 CI| ), and the abbreviation: 

Xir 2 ' 


Hi{xi,X2) = -H cosi/j, 


X 2 ri 


(B, 15 ) 


The second candidate of ( [B.12| ) does not contribute to the naive factorization 
(|B. 2 |) at large angles. Therefore no ‘subtraction’ is necessary, and we have: 

xir 2 Co 


i} 


XiX 2 r 2 

1 

XiX 2 r 2 

1 


X2>'lCo 


<^ri 2 


Vi 


xir 2 


V^iVuV 2 ) 


In r]i In t]2 
H 2 {Xi,X 2 ) , 


/ ^2 
cosip, — 

Vi 


^1^2 Co 


rii =xir 2 


XiX 2 r 2 

with (x = 1 — — 0:2): 

H, = + ^ + 

2 X2X 6 

Similarly, the third candidate yields: 

xir 2 Co 

2:2 


- »? 2 =X 2 riCo 

(B. 16 ) 


X1X2 


I—- 

\ziA2 zi{riZ 2 y 


XiX 2 ri 


X 

hi 


— S cos ijj, 


X 2 ri 

Xir 2 


(B. 17 ) 


dT ]2 


xir 2 


XiX 2 rl 

-2 -^3) 

XiX 2 ri 


X 2 riCo 
hi 


(hi + h2)h2 _ J 

V^(hl,h 2 )' 

xir 2 Co 


COS i/j,— + T(cos' 0 ;hnh 2 ) 

h2 


r?l=xir -2 

(B.I8) 


where 


1 + 2 cosV', h2 + hi cosi/j + a/T>(i7i,?72) 

W(cos' 0 ;hnh 2 ) = -r---In- 

1 + cos ip rji 

\n[r]i + r]2 cosijj + A/D(h^ 7 ^] 


1 + cos ijj 


(B. 19 ) 


One hnds: 


H. = H, 


2 cos ifj 


1 + cos Ip 
1 + 2 cos Ip 


111(20:2+) 


1 + cos Ip 
1 

1 + cos Ip 


1 + cos Ip 
ln(i7 + X2ri + a:ir2 cosip) 


In 


1 + cos Ip 


\n{ri + o:ir2 + X2ri cosip) . 


(B. 20 ) 
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For the last candidate, the same arguments apply as to the second, and we 
obtain: 


I— 

\z2A^ 


oo xir 2 (o 

^ ^ / dm 

XiX2r^ J ri2 J 

a;2»'lCo ^1^2 


(hi +V 2 )Vl 

.VD(Vl:V2y 


Inr/i lnr72 + 2 ( cos'i/', — ) + 'h(cos'0; h 2 ,hi) 

hi 


with 


Ha = 


XiX2r2 
-2 -^4 

a:ia:2i’2 


2 COS '0 , n 

H 2 + :r-^-r li4[(l + cos'ip)x2ri] + 


3^1F2Co 

??l=a;ir2 

(B. 21 ) 


1 + cos 'tp 


1 + cos p) 


In 


1 + cos ^p 


1 


1 + cos ^p 
1 + 2 cos ip 
1 + cos tp 


ln(?7 + X2ri + Xir2 cos -0) 
ln(r 7 + Xir 2 + X 2 ri cos'ip) . 


(B.22) 


As a representative of the harmless terms, we shall investigate whether 
the contribution of m‘^/{z‘fA) at large angles of the second photon deviates 
from the naive approximation. We hnd: 


I 


xir2(o 


r2 


xiX2ri 


dhi f , 
-+ / dh2 


hf 


xir2 

00 


a^2’^iCo 


\^D{m~y2) 


1’2 


XiX2ri 


X2riCo 


dh2 

hi 


h 2 - V'^(hnh 2 ) 


hi 


— cos "0 In 


h2 + hi cos A) (771,772) 


hi 


1 

h 2 


(B.23) 


3:i?’2Co 


?7l=xir2 


One can see even without explicit evaluation of the last integral that both 
terms in the integrand actually contribute only at order 0{1/<Pq) and are 
thus negligible. We leave it as an exercise to the reader to verify that the 
contributions of all other terms beyond the factorization ( [B.2| ) in the semi- 
collinear limit are suppressed by a power of 1/Co for each power of m? in 
the numerator, except to those given in (|B.12|) which after integration are 
leading to Hi.. .Ha. 
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To summarize, the integral over the non-factorizing contribution is found 
to be: 


Et 


vr^ 


dhli dhl2 Riiv = [-giivQf + 4(ri - X 2 ) {p^lPv)\ - H{xi, X 2 ) , 

X 1 X 2 


with (x = 1 — Xi — 0 : 2 ): 


(B.24) 


H{xi,X2) 


rf + Xr2 ^ , 

--ni + 

XiX2ri 


rj + Xri 

XiX2r2 


H 2 -X 



(B.25) 


This quasi-collinear contribution corresponding to one lost photon, i.e., being 
emitted outside the PD, is to be treated in collinear kinematics and only 
depends on the energies of the photons. Obviously, this can only be true as 
long as the strong hierarchy m/Ee -C "do 6* exists, so that the indicated 
approximations remain valid. 

As a check on the quasi-collinear contribution, let us finally investigate 
the leading behavior for X 2 —> 0. We find: 


Hi ~ 

H 2 ~ 


2 X2ri 6 


+ O {X2) , 



Xl 

X2ri 


TT 


+ - + 0 ( 0 : 2 ) 
0 




(B.26) 


leading to a cancellation of the apparent 1 / 0:2 singularity in H{xi,X 2 ) and 
leaving an expression that is integrable for X 2 —> 0. The entire soft behavior 
for photon 2 is thus contained in the factorizing piece (|B.2|). 
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